Naval  Research  Laboratory 

Stennis  Space  Center,  MS  39529-5004 


NRL/MR/7322-95-  7584 


Proceedings  of  the  ONR/NRL  Workshop 
on  Modeling  the  Dispersion  of  Nuclear 

Contaminants  in  the  Arctic  Seas 

October  18-19,  1994  —  Monterey,  California 


Ruth  H.  Preller 

Ocean  Dynamics  and  Prediction  Branch 
Oceanography  Division 


99io»\  Y 

All  JSSIC 

@1.41  1^  te  Midi  -* 


Robert  Edson 

Office  of  Naval  Research 
Arlington,  VA  22217-5000 


OCEAN 
NUCLEAR 
WASTE  SITES 

O  LIQUID 
★  REACTOR 
□  SOLID 


February  22,  1995 


1 9950404  057 


Approved  for  public  release;  distribution  is  unlimited. 


THIS  DOCUMENT  IS  BEST 
QUALITY  AVAILABLE.  THE  COPY 
FUENISHED  TO  DTTC  CONTAINED 
A  SIGNIFICANT  NUMBER  OF 
COLOR  PAGES  WHICH  DO  NOT 
REPRODUCE  LEGIBLY  ON  BLACK 
AND  WHITE  MICROFICHE. 


REPORT  DOCUMENTATION  PAGE 

Form  Approved 

OBM  No.  0704-0188 

Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources, 
gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  or  any  other  aspect  of  this  collection 
of  information,  including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  information  Operations  and  Reports,  1215  Jefferson  Davis 
Highway,  Suite  1204,  Arlington,  VA  22202-4302,  and  to  the  Office  of  Management  and  Budget,  Paperwork  Reduction  Project  (0704-0188),  Washington,  DC  20503. 

1.  AGENCY  USE  ONLY  (Leave  blank) 

2.  REPORT  DATE 

February  22,  1995 

3.  REPORT  TYPE  AND  DATES  COVERED 

Final 

4.  TITLE  AND  SUBTITLE 

Proceedings  of  the  ONR/NRL  Workshop  on  Modeling  the  Dispersion  of  Nuclear 
Contaminants  in  the  Arctic  Seas  October  1 8-1 9,  1 994  —  Monterey,  California 

5.  FUNDING  NUMBERS 

Job  Order  No.  573571 504 

Program  Element  No.  060371 1 H 

Project  No. 

Task  No.  DNA 

Accession  No.  DNI 54-21 7 

6.  AUTHOR(S) 

Ruth  H.  Preller  and  Robert  Edson* 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Naval  Research  Laboratory 

Oceanography  Division 

Stennis  Space  Center,  MS  39529-5004 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

NRL/MR/7322-95-7584 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

Defense  Nuclear  Agency 

6801  Telegraph  Road 

Alexandria,  VA  22310-3398 

10.  SPONSORING/MONITORING 

AGENCY  REPORT  NUMBER 

11.  SUPPLEMENTARY  NOTES 
*Office  of  Naval  Research 

800  N.  Quincy  St. 

Arlington,  VA  2221 7-5000 

12a.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited. 

12b.  DISTRIBUTION  CODE 

13.  ABSTRACT  (Maximum  200  words) 

A  workshop  focusing  on  modeling  the  dispersion  of  radionuclides  which  have  been  dumped  into  the  Arctic  Ocean 
was  held  in  Monterey,  CA,  in  October  1994.  This  workshop  was  sponsored  by  the  Office  of  Naval  Research  and 
hosted  by  the  Naval  Research  Laboratory.  Over  40  participants  attended  this  meeting  during  which  23  oral  presentations 
were  given. 

The  focus  of  this  workshop  was  on  the  existing  modeling  effort  within  the  Arctic  Nuclear  Waste  Assessment 
program,  complementary  data  sets  and  the  future  direction  of  these  efforts.  The  goal  of  the  workshop  was  to  foster 
communication  through  presentations  of  ongoing  work.  The  intent  was  that  discussions,  brought  about  by  these 
presentations  would  help  the  most  important  processes,  effects,  and  issues  that  should  be  addressed  by  present  and 
future  modeling  efforts. 

tea 

'  ’ . *  '  '' 

14.  SUBJECT  TERMS 

Kara  Sea,  nuclear  waste,  nuclear  contamination 

15.  NUMBER  OF  PAGES 

415 

16.  PRICE  CODE 

17.  SECURITY  CLASSIFICATION 

OF  REPORT 

Unclassified 

18.  SECURITY  CLASSIFICATION 

OF  THIS  PAGE 

Unclassified 

19.  SECURITY  CLASSIFICATION 

OF  ABSTRACT 

Unclassified 

20.  LIMITATION  OF  ABSTRACT 

SAR 

NSN  7540-01 -280-5500 


Standard  Form  298  (Rev,  2-89) 
Prescribed  by  ANSI  Std.  Z39-18 
298-102 


Proceedings  of  the  ONR/NRL 

Workshop  on  Modeling  the  Dispersion  of  Nuclear  Contaminants 

in  the  Arctic  Seas 

October  18-19,  1994,  Monterey  California 


Introduction  1 

Workshop  Agenda  3 

Abstracts/Figures 

Data  Presentations  7 

Modeling  Presentations  197 

Workshop  Summary  and  Conclusions  403 

List  of  Attendees  409 


Editors  of  Workshop  Proceedings 

Ruth  H.  Preller  and  Robert  Edson 

Workshop  Support  Staff 

Charlene  Parker,  Thomas  Kozo 


Introduction 


Ruth  H.  Preller 

From  October  18-19,  1994,  the  Office  of  Naval  Research  (ONR)  sponsored  a 
workshop  on  Modeling  the  Dispersion  of  Nuclear  Contaminants  in  the  Arctic.  This  meeting 
was  hosted  by  the  Naval  Research  Laboratory  and  held  at  the  Naval  Research  Laboratory 
Facilities  in  Monterey,  California.  Over  40  participants  attended  the  workshop,  the  majority 
of  them,  scientists  working  on  the  ONR  Arctic  Nuclear  Waste  Assessment  Program 
(ANWAP). 

The  focus  of  this  modeling  workshop  was  to  be  on  the  existing  modeling  efforts, 
complementary  data  sets  and  the  future  direction  of  these  efforts  within  the  ANWAP  program. 
The  goal  of  the  workshop  was  to  foster  communication  between  the  modelers  and  the 
observationalists  through  presentations  of  ongoing  work.  The  intent  was  that  discussions, 
brought  about  by  these  presentations,  would  help  determine  the  most  important  processes, 
effects  and  issues  that  should  be  addressed  by  present  and  future  modeling  efforts. 

Twenty  three  presentations  were  given  during  this  workshop.  Ten  of  these 
presentations  dealt  with  the  gathering  and  analysis  of  data.  Thirteen  presentations  dealt  with 
various  issues  of  modeling  processes  important  to  the  study  of  Arctic  Nuclear  Waste 
Assessment. 

The  workshop  was  organized  in  such  a  way  that  the  data  presentations  were  given 
first,  taking  up  the  majority  of  the  first  day.  Modeling  presentations  began  in  the  afternoon  of 
the  first  day  and  continued  into  the  afternoon  of  the  second  day.  The  data  presentations  were 
organized  such  that  an  overview  of  the  Naval  Research  Laboratory’s  development  of  a 
Geographical  Information  System  (GIS)  data  base,  containing  a  number  of  different  sources  of 
published  data,  was  given  first.  Following  that  presentation  was  an  overview  of  the  1993  and 
1994  Geochemical  and  Environmental  Research  Group  (GERG)  field  experiments  in  the  Kara 
and  Pechora  Seas  and  their  adjacent  rivers.  The  four  following  presentations  dealt  with 
observations  of  sea  ice  and  role  sea  ice  might  play  in  the  transport  of  nuclear  contaminants 
via  sediment  uptake.  The  seventh  presentation  was  centered  on  river  data  for  the  Ob  and 
Yenisei  rivers.  The  next  presentation  was  a  discussion  of  the  circulation  of  the  Arctic  Seas 
based  on  Russian  data  sources.  The  ninth  presentation  discussed  seabed  and  fluid  flow 
conditions  from  observations  in  the  Barents  Sea.  The  last  data  presentation  was  given  on  the 
second  day  and  discussed  the  use  of  AVHRR  data  in  studying  river  outflow  and  ice 
characteristics  in  the  regions  near  river  mouths. 

The  modeling  papers  were  organized  in  such  a  way  as  to  go  from  larger  to  smaller 
scale  models.  The  first  three  presentations  discussed  large  scale  modeling  efforts  of  the 
Arctic  and  its  adjacent  seas.  The  next  four  presentations  looked  at  basin  modeling  efforts  in 
regions  such  as  the  Kara  and  Laptev  Seas.  The  eighth  presentation  discussed  possible 
methods  for  modeling  the  release  of  radioactive  contaminants  into  the  environment  from  the 
Icebreaker  Lenin.  The  two  following  papers  discussed  river  and  estuarine  modeling  efforts. 
The  next  paper  switched  focus  from  ocean  modeling  to  atmospheric  modeling  and  discussed 
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the  development  of  a  regional  atmospheric  model  for  the  Kara  Sea  region.  The  final  two 
papers  discussed  modeling  the  Sea  of  Okhotsk,  another  source  site  for  nuclear  contamination 
and  the  modeling  of  the  motion  of  dense  bottom  currents  which  could  affect  the  transport  of 
contaminants  trapped  in  bottom  sediments. 

The  abstracts  and  figures  presented  here  follow  the  order  of  presentation  at  the 
workshop.  The  text  of  these  proceedings  will  contain  the  following 

1 .  Introduction 

2.  Meeting  Agenda 

3.  Abstracts  and  Figures 

4.  Summary  and  Conclusions 

5.  List  of  Attendees 
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WORKSHOP  ON  MODELING  THE  DISPERSION  OF  NUCLEAR 
CONTAMINANTS  IN  THE  ARCTIC  SEAS 


Hosted  by  the  Naval  Research  Laboratory 

Sponsored  by  the  Office  of  Naval  Research 

October  18-19,  1994 

Naval  Research  Laboratory 
Monterey,  CA 


Modeling  Workshop  Agenda 
Tuesday,  October  18, 1994 

8:30  Opening  remarks  -  Edson/Preller 

Data  Presentations 

9:00  Kathy  Crane  -  Data  Base  Needs  and  Efforts 

9:20  James  Brooks  -  Overview  of  1993  and  1994  Expeditions  in  the 

Kara  Sea  and  Pechora  Seas  and  Adjacent  Rivers 

9:40  Terry  Tucker  -  Radionuclide  Contaminants  in  Central  Arctic  Sea  Ice 

10:00  Coffee  Break 

10:20  Roger  Colony  -  Patterns  of  Monthly  Wind  and  Ice  Motion  in  the 

Arctic  Basin 

10:40  Ignatius  Rigor  -  Russian  Historical  Ice  Charts  and  Dataset 

1967-1990 


11:00  Erk  Reimnitz  -  Transpolar  Drift  -  The  Garbage  Scow  from  Siberia? 
11:30  Lunch 
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1:30  Peter  Becker  -  Data  on  the  Ob  River  System,  Yenisei  River  System 

and  Kara  Sea  Region 

1:50  Andrey  Proshutinsky  -  Circulation  of  the  Arctic  Seas  Based  on 

Russian  Sources  of  Information 

2:10  Richard  Sternberg  -  Time  Series  Observations  of  Seabed,  Fluid,  and 

Flow  Conditions  in  the  Barents  and  Norwegian 
Seas 

2:30  Short  wrap  up  discussion  on  data 

2:50  Break 

Modeling  Presentations 

3:10  Ruth  Preller  -  Overview  of  the  NRL  Large  Scale  Modeling  Effort 

Studying  the  Dispersion  of  Radioactive 
Contaminants  in  the  Arctic 

3:30  Abe  Cheng  -  Model  Results  and  Data  Comparison  for  the  Kara  and 

Barents  Seas 

3:50  Rick  Allard  -  The  Effect  of  River  Discharge  in  a  Large  Scale  Ice/ 

Ocean  Model 

4:10  Ignatius  Rigor  -  Sea  Ice  Transport  of  Pollutants  in  the  Laptev  Sea 

4:30  Meeting  concluded  for  the  day 

Wednesday,  October  19, 1994 

8:30  Max  Coon  -  An  Oriented  Sea  Ice  Dynamics  Model  with  Contami¬ 
nation  Incorporation,  Transport,  and  Release 

8:50  David  Smith  -  A  Numerical  Study  of  Radionuclide  Dispersal  in 

the  Kara  Sea 

9:10  David  Brooks  -  Multiple-River  Interactions  in  Buoyancy-Driven 

Coastal  Currents 
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9:30 

9:50 

10:10 

10:30 

10:50 

11:10 

11:30 

1:30 

4:30 


Mark  Mount  -  Modeling  the  Release  to  the  Environment  in  the 

Kara  Sea  from  Radioactive  Waste  in  the  Dumped 
Reactor  Compartment  of  the  Icebreaker  Lenin 


Break 

Terry  Paluszkiewicz  -  Overview  of  Modeling  from  Land  Based 

Sources  Through  the  Ob  and  Yenisei 
Rivers  and  Estuaries  to  the  Kara  Sea 

Lyle  Hibler  -  Preliminary  Assessment  of  Radioactive  Contaminant 
Transports  from  Catastrophic  Release  Scenarios 
in  the  Ob  and  Yenisei  River  Systems 

William  Thompson  -  The  NRL  NORAPS  Model 

Robert  Fett  -  Use  of  Multi-channel  NOAA  AVHRR  Data  in 
Studying  River  Outflow  Effects  and  Ice 
Thickness  Patterns 

Jiayan  Yang  -  Modelling  the  Sea  of  Okhotsk 

Lin  Jiang  -  Topographic  Effects  on  Dense  Bottom  Currents 

Lunch 

Modeling  discussions 

Group  discussion  of  the  processes,  effects  and  issues  important 
to  the  study  of  the  dispersion  of  radioactive  contaminants  in 
the  Arctic  seas.  What  are  we  doing,  how  can  we  collaborate 
and  what  else  do  we  need  to  be  doing  to  address  the  issue  correctly. 

Meeting  adjourns 
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Data  Base  Needs  and  Efforts 


Kathy  Crane 
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DEVELOPMENT  OF  CIS  DATA  BASES  TO  DETERMINE  THE  STATE  OF  THE 
RADIONUCLIDE  CONTAMINATION  IN  THE  ARCTIC 

Kathleen  Crane^  and  Clare  Brown^ 

^Marine  Physics 

Naval  Research  Laboratory 

Washington  DC  20375 

kathyc0hp8c.nrl.navy.mil,  202-767-0522  (phone),  202-767-0167  (fax) 
^cbrown@boudin. nrlssc.navy.mil,  601-688-4422  (phone) 

LONG  TERM  GOALS 

Until  recently,  the  Arctic  has  been  thought  of  as  remote  and 
pristine,  far  from  the  environmental  problems  associated  with 
industrial  and  agricultural  development  of  lower  latitudes. 
The  Cold  War  cloaked  many  activities  in  the  region  under  a 
curtain  of  secrecy  and  for  most  of  the  world,  the  Arctic 
remained  largely  out  of  sight  and  out  of  mind. 

Information  released  in  1992  on  deliberate  dumping  of  nuclear 
materials  in  shallow  Arctic  Seas  raised  even  more  disquiet 
about  pollution  of  the  Arctic  marine  environment.  However, 
adequate  data  sets  for  the  distribution  of  anthropogenic 
radionuclides  in  the  Arctic  Ocean  and  its  surroundings  did 
not  yet  exist  because  of  the  paucity  of  data  particularly  in 
the  western  Arctic  Ocean  near  North  America,  in  the  Central 
Arctic  Ocean,  and  north  of  Siberia. 

Given  the  presently  unstable  situation  in  Russia  (with  regard 
to  the  state  of  radioactive  dumping  in  the  seas,  the  sorry 
condition  of  its  nuclear  reactors  and  the  uncertain  future  of 
its  stockpile  of  weaponry) ,  it  became  necessary  to 

investigate  all  avenues  by  which  airborne  radioactive  fallout 
and  water  transported  waste  may  reach  the  shores  of  North 
America.  Given  the  lack  of  adequate  radionuclide  data,  it  is 
necessary  to  observe  how  winds,  currents  and  ice  move  and  to 
what  degree  they  may  transport  radionuclide  pollution  from 
the  former  Soviet  Union  across  the  Arctic  Ocean  to  North 

America . 

Existing  time  series  data  needed  to  be  collected  and  compiled 
from  Russian  monitoring  programs  and  integrated  into  western 
data  sets  including  information  on  concentrations  of 

radionuclides  in  the  water,  sediment,  ice,  snow,  flora  and 
fauna  and  how  these  concentrations  have  varied  over  space  and 
time.  It  was  also  important  to  gather  data  from  the 

neighboring  seas  to  provide  us  with  a  means  by  which  we  can 
measure  the  degree  of  radionuclide  pollution  in  the  Arctic 
relative  to  the  rest  of  the  world's  population  centers. 

The  Naval  Research  Laboratory  was  tasked  with  the  function  of 
setting  up  a  GIS  to  address  the  Arctic  radionuclide 


contamination  issue  and  to  create  a  data  base  of  information 
collected  from  ONR  funded  expeditions  in  1993  and  1994 

OBJECTIVES 

The  goals  of  the  1994  NRL  radionuclide  GIS  are: 

1)  to  create  a  data  base  of  existing  radionuclide  data  in  the 
water,  sediment,  ice  and  biota; 

2)  to  create  data  bases  of  bathymetry,  rivers,  sedimentation, 
biota  and  physical  and  chemical  oceanographic,  riverine  and 
estuarine  processes  and 

3)  to  provide  the  platform  for  prediction  and  determination 
of  the  degree  of  risk  these  radionuclides  pose  to  the  Arctic 
environment  and  its  inhabitants  and  to  those  who  depend  upon 
the  Arctic  marine  life  for  sustenance  (this  includes  many 
people  in  the  mid  latitudes) . 


APPROACH 

Our  approach  has  been  multifaceted,  incorporating  a 
compilation  of  preexisting  radionuclide  data,  the 
digitization  of  preexisting  bathymetric,  sediment  and 
physical  and  chemical  oceanographic  data,  the  development  of 
connections  with  Russian  colleagues  to  further  the 
compilation  of  preexisting  data  and  the  development  of 
collection  efforts  for  new  radionuclide,  sediment  and 
bathymetric  data.  In  addition,  some  of  our  efforts  have  been 
directed  towards  developing  interfaces  with  the  ARC/INFO 
system  to  enable  individuals  to  query  the  data  bases  and  to 
gather  statistical  information  related  to  the  distribution  of 
and  correlation's  between  the  parameters  entered.  Efforts 
have  also  been  placed  in  developing  a  more  inexpensive  and 
user  friendly  GIS  operating  system  with  pertinent  subsets  of 
the  data  bases  upon  which  individuals  may  perform  their  own 
analyses  of  the  information  provided. 


STATUS 

As  of  October,  1994,  data  bases  have  been  constructed  for: 

1)  the  distribution  of  nuclear  events  (e.g.  bomb  tests, 
reactor  accidents,  etc.)  which  have  occurred  around  the 
Arctic,  and  the  location  of  radionuclide  dump  sites 
in  the  Arctic,  Atlantic  and  Pacific  Oceans  (figures 
1,2) .  Parameters  entered  include:  ID#,  country,  site. 
Start  Year  of  Dumping,  End  Year  of  Dumping,  Package?, 
Nature  of  Waste,  #  of  Reactors  with/without  fuel,  #  of 
containers.  Container  Type,  Container  Matrix,  Total 
Weight,  (Tons),  Total  Volume  (Liters),  Total  GBq,  Total 


Alpha  GBq,  Total  Beta-gamma  GBq,  Total  H-3  GBq,  Total 
Ra-226  GBq,  Comments,  References 

2)  the  location  of  stations  and.  ship  tracks  of  the  1993 
and  1994  expeditions  in  the  Arctic  Seas  (figures  3,4), 

3)  the  distribution  and  concentration  of  radionuclides 

in  the  water  column  in  the  Arctic  and  its  neighboring 
seas  from  the  surface  to  the  seafloor  and  from  1950  to 
the  present.  These  include,  ^^'^Cs,  ^®Sr,  226j^a,  Pu, 

239/240pu^  241;^n\^  99tc  and  ^29];  (figures  5-9)  . 

4)  the  distribution  and  concentration  of  radionuclides 

in  marine,  lacustrine  and  riverine  sediments  in  the 
Arctic  and  in  its  neighboring  regions (data  are  still 
being  added) .  These  isotopes  include,  127qs,  ^O^o, 

90sr  228Th,  232Th,  234Th,  238u,  55Fe,  226Ra,  210pb, 

212pf|^  214pf)^  40k,  238pu^  239/240pu,  241^^  1440e^ 

95Nb,  952n,  ^^Tc,  and  "^Be  (figures  10-13)  . 

Other  parameters  entered  into  the  sediment-radionuclide 
data  base  include:  ID#,  Cruise,  Station,  Year, 

Day /Month,  Bottom  Depth,  (m) ,  depth  in  core  (cm), 
sediment  type,  %sand,  %silt,  %clay,  %carbonate,  -51120, 
guain  size.  Porosity,  Sedimentation  rate,  macrofauna 
/m2 ,  biomass  (mg/m2) ,  year  of  deposition,  comments  and 
references . 

5)  the  distribution  of  nuclear  power  plants,  weapons 
factories,  and  labs,  sites  of  plutonium  and  uranium 
production  and  enrichment,  nuclear  test  sites  and 
military  sites  around  and  in  the  Arctic. 

6)  digitized  Arctic  bathymetry  from  500  m  to  ocean 
depths 

7)  digitized  Arctic  Rivers,  and 

8)  digitized  distribution  of  important  fish  stocks, 
marine  mammals,  birds  (used  as  food  sources  in  the 
Arctic) 

In  addition,  efforts  are  being  taken  to  gather  up-to-date 
bathymetry  shallower  than  500  m  along  the  Russian  Shelves 
with  the  assistance  of  Russian  oceanographers.  In  addition, 
sediment  type,  physical  properties,  physical  and  chemical 
oceanography  for  the  Kara  Sea  and  the  Arctic  Ocean  are  being 
entered  into  the  GIS. 

Further  development  of  ARC/INFO  interfaces  have  progressed  as 
well  as  the  development  of  a  MAPINFO  GIS  with  subsets  of  the 
data  (designed  for  use  by  scientists  who  may  not  have  the 
resources  to  purchase  or  access  the  full  ARC/INFO  system) . 


PROPOSED  DEVELOPMENT 


In  addition,  we  are  planning  to  add  to  the  GIS  the  following 
data  bases,  depending  on  levels  of  funding: 

BIOTA 

1)  distribution  and  contamination  of  benthic  biota; 

2)  distribution,  concentrations,  seasonal  variation  and 
contamination  of  fish  in  the  Arctic  and  near-Arctic 

3)  distribution  and  contamination  of  other  important 
food  sources  such  as  seals,  walrus,  beluga,  narwhal, 

4)  distribution  of  indigenous  and  other  Arctic  people; 

5)  statistics  on  food  sources  of  the  peoples  of  the 
Arctic  and  dose  assessments  where  available. 

MARINE  GEOLOGY 

1)  seafloor  sediment  type; 

2)  seafloor  sedimentary  thickness 

3)  recent  sedimentation  rates 

4)  detailed  bathymetry  of  the  Arctic  Seas; 

5)  detailed  acoustic  backscatter  information; 

6)  characterization  of  the  transition  from  the  riverine 
environment  to  the  oceanic  environment . 

PHYSICAL  OCEANOGRAPHY 

1)  compilations  of  temperature,  salinity,  density,  and 
oxygen . 

2)  seasonal  oceanographic  and  riverine  information 
pertinent  to  the  spring  thaw  and  flushing  of  the  Arctic 
rivers (including  salt  wedge  and  flocculation  zones  in 
estuaries) ;  rare  events  (such  as  unusually  large 
riverine  fluxes) 

3)  "remote  sensing"  compilations  of  ice  movement  and 
transport  on  a  daily,  monthly  and  yearly  basis, 

IMPACT  FOR  APPLICATIONS  AND  RELATIONSHIP  TO  OTHER 
PROJECTS 

The  NRL  GIS  is  an  integral  part  of  an  international  attempt 
to  provide  base  line  information  to  oceanographic  modeling 
and  monitoring  efforts.  Data  from  the  GIS  have  been  used  by 
the  Norwegian  Radiation  Protection  Agency,  Tokai  University 
in  Japan,  the  University  of  Edinburgh,  the  Norsk 
Polarinstitutt ,  KEMA,  Netherlands,  the  IAEA,  Vienna,  the  IAEA 
Marine  Environmental  Lab  in  Monaco,  and  KORDAI,  from  the 
South  Korean  Institute  of  Oceanography,  and  the  German 
Hydrographic  Service  under  the  auspices  of  the  IAEA  and  the 
AMAP  program. 
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Overview  of  1993  and  1993  Expeditions  in  the  Kara  and  Pechora 

Seas  and  Adjacent  Rivers. 

James  M.  Brooks,  Michael  A.  Champ  and  M.  Baskaran 

Presented  at  Workshop  on 
"Modeling  of  Nuclear  Contaminants  in  Arctic  Seas" 

18-19  October  1994 
Monterey,  California 

Abstract 

GERG  conducted  two  extended  sampling  expeditions  in  1993  and  1994  in 
the  Russian  Arctic.  The  first  yecU-  expedition  was  conducted  in  August  and 
September  1993  using  the  M/V  OLFSON  from  the  Hydrobase  in  Igarka  on 
the  Yenisey  River.  This  cruise  successfully  sampled  ~80  stations  in  the  Ob 
and  Yenisey  Rivers  and  Kara  Sea.  Figure  1  shows  the  station  locations 
sampled  during  the  1993  expedition. 

The  1994  expedition  was  conducted  using  the  M/V  YAKOV  SMIRNITSKIY 
out  of  Arkhangelsk,  Russia  on  the  White  Sea.  The  vessel  left  Arkhangelsk 
on  10  July  with  a  scientific  party  of  18  (10  Americans  and  8  Russian 
participants).  Five  graduate  students  participated  in  the  cruise  collecting 
samples  for  M.S.  and  Ph.D.  degrees.  The  Russian  participants  included 
groups  from  the  Zoological  Institute  of  St.  Petersburg  and  the  Arctic  and 
Antarctic  Research  Institute.  The  M/V  SMIRNITSKIY  expedition 
successfully  sampled  stations  in  the  Pechora  and  Kara  Seas,  and  Ob,  Yenisey 
and  Pechora  Rivers  (Figures  2  and  3)  during  July  and  August  1994.  A  total  of 
88  stations  were  occupied  during  the  2-month  expedition.  Approximately 
500  miles  of  geophysical  data  were  collected  using  the  Chirp  subbottom 
profiler.  Fifteen  (15)  cores  were  collected  and  logged  with  the  multi-sensor 
core  logger.  Biological  tissue  samples  were  collected  from  37  sites.  Permits 
for  these  cruises  were  provided  by  Mr.  Danilov-Danilian  Minister  of  the 
Protection  of  the  Environment  and  Natural  Resources  and  closely 
coordinated  with  the  Research  Institute  for  Nature  Conservation  of  the 
Arctic  and  North  (RINCAN). 

The  M/V  SMIRNITSKIY  was  equipped  to  provide  full  oceanographic 
sampling  capabilities  for: 

•  Chemical  Contaminant  Studies  (full  sampling  capabilities  for  bottom  and 
suspended  sediment,  biota  and  water); 

•  Physical  Oceanographic  Studies  (hydrographic  and  CTD  instrumentation); 

•  Biological  Studies  (trawls,  box  cores  and  traps);  and 

•  Geological  Studies  (side  scan  sonar,  subbottom  profiler,  piston  cores  and 
box  cores). 

Sampling  efforts  consisted  of  sampling  water,  surficial  sediment  and 
organisms  for  radionuclide  analysis.  All  sediment  and  tissue  samples  were 


placed  in  clean  containers  and  returned  to  the  US  for  analyses.  Surfieial 
sediments  (upper  2-cms)  were  collected  with  a  small  box  corer.  Benthic 
sentinel-indicator  organisms,  including  filter  feeder  bivalve  molluscs  were 
collected.  Selected  water  samples  were  collected  and  processed  through 
special  extraction  columns  to  remove  Cs  and  Pu.  Samples  of  particulates 
and  colloidal  material  were  collected  for  radionuclide  analysis  to  determine 
transport  phases  and  partitioning  of  isotopes.  One  of  the  major  routes 
through  which  transportation  of  many  of  the  reactor-derived  radionuclides 
in  the  coastal  waters  take  place  are  through  the  association  with  the 
colloidal  material.  Rates  of  deposition  and  transport  will  be  determined. 
The  samplings  were  closely  coordinated  with  the  sedlmentological  and 
geological  component  of  this  study  that  identified  the  depositional  areas 
(i.e.,  troughs,  fine  grained  sediment  depositional  areas,  etc.)  of  the  marginal 
seas. 

GERG's  1993  expedition  to  the  Ob  and  Yenisey  Rivers  and  Kara  Sea  provided 
some  of  the  first  samples  taken  by  USA  scientists  of  the  sediments,  waters 
and  biota  from  this  Arctic  region.  One  of  the  most  interesting  findings  of 
the  initial  survey  was  the  diseoveiy  of  a  depositional  lobe  of  radionuclides  in 
sediments  in  the  Yenisey  River  below  the  salt  wedge  (Figure  4).  Before  our 
survey,  the  general  expeetatlon  was  that  the  Ob  River  would  have  higher 
levels  of  radionu elide  contamination  in  its  estuarine  waters. 

An  initial  evaluation  of  plutonium  coneentrations  in  samples  from  the  1993 
expedition  to  quantitatively  evaluate  the  contribution  of  Pu  from  the  nuele^ 
reactors  dumped  in  the  Kara  Sea  has  been  conducted  and  published  in 
Applied  Radiation  and  Isotopes.  The  distribution  of  Pu  concentrations  and 
inventories  in  surfieial  sediments  (upper  3  cm)  depends  on  many  factors 
besides  local  inputs,  including  sedimentation  rate,  organic -matter  content, 
variations  in  miner^ogy,  grain  size  parameters,  and  water  content  of  the 
sediments.  The  239,240pu  concentration  in  the  surfieial  sediments  of  the 
rivers  and  the  Kara  Sea  varied  between  9.4  and  627  mBq  kg  ^,  with  a  mean 
of  250  mBq  kg  U  In  the  five  sediment  eores  collected  from  the  Ob  and 
Yenisey  Rivers  and  Kara  Sea,  the  maximum  depths  at  which  239,240pu  was 
detected  varied  between  6  and  12  cm.  Assuming  that  the  Pu  was  introdueed 
to  the  sedimentary  layer  in  the  year  1952  (Koide  et  ah,  1975),  the 
penetration  depths  of  239,240pu  can  be  utilized  to  obtain  the  average 
apparent  sedimentation  rates.  These  estimates  of  sedimentation  rates  vary 
between  1.5  and  3.0  mm  yr'k  The  238pu  concentration  varied  between 
below  detection  limit  and  26.5  mBq  kg'l,  with  a  mean  of  8.5±6.1  mBq  kg  U 
The  distribution  of  239,240pu  concentrations  in  the  surfieial  sediments  is 
shown  in  Figure  5.  There  is  no  systematic  trend,  either  with  distance  from 
the  coast  line  or  Avith  water  depth,  on  the  239,240pu  concentrations. 

The  Pu  inventory  obtained  from  the  ocean  water  column  as  well  as  soil 
samples  from  the  Arctic  region  are  given  in  Table  1.  Since  a  major  portion 
of  Pu  in  the  world  ocean  still  remains  in  the  water  column  and  the 
scavenging  of  particle-reactive  radionuclides  is  poor  in  the  open  oeean 
Aretlc  waters  where  the  particle  concentrations  are  low,  most  of  the  Pu  is 


likely  to  remain  in  the  water  column.  The  Pu  inventory  for  five  sediment 
cores  from  the  Ob  and  Yenisey  Rivers  and  the  Kara  Sea  varied  between 
2.69±0.45  Bq  m'^  and  24.42±2.21  Bq  m'2,  with  a  mean  of  8.86  Bq  m‘2.  In 
the  three  sediment  cores  from  Kara  Sea,  the  Pu  inventories  varied  between 
2.6110.63  and  10.911.3  Bq  m'2.  Among  the  two  sediment  cores  collected  at 
adjacent  sites  in  Kara  Sea  (water  depth  =  39  m)  within  a  distance  of  less 
than  1  km,  Pu  inventories  varied  by  a  factor  of  4.2  (2.61  and  10.9  Bq  kg  l). 
These  values  are  comparable  to  the  soil  inventory  values  of  239,240pu^  14  8  Bq 
m*2,  measured  at  Barrow  (71.6°  N),  Alaska  (Hardy  et  al.,  1973).  Also,  these 
values  are  considerably  lower  than  the  water  column  inventory  of  239,240pu^ 
24.2  and  35  Bq  m-2,  near  the  North  Pole  (89°  N)  reported  by  Livingston  et 
al.  (1984).  If  all  the  Pu  from  the  dumped  naval  reactors  (i.e.,  2-5  x  10i3  Bq) 
would  be  released  to  the  sediments  uniformly  over  the  entire  Kara  Sea,  then 
that  could  result  in  an  additional  Pu  inventory  of  19.6-48.3  Bq  m‘2.  From 
the  average  vsdue  of  the  measured  inventory  of  Pu  (8.5  ±6.1  Bq  m'2)  in  the 
five  sediment  core  samples  collected  from  the  Ob  and  Yenisey  Rivers  and 
Kara  Sea,  and  taking  the  area  of  Kara  Sea  to  be  9.95  x  10^  km^,  the  total 
inventory  of  Pu  in  the  Kara  Sea  is  estimated  to  be  (8.5  +  6. 1)  x  lO^^  Bq. 

The  different  sources  of  radionuclides;  such  as  global  fallout  due  to  nuclear 
weapons  testing,  nuclear  reprocessing  plant  effluents  such  as  from 
Sellafield,  leachates  from  dumped  nuclear  reactors  in  the  seabed,  etc.;  have 
distinct  238pu/239,240pu  activity  ratios.  The  published  data  on  the 
238pu/239.240pu  activity  ratio,  decay-corrected  for  the  year  1993,  for  the 
Arctic  region  is  given  in  Table  2.  The  238pu/239,240pu  activity  ratio  in  the 
surficial  sediments  varied  between  0.009  and  0.065  with  a  mean  value  of 
0.034.  If  most  of  this  Pu  is  derived  from  the  global  fallout  of  nuclear  tests, 
then  during  the  past  30  years,  the  238pu  has  undergone  radioactive  decay,  by 
about  27%.  Thus,  variations  in  the  sedimentation  rate  could  affect  the 
activity  ratio  in  the  surficial  sediments.  The  best-fit-line  between  238pu  and 
239,240pu  concentration  in  these  samples  yields  an  activity  ratio  of 
0.034±0.003.  Comparing  this  value  with  the  published  values  on  the 
European  nuclear  effluents  discharged  into  the  coastal  waters,  feillout  values 
of  the  nuclear  weapons  tests,  and  the  estimated  Pu  inventories  in  the 
dumped  reactors  in  Kara  Sea,  we  conclude  that  there  is  virtually  no 
detectable  input  from  either  the  European  nuclear  effluents  or  the  dumped 
nuclear  reactors  in  the  Kara  Sea. 


Figure  1.  1993  Kara  Sea  sampling  stations. 
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Figure  2.  1994  Kara  Sea  sampling  stations. 
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Figure  4.  distributions  from  1993  Kara  Sea  sampling  stations. 
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Figure  5.  239,24opu  concentrations  in  surficial  sediments  at  1993  sampling  stations. 


Table  1.  239, 240pu  Inventories  in  the  Water  Samples  and  Sediments  from  the  Arctic  Regions. 


Location  and  Nature 
of  Sample 

Water 

Depth 

(m) 

Year  of 
Collection 

239,240pu  (Bq  m-2)* 

Source 

76.6°  N;  Thule,  Greenland 
Soil  Sample 

- 

1971 

12.2 

Hardy,  1973 

71.3°  N;  Barrow,  Alaska 
Soil  Sample 

- 

1971 

14.8 

Hardy,  1973 

64.8°  N,  Fairbanks,  Alaska 
Soil  Sample 

- 

1971 

31.5 

Hardy,  1973 

65°55.7  N;  27°27.0  W 
Water  Samples 

646 

1972 

20 

Livingston,  1985  Data 
reported  inAarkrog,  1989 

74°56.2N;  r7.2  W 

Water  Samples 

3740 

1972 

54 

Livingston,  1985  Data 
reported  inAarkrog,  1989 

89°  N;  89°  W 

Water  Samples 

2497 

1979 

35.0 

Livingston  etal,  Pr  i4 

89°  N;  111°W 

Water  Samples 

3000 

1979 

24.2 

Livingston  et  al,  1984 

89°  N 

(Expected  from  Fallout) 

- 

- 

5.13 

Livingston  et  al,  1984 

68°  -  76°  N;  67°  -  84°  E 
Sediment  Cores  from  Ob  and 
Yenisey  Rivers  and  Kara  Sea 

7-290 

1993 

8.85  (2.61  -  24.4) 

this  study 

♦ 


Numbers  in  parenthesis  denote  the  range  of  values  reported. 


Table  2.  238pu/239,240pu  Activity  Ratios  in  Nuclear  Effluents  from  Sellafield  and  Cap  de  la  Hague,  Dissolved,  Particulate 
Phases,  Surficial  Sediments  and  Terrestrial  Samples  in  the  Arctic  Region. 


VO 

VO 

VC 

VO 

00 

00 

00 

00 

Os 

On 

On 

On 

a 

a 

e*5 

tu 

s 

s 

.§ 

■s 

"o 

"o 

"o 

O 

:r: 

a: 

Ji: 

ON 

1  = 
0\  'rr 

X)  > 

Xj 

< 


c 

0.2 
U  o 
C3  (D 

O  — 


I  I 

C 


!  zn 


4) 

a 

E 

CO 

c/5 

4x 

o 

(D 

Wh 

3 

cd 


(N 

O 

d 

+1 

ON 

CN 

d 


VO 

d 


On 

d 

I 

VO 

(N 


O 

o 


CN 

O 


VO 

o 

d 

+!  . 


\  I  ^ 

CN 

ro 

o  o 


VO 

CN 

o 

d 

+i 

ON 


VO 

VO 

o 


d  o 

£r2 

oq 

do 

+i 

o 

VO 

O 


o 


+1 

CN 

O 

o 


<u 

to 

I 


c 

cO 


c 

E 

’5b 


o 

c: 

•E 

E 

p 


c 

o 

3 

E 

4-1 


Cu 

CO 

U 

E 

p 


3 

OX) 

CO 

E 


c 

0) 

3 

E 

UJ 


eO 

O 

C/5 

CO 

CO 


o 

CO 

o 

■o 

O) 

Q. 

E 

3 

Q 


a. 

< 

c/5 

+ 


3 

'cO 

"cO 

X) 

O 


T3 

> 


I  O 
(/5  t/5 
D  CO 

7=5 

E 

CO 


CO 


CO 


to 

o 

o 

d 

+1  ^ 
VO  to 
m  Q 

o  o 
d  d 

»o 
O  tq 
o  q 
o  o 
+l 
CN 

O 


CN 


4> 

4) 

(D 

X 

o  ^ 

c 

0  C£ 

:z; 

C 

cO 

(D 

C/5 

c 

3 

CO 

X 

eO 

4) 

C/5 

g- 

E 

D 

►7  CO 

--T  ^ 

3  £ 

cO 

C:S 

TD 

Im 

0 

-o 

3 

i2 

cO 

o 

o 

X 

c« 

'C 

U 

X 

C/D 

1 

a: 

E 

CO  t: 

*3  0 

Sz 

CO  t: 
*3  ® 

c/s 

? 

*3 

0 

E 

J 

W) 

*>H 

4) 

E  *0 

CO 

0 

c 

X 

t: 

o 

CA 

c/T 

X 

> 

z. 

3 

3 

c/5 

E 

E 

3 

cO 

PQ 

2 

3 

O 

£ 

CO 

Qu 


Oh  Cl 

E 


C 

(D 

X 

u 

c/a 

p2  o 
a- CO 

E  ^ 

CO  c/a 

C«  Q 

•c 

'4-I 

C/5 

0) 

H 


VO 

O 

o 

? 

CN 

VO 

O 

d 


0 

0 

0 

CN 

00 

00 

ON 

On 

CN 

0 

0 

0 

0 

— .H 

—1 

.  00 

00 

00 

00 

00 

1 

00 

r- 

1 

CN 

VO 

'  ON 

ON 

ON 

ON 

ON 

ON 

ON 

m 

CO 


to 

CN 


CO 

CN 

00 

I 

CN 

ON 


0) 

E 


2 

'5 

E 

3 

C/5 


W5 

VO 

O 

o 

O 

o 

d 

'w' 

o 

o 

d 

+1 

VC 

o 

d 


a\ 

o\ 


in 

in 


2 

es 

iJt 

c 

08 

CA 

ii 

>> 

4) 

(#3 

*s 

4^ 

x" 

O 


i2 

c 

s 

sg 

4> 

C/5 

*U 

E 

3 


Numbers  in  parenthesis  denote  the  range  of  values  reported. 
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Radionuclide  Contaminants  in  Central  Arctic  Sea  Ice 


Debra  A.  Meese  and  Walter  B.  Tucker  III 
U.S.  Army  Cold  Regions  Research  and  Engineering  Laboratory 

Hanover,  NH 

The  U.S./Canada  Arctic  Ocean  Section  (AOS-94),  which  took  place 
during  July  and  August,  1994,  was  the  first  science  expedition  in  100  years  to 
make  a  complete  transect  of  the  Arctic  Ocean.  Investigations  of  oceanography, 
biology,  geology  and  geophysics,  ice  physics,  meteorology  and  contaminants 
were  conducted  by  researchers  on  board  the  USCGC  Polar  Sea  and  the  CCGS 
Louis  S.  St.  Laurent.  The  overall  goals  of  the  expedition  were  to  establish  a 
baseline  data  set  of  ocean  structure  and  ice  properties  relating  to  global  climate 
change  research  in  a  very  sparsely  sampled  part  of  the  Arctic,  to  investigate 
paleoclimate  and  plate  tectonics  in  the  central  Arctic,  and  to  determine  the 
pervasiveness  of  radionuclide  and  persistent  organic  contaminants. 

Our  objective  was  to  determine  whether  radionuclide  contaminants  are 
currently  present  in  central  Arctic  sea  ice.  Should  contaminants  be  observed,  it 
would  indicate  that  ice  is  an  effective  mechanism  for  transporting  and  releasing 
radionuclides  far  from  source  areas.  A  total  of  5  ice  samples  and  12  sediment 
samples  were  collected,  essentially  from  the  ice  edge  in  the  Chukchi  Sea  to  the 
North  Pole,  as  shown  in  Figure  1.  Radionuclide  analysis  of  ice  requires  a 
minimum  of  200  I  of  meltwater,  which  translates  to  about  25  m  of  ice  core.  The 
cores  for  a  single  site  filled  4  large  sample  containers  (plastic  garbage  pails) 
which  were  melted  for  further  processing.  Subsequent  processing  required 
treating  the  meltwater  with  various  tracers  and  co-precipitating  agents,  and 
decanting  the  precipitates.  The  entire  procedure  required  5  -  7  days.  Sediment 
samples  were  simply  scraped  from  the  ice  surface,  collected  from  cryoconite 
holes,  and  from  the  bottom  of  melt  ponds. 

Somewhat  surprising  was  the  fact  that  sediment  was  observed  on  the  ice 
from  the  ice  edge  to  the  Pole.  In  some  areas,  sediment  covered  an  estimated 
10%  of  the  ice  surface.  The  amounts  of  sediment  observed  varied,  yet  on  this 
part  of  the  transect,  there  was  nearly  always  some  in  evidence.  The 
photographs  in  Figure  2  show  examples  of  the  sediment  observed  on  the 
expedition.  In  stark  contrast,  virtually  no  sediment  was  observed  from  the  Pole 
through  the  entire  eastern  Arctic.  This  was  unusual  because  it  has  frequently 


been  observed  in  the  eastern  Arctic  on  previous  expeditions.  It  indicates  that 
sediment  incorporation  is  episodic,  likely  occurring  on  the  shallow  circum-Arctic 
shelves  during  storms  in  the  freezing  season.  The  mechanisms  of  sediment 
incorporation  in  ice  are  not  well  understood.  It  is  thought  that  the  sediment 
suspended  in  the  water  column  is  entrained  into  frazil  by  various  mechanisms. 
The  frazil  later  consolidates  into  a  solid  ice  cover  with  its  sediment  load.  Anchor 
ice  formation  on  the  shallow  seafloor,  may  also  incorporate  sediment  into  the 
consolidating  ice  cover  when  it  rises  from  the  bottom.  It  is  clear  from  our 
observations  that  incorporated  sediment  is  capable  of  being  moved  large 
distances  from  the  areas  of  initial  entrainment  by  the  ice.  This  implies  that 
contaminants  may  be  transported  and  released  throughout  the  Arctic  if  they  are 
preferentially  sorbed  by  the  sediment,  or  incorporated  directly  into  the  ice. 

The  collected  samples  are  now  being  analyzed  for  radionuclides  by  Dr. 
Lee  Cooper  and  Dr.  Jacqueline  Grebmeier  at  Oak  Ridge  National  Laboratory. 
We  are  collaborating  with  Roger  Colony,  University  of  Washington  and  Dr. 
Stephanie  Pfirman,  Barnard  College  to  determine  probable  trajectories  of  the 
ice  and  ice-borne  sediment  samples.  The  sedimentological  analysis  is  being 
conducted  out  by  Dr.  Erk  Reimnitz,  an  expert  in  ice-borne  sediment  from  the 
USGS,  who  was  also  a  member  of  the  AOS-94  expedition.  In  conjunction  with 
the  ice  physical  properties  measurements,  we  should  be  able  to  delineate  likely 
source  regions. 


Figure  1 .  Locations  of  the  ice  and  sediment  samples  collected  for 

radionuclide  testing  on  the  1994  U.S./Canada  Arctic  Ocean 
Section. 


Surface  sediment  on  sea  i 


observed  during  AOS-94. 
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Patterns  of  Monthly  Wind  and  Ice  Motion  in  the  Arctic  Basin 

Roger  Colony,  Polar  Science  Center,  University  of  Washington 

Monthly  fields  of  surface  atmospheric  pressure,  geostrophic  wind,  and  ice  motion 
have  been  analyzed  for  the  years  1979  through  1993.  We  apply  a  principle 
component  analysis  (PCA)(also  known  as  empirical  orthogonal  function  analysis)  to 
the  aggregate  of  the  180  monthly  fields  of  ice  motion.  This  type  of  analysis  extracts 
the  spatial  pattern  most  resembling  the  aggregate  of  monthly  fields.  By  continuing 
the  analysis,  a  sequence  of  patterns,  Pj(x),  may  be  determined  such  that  each 
monthly  field  can  be  constructed  as 


i 

Our  analysis  shows  that  only  three  spatial  patterns  are  required  to  capture  96%  of  the 
variance  of  monthly  fields  of  ice  motion  in  the  central  Arctic  Basin.  Applying  the 
same  analysis  to  the  monthly  fields  of  geostrophic  wind,  we  find  that  another  three 
patterns  capture  97  %  of  the  variance  of  the  monthly  fields.  Furthermore,  the  patterns 
of  ice  motion  and  wind  bear  a  remarkable  resemblance  and  are  temporally  correlated. 

Borrowing  from  recent  research  in  the  atmospheric  sciences,  we  apply  the  single 
valued  decomposition  (SVD)  method  to  the  pairs  of  monthly  wind-  and  ice-motion 
patterns.  The  SVD  analysis  identifies  pairs  of  patterns  which  explain  the  covariance 
between  spatial  patterns  of  wind  and  spatial  patterns  of  ice  motion.  This  analysis 
verifies  that  monthly  wind  anomaly  patterns  are  primarily  responsible  for  the 
monthly  anomalies  in  large-scale  ice  motion.  Initially  used  in  forecasting  and 
predictability  studies,  the  PCA  and  SVD  analyses  are  powerful  tools  in  identifying 
the  large-scale  spatial  patterns  of  several  variables  and  in  synthesizing  particular 
realizations. 

The  findings  are  relevant  to  the  design  of  a  buoy  array  to  monitor  monthly  wind  and 
ice  motion.  We  conclude  that  relatively  few  buoys  are  required  to  characterize  the 
monthly  patterns  in  the  central  Arctic  Basin. 


Figure  !•  Illustration  of  buoy  trajectories,  1979^1993  (the  analysis  will  include  manned  station 
data  for  this  period)  (data  set  is  collection  of  monthly  wind  and  ice  motion  fields) 


Optimal  Monthly  Ice  Velocity  •  May  1980 
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Optimal  Monthly  Ice  Velocity  -  May  1990 
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Statistics  of  Monthly  Wind  and  Ice  Motion 

denote  u{x,  t)  as  monthly  ice  motion  for  each  month. 
t=l,2,3,4,...,180  (Jan.  79  -  Dec.  93) 


•  Annual  Mean 


180 

ll  (x)  —  U  f) 
r=  1 


•  Monthly  Mean 

15 

n(x,x)  =  X  “  3') 

j.  =  i 


Variance  (several) 


var{uix)}^  =  ^  [(u(x,t)  -uix))] 


/=  1 


var{u(x,T)}^=  Y,  y)  -  ^  ^ 

y=  1 

180  ^ 

var{u(x)  } ^  =  [u(x,t)  -  u  (x,  T)  ] 

r  =  1 


var{u  (x)}  =  [uix^z) -uix)] 


i=l 


Figure  4.  Equations  for  means  and  variances,  annual  and  seasonal  cycle. 


Annual  Mean  Pressure  1979-1993 


Annual  Mean  Geostrophic  Wind  1979-1993 


Figure  5b 


Optimally  Interpolated  Current 
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Mean  Sea  Ice  1979-1993 
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Figure  5d 


Pressure  Total  Variance  about  Annual  Mean  in  mb''2 


Figure  6.  Variance  of  monthly  pressure  (departure  from  annual  mean). 


Pressure  Variance  of  Climatological  Annual  Cycle  in  mb''2 


Figure  ?•  Variance  of  mean  seasonal  cycle  of  pressure  (from  annual  mean). 


Geostrophic  Wind  Total  Variance  about  Annual  Mean  in  m''2/s''2 


Figure  8.  Variance  of  monthly  geo.  wind  (departures  from  annual  mean). 


Geostrophic  Wind  Variance  of  Climatological  Annual  Cycle  in  m'^2/s''2 


Figure  9.  Variance  of  the  mean  seasonal  cycle  of  geostrophic  wind. 


Ice  Velocity  Total  Variance  about  Annual  Mean  in  cm''2/s''2 


Figure  10.  Variance  of  the  monthly  ice  motion  (from  annual  mean). 


Ice  Velocity  Variance  of  Climatological  Annual  Cycle  in  cm'^2/s''2 


Figure  11.  Variance  of  the  mean  seasonal  cycle  of  ice  motion. 
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Figure  1.  Represent  U  as  sums  of  products  (separation  of  time  and  space). 


Familiar  Patterns 


Vibrating  String 
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dx 


X  (JC,  t) 
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k=  1 


a^sin  {kx)  exp  {-ikat) 


T  (jc,  0  =  X  {-k^oP'i) 


Figure  2.  Wave  equation,  diffusion  equation. 


Principle  Component  Analysis 
•  remove  mean  annual  field  from  each  monthly  field 
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Figure  3.  Principle  component  analysis. 
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Annual  Ice  Motion 
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Figure  5.  Annual  ice  motion  (show  study  region,  vector  of  600  elements). 


Annual  Geostrophic  Wind 


Figure  1.  Begin  with  the  wind,  first  subtract  the  mean  annual  wind  vector. 


Wind  Vfilocitv  EOF.  Mode  1 


Figure  2a,  Patterns  of  wind,  mode 


Wind  Velocity  EOF.  Mode  2 


Figure  2b.  Patterns  of  wind,  mode  2. 


Wind  Velocity  EOF,  Mode  3 


Figure  2c-  Patterns  of  wind,  mode 


EOF  Geostrophic  Wind,  Departures  from  Annual 


Mode 

Figure  3.  Scalar  multiplier  ALPHA,  modes  1,  2,  and  3. 


EOF,  Multiplier  of  Wind,  mode  = 
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Figure  5a.  Patterns  of  ice  motion. 
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Figure  5b,  Patterns  of  ice  motion,  mode 


Ice  Velocity  EOF.  Mode  2 


Figure  5c.  Patterns  of  ice  motion,  mode 


Ice  Velocity  EOF.  Mode  3 
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Figure  5d.  Patterns  of  ice  motion,  mode 


EOF  Ice  Velocity,  Departures  from  Annual 


Figure  6.  Scalar  multiplier  ALPHA,  modes  1,  2,  and  3. 
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Figure  7.  Time  series  of  pattern  multiplier,  modes  1,  2,  and 


CONCLUSION 


1992  EOF,  Multiplier  of  Wind,  mode  =  1 ,  Departures  From  Monthly  Means 
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Figure  9a,  Comparison  of  patte 
(these  multipliers  are  for  depa 


1992  EOF,  Multiplier  of  Velocity,  mode  =  1,  Departures  From  Monthly  Means 
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Figure  9b.  Comparison  of  pattern  multipliers,  ice  motion 
(these  multipliers  are  for  departures  from  mean  seasonal  cycle) . 


COVARIANCE  OF  PAIRS  OF  FIELDS 


Correlation  Coefficients 


Covariance  of  x  and  y  (scalars) 


cov{x,y}  =E{{x-x)  iyi-y) 
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Covariance  of  x  and  y  (vectors) 


C  =  £  {  (X  -  X)  ( J  -  J')  }  =  matrix  with  elements 
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Frobenius  norm 


m  n 
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Figure  1.  Define  covariance  for  a  pair  of  vectors. 


Wind  Velocity  SVD,  Mode  1 


Figure  2a.  Pattern  pairs  for  wind  and  ice  motion,  mode  1. 


Wind  Velocity  SVD.  Mode  2 


Figure  2b.  Pattern  pairs  for  wind  and  ice  motion,  mode  2. 


Ice  Velocitv  SVD,  Mode  3 


Wind  Velocitv  SVD.  Mode  3 
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Pattern  pairs  for  wind  and  ice  motion,  mode  3. 


Wind-Ice  Multiplier.  Mode  1 _ 


Figure  3a.  Interpretation  as  U  =  AG  for  Mode  1. 


Turning  Angle.  Mode  2 


Figure  3b.  Interpretation  as  U  =  AG  for  Mode  2. 


Figure  3c.  Interpretation  as  U  =  AG  for  Mode  3. 


SVD  Velocity  -  Geostrophic  Wind,  Departures  from  Annual 
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Figure  h.  Scalar  multiplier  ALPHA,  modes  1,  2,  and 


SVD,  Multiplier  U-G,  mode  = 
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Figure  5.  Time  series  of  pattern  multiplier,  modes  1,  2,  and 


1992  SVD,  Multiplier  U-G,  mode  =  1,  Departures  From  Monthly  Means 
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CONCLUSIONS 


Things  to  Do 

•  EOF  and  SVD  analysis  of  monthly  ice  motion  strain 

•  EOF  and  SVD  analysis  of  microwave  data 

•  Use  monthly  pressure  fields  1953-1978  to  synthesize  ice  motion  fields 
(validate  with  obs  from  NP,  DARMS,  ARLIS) 

•  Use  coastal  wind-pressure  data  1920’s  -  1990  to  synthesize  ice  motion  fields 

•  Investigate  systematic  differences  between  wind  driven  patterns  and  obs.  patterns 

•  Investigate  coastal  regions  -  Fram  Strait 

•  Spectral  analysis  of  A(t),  B(t) 

(predictability) 


Russian  Historical  Ice  Charts  and  Dataset 

1967-1990 


Ignatius  Rigor 
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Russian  Historical  Ice  Chart  Dataset,  1967  - 1990 

Analysis  of  ice  conditions  for  the  arctic  regions  have  been  prepared  to  support  high  lati¬ 
tude  sealing  and  whaling  activitied  since  the  17th  century.  These  early  charts  docu¬ 
mented  the  position  of  the  ice  edge,  conditions  within  the  ice  pack,  and  access  to  arctic 
ports.  Ice  charts  for  specific  regions  of  the  Arctic  Basin  began  to  be  routinely  prepared  in 
the  1 940’s  in  support  of  commercial  and  military  requirements.  Presently,  operational  ice 
charts  for  the  arctic  are  prepared  at  the  Arctic  and  Antarctic  Research  Institute  (AARI), 
the  US  National  Ice  Center  (NIC),  and  the  Norweigian  Meteorological  Institute. 

Modern  analyses  of  ice  conditions  are  based  on  a  combination  of  observations,  ice  fea¬ 
ture  persistence,  mathematical  models,  and  climatology.  The  most  important  ingredients 
are  the  observational  data-  ship  and  shore  station  reports,  aerial  reconnaissance,  and 
satellite  imagery.  Current  satellite  imagery  makes  global  ice  analysis  possible.  Much  of 
the  earlier  (1966-1972)  satellite  data  were  low  resolution  visual  and  infrared  imagery  lim¬ 
ited  by  the  polar  darkness  and  cloud  cover.  In  1972  the  Advanced  Very  High  Resolution 
Radiometer  (AVHRR)  provided  one  km  resolution  in  the  visual;  this  sensor  is  still  used  in 
the  preparation  of  ice  charts.  In  1973  the  US  launched  the  Electronically  Scanning  Micro- 
wave  Radiometer  (ESMR)  which  provides  for  all  weather  remote  sensing  capability.  The 
Russians  have  placed  great  emphasis  on  active,  all  weather,  microwave  systems,  such 
as  airborne  and  satellite  based  Sideways  Looking  Airborne  Radar  (SLAR). 

The  Arctic  and  Antarctic  Research  Institute  began  preparing  routine  ice  charts  in  1947, 
to  support  shipping  along  the  Northern  Sea  Route.  Since  1964  AARI  has  prepared  10 
day  ice  charts  throughout  the  Arctic  Basin,  emphasizing  the  Kara,  Laptev,  and  East  Sibe¬ 
rian  Seas. 

A  digital  form  of  these  ice  charts  is  available  as  a  sea-ice  dataset  in  Sea-ice  Data  in  Dig¬ 
ital  Form  (SIGRID)  format.  The  Commission  for  Marine  Meteorology  of  the  World  Meteo¬ 
rological  Organization  (WMO)  recognized  the  need  for  a  uniform  data  base  to  facilitate 
the  archival  and  distribution  of  sea-ice  data.  In  1989,  the  WMO  adopted  the  SIGRID  for¬ 
mat,  see  WMO  report  No  558  (19ZZ).  MANUAL  ON  MARINE  METEOROLOGICAL 
SERVICES,—  The  SIGRID  format  was  adopted,  in  part,  because  it  was  able  to  focus  on 
conventional  sea-ice  variables  while  accommodating  flexibility.  SIGRID  data  are 
encoded  on  a  geographical  grid.  For  the  above  marginal  seas,  the  standard  resolution  is 
1/4  degree  of  latitude  and  1  degree  of  latitude  north  of  75N  latitude  and  1/4  degree  of  lat¬ 
itude  and  1/2  degree  of  longitude  south  of  75N.  In  the  Laptev  Sea,  for  instance,  this 
amounts  to  about  1200  grid  points. 

For  each  grid  point  the  sea-ice  variables  are  divided  into  ten  categories.  Each  data  cate¬ 
gory  may  contain  specific  and  detailed  information,  or  it  may  be  flagged  as  containing  no 
reported  data.  Most  of  the  information  in  the  AARI  SIGRID  charts  is  found  in  the  concen¬ 
tration/stage/form  category.  This  information  includes  total  ice  concentration  and  partial 
concentration,  stage,  and  form  of  the  three  thickest  ice  types.  The  encoded  form  informa¬ 
tion  includes  floe  size  estimates,  bands  and  patches  of  ice,  and  fast  ice. 


Concentration  /  Stage  /  Form  category. 

This  information  includes  total  ice  concentration  and  partial  concentration,  stage, 
and  form  of  the  three  thickest  ice  types. 
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Transpolar  Drift  -  The  garbage  scow  from  Siberia? 

Erk  Reimnitz,  U.S.  Geological  Survey 

A  video  by  the  dirty-ice  group  of  the  USGS  demonstrates  drift-ice  related  processes 
on  the  shallow  shelf  and  in  the  coastal  zone.  That  work  has  shown  that  the  classic  fast-ice 
entrainment  mechanisms:  aeolian  deposition,  river  flooding,  slumping  from  coastal 
cliffs,  bottom  adfreezing,  and  ice  plowing,  are  insignificant  compared  to  the  mechanism 
of  suspension  freezing  in  supercooled  water.  That  mechanism  is  most  effective  on 
shallow,  circum- Arctic  shelves  landward  of  the  30-m  isobath  (Figure  1).  Comparing 
sediments  in  sea  ice  of  the  Beaufort  Gyre  (Figure  1)  to  those  of  the  bordering  <30-m  deep 
shelf  surface  to  assess  sorting  by  entrainment,  we  find  them  much  finer  (10%  vs.  48% 
sand),  having  a  much  higher  organic  carbon  content,  and  lower  carbonate  content  than 
those  of  the  source  surface.  Clay  mineral  composition  and  lithology  of  sand  grains  do  not 
show  differences  that  suggest  suspension  freezing  is  selective.  No  data  are  available 
about  radio-nuclide  entrainment,  but  we  may  assume  an  affinity  for  fine  particulate  and 
organic  matter. 

In  the  Beaufort  Gyre  off  North  America,  sediment  loads  are  highly  variable,  ranging 
upward  to  over  1,000  tonnes/km^,  depending  on  whether  the  necessary  strong,  freezing 
winds  occurred  during  the  generally  short  open-water  period  in  late  fall.  The  loads  can  be 
much  higher  than  the  sediment  supply  from  local  rivers  feeding  the  area,  and  therefore 
lead  to  shelf  erosion.  To  emphasize  the  importance  of  sediment  transport  by  sea  ice  for 
the  Arctic  Ocean  sediment  budget  we  compare  the  above  sediment  load  in  sea  ice  to  the 
low  annual  denudation  rate  of  about  10  tonnes/km^  on  continents  surrounding  the  basin. 
In  years  between  such  sediment  pulses  spaced  5  to  10  yrs  apart,  the  ice  of  the  Beaufort 
Gyre  has  been  observed  and  measured  to  be  very  clean. 

In  the  Laptev  Sea  of  Siberia,  which  is  dominated  by  a  major,  perennial  polynya 
(Figure  2),  conditions  for  suspension  freezing  to  occur  apparently  are  met  yearly, 
resulting  in  a  rather  continuous  sediment  flux  to  Fram  Strait  and  the  North  Atlantic.  Very 
preliminary  calculations  of  sediment  flux  in  sea  ice  actually  exported  from  the  Laptev  Sea 
are  in  the  same  order  of  magnitude  as  sediment  supply  by  the  h&na.  River  feeding  it. 

Very  little  is  known  about  pollutants  in  Arctic  pack  ice.  Many  mud  samples 
collected  from  sea  ice  by  the  USGS  in  different  expeditions  to  the  Beaufort  Gyre  and  the 
Siberian  Branch  of  the  Transpolar  Drift  over  the  last  5  years  are  undergoing 
sedimentological,  mineralogical,  and  paleontological  investigations.  These  samples  are 
also  being  compared  to  surficial  deposits  in  the  deep  basins.  Portions  of  26  of  these 
valuable  samples  (Figure  3)  were  originally  offered  to  the  ONR  program  for  radio  nuclide 
analyses,  but  are  presently  being  analyzed  through  funding  from  other  sources. 


Figure  1.-  Generalized  Arctic  Ocean  ice  circulation  pattern,  showing  the 
sediment  sources  on  shelf  areas  <30  m  deep. 


March  1992 


Figure  2.  -  The  fast  ice  and  perennial  flaw  lead  with  ice  thinner  than  15  cm  as 
mapped  for  the  period  March  20  to  24,  1992  by  the  Arctic  and  Antarctic 
Research  Institute  in  St,  Petersburg. 
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Figure  3.  -  Locations  of  26  mud  samples  collected  from  pack 
ice  by  the  U.S.G.S.  from  different  ships  during  a  4-year  time 
span,  and  presently  being  analyzed  for  radio  nuclide  contamination. 
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Data  on  the  Ob  River  System, 
Yenisei  River  System  and  Kara  Sea  Region: 

What  we  have,  don’t  have...  and  need. 

Peter  Becker 

CCPO,  Old  Dominion  University,  Norfolk,  Virginia 

Abstract 

At  the  outset  we  were  concerned  that  a  dearth  of  actual  field  data 
on  the  Ob  and  Yenisei  Rivers  and  Kara  Sea  would  limit  our  ability  to 
establish  baseline  values  for  the  various  parameters  necessary  to  help 
model  the  transport  of  radionuclides  by  the  large  Arctic  rivers  to  the 
Kara  Sea  and  Arctic  Ocean.  Instead  we  find  that  there  are  both  many 
different  data  sets  and  a  rich  literature  to  assist  us  and  that  much  is 
in  english  or  english  translation. 

These  data  sets  include,  but  are  not  limited  to: 

•  Drift  Bouy 

•  Hydrographic 

•  Hydrological 

•  Meteorological 

•  Radionuclide  Surveys 

•  Tidal  Constants 

The  Drift  Buoy  data  encompass  both  US  and  Russian  sets  (eg.  Fig. 
la,b  with  two  of  the  Russian  sets  indicated  by  Fs  in  the  Kara  Sea.) 

The  Hydrographic  data  sets  include  US  icebreaker  cruises  in  1965 
and  1967.  Other  data  sets  include: 

•  Ostliind  et  al.  (Tritium) 

•  Anderson  et  al.  (salinity/temperature/etc. 

ODEN  data) 

•  Schlosser  et  al.  (^®0/Salinity;  ^i7:He 
ARKTICA-IV  Data) 

•  Gloersen  et  al.  (Seaice  passive  microwave 
1978-1987) 

•  Wallace  et  al.  (CFC’s) 

The  Hydrological  and  Meteorological  and  Radionuclide  data  sets 
include  Russian  data  through  1992  at  stations  in  Figure  2, 

Tidal  constants  are  available  for  all  the  stations  in  Figure  3  as  well 
as  200  other  stations  about  the  Arctic. 

A  Bibliography  of  over  600  citations  has  been  compiled  and  placed 
along  with  all  open  data  sets  in  a  anonymous  FTP  archive  at: 

ftp.ccpo.odu.edu  cd/pub/becker/IAEA 


The  Bibliography  contains  all  of  the  relevant  translated  Russian 
Literature  citations  from: 

•  Oceanology 

•  Problems  of  the  North 

•  Soviet  Hydrology 

•  Soviet  Geography 

•  Soviet  Geology  and  Geophysics 

•  Water  Resources 

as  well  as: 

•  Annual  Review  of  Energy 

•  Health  Physics 

•  J.G.R. 

•  J.P.O. 

•  Deep  Sea  Research 

•  Others... 

The  literature  citations  cover  a  range  of  topics  as  seen  in  Figure  4. 

Intensive  Russian  studies  of  the  Kara  Sea  resulted  from  the  eco¬ 
nomic  necessity  of  predicting  the  ice  conditions  on  the  Northern  Sea 
Route  (Figure  5)  and  resulted  in  numerous  expeditions  starting  before 
1928  and  covering  the  entire  Arctic  Ocean.  This  included  ice  stations 
(Figure  6)  and  numerical  circulation  models  (Figure  7).  and  coastal 
observations  (Figure  8). 

The  literature  also  includes  citations  on  the  effect  of  anthropogenic 
influences  on  the  river  discharge  of  the  Ob  and  Yenisei  (Figure  9) 

We  have  recently  succeded  in  obtaining  higher  resolution  data  sets 
for  the  Ob  and  Yenisei  Rivers  including  the  following: 


•  Hydrographic  Data 

•  Hydrologic  Data 

•  Kara  Sea 

•  Ob  River 

•  Radioisotopes 

•  Suspended  Sediments 

•  Yenisei  River 


These  data  sets  will  assist  in  improving  our  understanding  of  the 
Ob  and  Yenisei  river  transports  and  their  interaction  with  the  Kara 
Sea  and  Arctic  Ocean. 


In  conclusion: 


•  Data  on  hand  is  extensive  at 

mean  monthly  time  and  basin  size  scales. 


•  Data  on  hand  is  now  adaqiiate  to  resolve 
short  term  flow  (<  30  days)  and  storms. 


•  Data  on  hand  is  inadaquate  to  resolve 

features  of  estuary  stratification. 

More  data  is  needed  hut  it  is  available. 

•  Radionuclide  data  available  do  not  yet  fully  resolve 

nuclide  transfer  from  river  to  estuary  to 
ocean. 


•  Impact  of  reservoirs  is  significant  and 
must  be  built  into  models. 
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Figure  1.  Map  of  the  Ob  and  Yenisei  Rivers  and  their  watersheds.  Hie  hatched  lines  denote 
the  watershed  boundaries.  Potential  sources  of  radionuclide  contaminants  are  marked  with 
triangles.  Many  of  the  meteorological,  hydrological,  and  hydrographic  stations  that  are  being  used 
in  the  riverine  modeling  are  shown. 


Figure  3 


Figure  1.  Map  of  the  Ob  and  Yenisei  Rivers  and  their  watersheds.  The  hatched  lines  denote 
the  watershed  boundaries.  Potential  sources  of  radionuclide  contaminants  are  marked  with 
triangles.  Many  of  the  meteorological,  hydrological,  and  hydrographic  stations  that  are  being  used 
in  the  riverine  modeling  are  shown. 


Literature  Sets: 


•  Estuaries 

•  Freshwater  Ice 

•  Geology 

•  Ground  Water 

•  Interannual  Variation 

•  Kara  Sea 

•  Lena  River 

•  Miscellaneous 

•  Ob  River 

•  Permafrost 

•  Radioisotopes 

•  River  Models 

•  Sea  Ice 

•  Suspended  Sediments 

•  Yenisei  River 


Map  3.1  The  Northern  searoute 
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Table  1.  Suspended  matter  flux  in  the  mouths  of  rivers 
regulated  by  reservoir  systems  [2] 


System 

Numberof 

reservoirs 

Average  annual  suspended  matter 
discharge  across  the  mouth, 
in  million  tons 

prior 

to  regulation 

after 

regulation 

Volga  j 

7 

25.9 

7-8 

Kama 

2 

10.0 

1.3-  1.9 

Dnieper 

6 

2.0 

0.4 

Irtysh 

2 

11.2 

2.5  -  3 

Yenisey 

2 

12.4 

1.8 -2.2 

Angara 

3 

2.4 

0.5 

(after  A/alyan,  1988)  . 
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CIRCULATION  OF  THE  ARCTIC  SEAS  BASED  ON 
RUSSIAN  SOURCES  OF  INFORMATION 

Short  review  of  sources  of  information  and  methods  of  simulation  of  the  Arctic  Seas  water 
circulation  is  presented.  Summary  currents  in  Arctic  Seas  axe  composed  of  permanent, 
\^^nd-induced  and  tidal  currents.  In  the  ICara,  East-Siberian  and  Chukchi  Seas  permanent 
currents  are  formed  as  a  result  of  fresh  water  inflow  from  rivers  and  water  exchange  with 
adjacent  areas.  These  currents  in  the  Arctic  Seas  have  a  relatively  stable  pattern  The  rate 
of  the  permanent  currents  is  0. 1-0.3  knots  in  most  of  the  Arctic  Seas  area.  In  some  regions, 
mainly  in  straits  and  adjoining  mouths  of  large  rivers,  the  rate  of  the  permanent  currents 
increase  to  0.4  -  1.2  knots. 

Schemes  of  wind-driven  currents  in  the  surface  layer  of  Arctic  Seas  for  steady  (i.e.  lasting 
not  less  than  12  hours)  northeast,  southeast,  southwest  and  northwest  winds  with  a  speed 
10  m/s  were  published  in  1992. 

In  the  Laptev,  East-Siberian  and  Chukchi  seas,  southeast  wind  drives  a  current  flowing 
mainly  northwest  at  a  magmtude  of  0.4-0. 8  knots.  In  the  Kara  Sea  similar  wind  drives 
a  current  flowing  mainly  northeast  at  the  same  magmtude.  In  Prohv  Dmitriya  Lapteva, 
Prohv  Sannikova  and  Proliv  Longa  currents  flow  west  and  in  Proliv  Kaxskiye  Vorota  and 
in  Bering  Strait  they  are  flow  northeast.  In  the  Bering  Strait  current  velocity  ranges  from 

1  to  1.2  knots.  In  Proliv  Vil’kitskogo  southeast  wind  drives  currents  flowing  east  and  and 
northeast  at  1.3  knots. 

Southwest  wind  drives  currents  flowing  predominantly  to  the  east  and  northeast.  Cur¬ 
rent  speed  in  the  open  sea  ranges  from  0.3  to  0.8  knots,  in  the  Bering  Strait  is  1.7  knots, 
in  Proliv  Sannikova  is  2.1  knots,  and  in  Prohv  Yugorskiy  Shax  is  2  knots.  In  Prohv 
Vil’kitskogo  this  wind  drives  a  current  east  at  1.3  knots. 

In  the  Laptev,  East-Siberian  and  Chukchi  seas,  northwest  wind  drives  currents  flowing 
southeast  at  0.4  to  0.6  knots.  In  the  Kara  Sea,  similar  wind  drives  currents  flowing  south 
and  southwest  at  0.2  to  0.5  knots.  Current  is  south-going  in  Bering  Strait  reaching  a  speed 
of  1.1  knots. 

Maximum  speed  of  the  wind  currents  is  generally  1.5  -  2  knots,  in  some  places  reaches 
2.5  to  3  knots,  and  may  be  as  much  as  5  to  6  knots  in  Prohv  Matochkin  Shax,  Prohv 
Yugorskiy  Shax  and  Prohv  Matisena.  Cmrents  of  such  speed  axe  usually  observed  dmdng 
winds  of  20  m/s  and  greater. 

Tidal  currents  in  the  Arctic  Seas  axe  dominantly  semidiurnal.  In  the  narrow  straits  and 
bays  they  have  a  reverse  character,  i.e.,  they  change  direction  every  6  hours  on  180".  In  the 
open  regions  of  the  Kara  and  western  Laptev  Seas,  the  spring  tide  velocities  range  from  0.4 
to  0.8  knots.  In  the  eastern  Laptev  and  in  the  East-Siberian  Seas,  tide  currents  are  about 
0.2-0.4  knots  and  in  the  Chukchi  Sea  velocities  range  from  0.4  to  0.6  knots.  In  several 
regions  of  the  Arctic  Seas  tidal  currents  reach  extremum  and  have  velocities  of  over  1.5  - 

2  knots.  For  example,  measurements  show  velocities  of  1.5-1. 8  knots  in  the  Baidaratskaia 
Bay,  to  the  north  of  Behi  Island,  in  the  Obskaia  Cuba,  and  in  the  region  of  Khatanga  Bay. 
In  the  straits  of  Karskie  Vorota,  Yugorskii  Shax,  Malygina  Strait  and  Nordvik  Bay,  they 
reach  2-2.8  knots. 


Project: 


fflGH  RESOLUTION  MODELING  AND  REMOTE  SENSING 

OF  THE  EAST-SIBERIAN  SEA  AND  ALASKAN  SHELF 

P.L  -  M.  Johnson,  Institute  of  Marine  Science  UAF 

CO-P.I.  -  A.  Proshutinsky,  Institute  of  Marine  Science  UAF 

(modeling) 

CO-P.I.  -  J.  Groves,  Geophysical  Institute  UAF 

(remote  sensing) 

I.  CIRCULATION  OF  THE  ARCTIC  SEAS  BASED  ON 
RUSSIAN  SOURCES  OF  INFORMATION 

1.  Sources  of  information 

2.  Observed  circulation 

3.  Numerical  models 

4.  Simulated  circulation 

5.  Two-dimensional  transport  of  pollutants 
in  the  East-Siberian  and  Chukchi  Seas 

6.  Plans 


1.  Objective: 

The  objective  is  to  understand  the  coastal  circulation  from  the  Kolyma  River 
on  the  eastern  Siberian  shelf  to  Pt.  Barrow  Alaska.  The  domain  includes  the 
eastern  Siberian  shelf,  the  Chukchi  Sea,  Bering  Strait,  Bering  Sea  and  part 

and  a  section  of  Alaska’s  shelf.  This  region  was  determined  to  be  an  important 
part  of  the  Arctic’s  circulation  as  it  relates  to  the  potential  pathway  for 
radionuclides  or  other  pollutants. 

2.  Motivation: 

Motivation  stems  from  the  hypothesis  that  the  eastern  Siberian  shelf  provides 
the  boundary  conditions  upstream  of  the  Chukchi  Sea  and  Alaska’s  coast. 
Speculation  is  that  most  of  the  eastern  Siberian  shelf  volume  transport  is  in 
the  eastward  flowing  Siberian  Coastal  Current,  a  narrow  boundary  current 
that  roughly  parallels  the  coastline  until  reaching  Bering  Strait. 

Yet  the  lateral  extent  of  the  Siberian  Coastal  Current,  its  alongstream 
coherence,  and  its  seasonal  variability  are  all  poorly  known. 

3.  Approach: 

The  dominant  physical  mechanisms  in  this  region  will  be  examined  by  using 
several  2  and  3-d  models  with  high  resolution  to  determine  main  features 
of  the  Siberian  Coastal  Current  (buoyancy  forcing  due  to  Kolyma  River 
discharge,  wind  and  ice  cover  effects). 


Sources  of  information: 


Numerical  models 

(time-dependent) 


a 

U  00 
J-H  ^ 

•1-H  cd 

o  ^ 
•S  ^ 

S  -5 
CZ 

O  o 


<1^ 

c3  >5 

"S 

>  S  o 


> 

o  o 


^  o  w  u  ^ 

^  D  Q  3  & 

ra  ^  -O  P 

^  S  ® 

O  O?  <  £  o 


cd  C 

<u  c3 

00  >  D 

•n  O  u 


> 


^  O  S3 

o  «  «  .,  c 


«iS  o  S  2 

3  5  -P  o  -g 

^  ^  P  C  o 

U  »-!  O  Vh 

<  CU  Oh 


S  C/3 
cd 

CO  0^ 

D  GO 
OJQ 

5:2  O 

^  -zs 

CO  Q 

n  H 


O  o 
CO  -S 


< 

H- 

•s  y 

o  M 

c/a  O 

0^ 

•1^ 


«.  CO 

>  >  s 

O  S 

fsi 

2^‘S  2 

pL,  Q  Oh 


c: 

o 

s 

•fH 

3 

<u 

o 

o 

o 

*zl 

.2 

o 

o 

p 

3 

c 

•s 

p 

o 

Oh 

D 

C4H 

CO 

o 

ill 

'o  “  2? 

bO  §  o 

.§  fa  'S 

S  3  « 

W  o  ^ 

o  >  « 

c2  .2 


(A 

<i)  'S' 

■n  C 
O  "O 
E  s=3 

c  <D 

OS 

CB  U 
0*2 
^  .2 


JS  c  ^ 

0}  U  V3 

S  .2 

Sp  d 

o.  2 

O  I  pC^ 

VO  CO 

c  »  2 

:z:  c 


O 

CO  ^ 
Vh  2 

(D  C 

S 


cd 

<D 

CO  > 


•c;  ^  2  is 

«l  6^ 

Ou 


S  -5  c 
c^42  "2 

ed  P  O 

S  2  O 

^■s  Q 


c: 

.2 

o 

•a 

cd 

.2 

CO 

cd 

ffl 

3 

.S 

a 

2 

cd 

JD 

> 

i 

o 

a 

a 

2 

a 

> 

<D 

CO 

o 

c 

*C 

3 

o 

< 

*T3 

PU 

Oh 

FIGURE  3-7.--  Surface  currents  in  the  Arctic  Seas  (Atlas  Arktiki,  1985 ). 


FIGURE  3-8. —Surface  currents  in  the  Arctic  Seas  (Sovetskaia  Arktika,  1970) 

I  -first  type  of.  water  circulation,  II  -  second  type  of  water  circulation; 
a  -  velocities  are  less  than  0.1  knots  (18.5  cmis); 
b  -  velocities  are  greater  than  0.1  knots. 

Currents:  1  -  Transarctic;  2  -  Nordkapskoe;  3  -  Pribrezhnoe  Murmanskoe;  4  -  Murmanskoe; 
5  -  Belomorskoe;  6  -  Yan-Mayenskoe;  7  -  Kolguevskoe;  8  -  Kaninskoe;  9  -  Litke; 

10  -  Zapadno-Novozeme’skoe;  11  -  Medvezgie;  12  -  Vostochno-Shpitsbergenskoe; 

13  -  Vostochno-NovozemeV  skoe;  14  -  Yamal’  skoe;  15  -  Ob’ -Yeniseyskoe;  16  -  Sv.  Anny; 

17  -  Zapadno-Taimyrskoe;  18  -  Vostochno-Taimyrskoe;  19  -  Lenskoe;  20  -  Novosibirskoe; 

21  -  Sedova;  22  -  Tikhokeanskoe;  23  -  Geral’ dovskoe . 


FIGURE  3-9.— Scheme  of  currents  in  0-10  m  layer. 

1  -  direction  of  a  current;  2  -  direction  of  a  supposed  current;  3  -  velocity  c"'  current,  knots 
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Adopted  from  Krutskikh,  1978  (First  type  of  ice  and  water  dynamics). 


Arctic  Seas  under  influence  of  anticyclonic  atmosphere  circulation  (Figure  3-12b). 
Adopted  from  Krutskikh,  1978  (Second  type  of  ice  and  water  dynamics). 


FIGURE  3-13.  Surface  permanent  currents  of  the  Barents  Sea  .  Adopted  from 
Tantsiura  (1959)  and  Novitsky  (1961). 

Currents: 

1  -  Near  coastal  current  of  the  Frants  Josef  Land ;  2  -  Barents  current  (East  Spitsbergen). 
3  -  South  Nord  Cap  current;  4  -  Bear  current  (Bear-hope  current);  5  -  Persei  current; 

6  -  Central  current;  7  -  Lithe  current;  8  -  South  Spitsbergen  current;  9  -  Nord  Cap  current, 
10  -  North  branch  of  Nord  Cap  current;  11  -  Central  branch  of  Nord  Cap  current; 

12  -  South  branch  of  Nord  Cap  current;  13  -  Murmansk  current;  14  -  Nurmansk  coastal 
current;  15  -  Kaninskoe  current;  16  -  Kolguevo-Pechorskoe  current;  17  -  Novozemelskoe 
current;  18  -Coastal  branch  of  Novozemelskoe  current;  19  -  East  branch  of  Novozemelskoe 
current;  20  -  West  branch  of  Novozemelskoe  current;  21  -  White  Sea  current; 

22  -  Pechora  current. 


geostrophical  balance  method .  Adopted  from  Potanin  et  al  (1989). 

Definition  of  the  currents  (numbers)  is  the  same  as  in  Figure  3-13. 

Numbers  near  arrows  are  calculated  velocities,  cm/s 
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FIGURE  3-15.  Surface  permanent  currents  of  the  Barents  Sea  simulated  using 
Burkov’s  model  (Burkov,  1980). 

Definition  of  the  currents  (numbers)  is  the  same  as  in  Figure  3-13. 

Numbers  near  arrows  are  calculated  velocities,  cm/s 


FIGURE  3-16.-  Surface  currents  of  the  Chukchi  Sea 
Adopted  from  Proshutinsky,  1986.  Numbers  are  velocities 


of  currents,  ends. 


FIGURE  3-17.—  Integral  currents  of  the  Chukchi  Sea. 
Result  of  numerical  simulation.  Adopted  from  Proshutinsky, 
1986. 


FIGURE  3-17a.--Scheme  of  maximum  velocity  of  summary  currents  in  0-10  m  layer. 

Navigation  period. 

(1)  -  maximum,  currents,  knots; 

(2)  -  local  maximum  of  current. 


-  direction  of  a  current;  2  -  direction  of  a  supposed  current;  3  -  velocity  of  a  current,  knots. 


FIGURE  3 -17c.  —Scheme  of  currents  in  0-10  m  layer  driven  by  NE  10  m/s  wind. 

1  -  direction  of  a  current;  2  -  direction  of  a  supposed  current;  3  -  velocity  of  a  current,  knots. 


-  direction  of  a  current:  2  -  direction  of  a  supposed  current;  3  -  velocity  of  a  current,  knots. 


20*  W 


FIGURE  3-17e.  -Scheme  of  currents  in  0-10  m  layer ,  driven  by  NW  10  mis  wind. 

1  -  direction  of  a  current;  2  -  direction  of  a  supposed  current;  3  -  velocity  of  a  current,  knots. 


Result  of  numerical  simulation.  Adopted  from  Kowalik  and  Proshutinsky, 
1994. 


FIGURE  3-19.-Mean  velocity  of  spring  tidal  currents,  knots. 


es  of  ice  particle: 


;  Perio' 


01.01 


Trajectories  of  ice  particles;  real  wind;  ice;  slope  1.1  E— 6 


Trajectories  of  Ice  particles;  real  wind;  ice;  slope  1.1  E— 6 
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Time  Series  Observations  of  Seabed,  Fluid,  and  Flow  Conditions 
in  the  Barents  and  Norwegian  Seas 


A.R.M.  Nowell,  R.W.  Sternberg,  K.  Aagaard 
School  of  Oceanography,  WB- 1 0 
University  of  Washington 
Seattle,  WA  98195 

D.  Cacchione 
U.S.  Geological  Survey 
345  Middlefield  Rd 
Menlo  Park,  CA  94025 

R.  Wheatcroft 

Woods  Hole  Oceanographic  Institution 
Woods  Hole,  MA  02543 

This  presentation  summarizes  long-term  observations  carried  out  from  moored  instruments 
and  benthic  tripods  located  in  the  vicinity  of  Svalbard  and  downslope  of  the  Bear  Island  Trough. 
Measurements  from  moored  instruments  include  (1)  current  speed,  direction,  conductivity, 
temperature,  and  depth  at  two  levels  in  the  water  column  at  locations  in  StonQorden  and  east  of 
Hopen  Island  in  the  Barents  Sea  and  (2)  current  speed,  direction,  conductivity,  temperature,  and 
light  transmission  at  20  m  above  the  seabed  at  two  locations  in  approximately  700  m  depth  west  of 
Bear  Island  in  the  Norwegian  Sea  (Figure  1).  Measurements  from  benthic  tripods  were  made 
within  2.2  m  of  the  seabed  and  include  current  speed,  direction,  pressure,  conductivity, 
temperature,  light  transmission  or  scattering,  and  bottom  photographs  at  stations  located  in 
StorQorden,  east  of  Hopen  Island,  and  Olgastretet.  Additionally,  a  benthic  camera  tripod  was 
located  in  Olgastretet  (Figure  1). 

Generally,  all  sensors  were  sampled  hourly  and  the  length  of  time-series  varied  from 
3-15  months.  All  instruments  were  not  deployed  simultaneously;  the  total  experiment  duration 
was  from  June  1991  through  July  1994.  At  present,  initial  data  reduction  has  been  carried  out  and 
detailed  analysis  and  interpretation  is  underway.  This  study  is  sponsored  by  the  Marine  Geology 
and  Geophysics  Program  of  the  Office  of  Naval  Research. 
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Figure  1.  Station  locations. 


SEDIMENT  TRANSPORT  MEASUREMENTS 
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SEDIMENT  TRANSPORT  MEASUREMENTS 
IN  A  REGION  OF  SEASONAL  ICE  COVER 

Objectives: 

Understand  the  transport  of  sediment  under  the 
seasonal  ice  cover  in  arctic  environments 

•  time  varying  concentrations  of  suspended 
sediment 

•  modifications  in  small-scale  bed 
roughness 

•  erosion  and  dispersal  of  water  and 
sediment 

•  forcing  mechanisms 

•  sediment  transport  modeling 


OBSERVATION  SUMMARY 


BARENTS  SEA 

•  STORFJORDEN 
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UW  tripod 
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Suspended  sediment  0-3  mab 

uses  GEOPROBE  tripod  (Summer  '92) 
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P  2  mab 

Sediment  concentration  2  levels  <1 .2  mab 


CTD  1 83  m 

UV  182  m 

UW  tripod 

UV  1  mab  2  mab 

CTD  1  mab 

Camera  2  mab 
P  2  mab 

Suspended  sediment  0-3  mab 

WHO!  tripod 

Bottom  photographs 

•  OLGASTRETET 

GEOPROBE  tripod 

UV  4  levels  <1.2  mab 
CTD  1  mab 

Camera 
P  2  mab 

Sediment  concentration  2  levels  <1 .2  mab 

Suspended  sediment  tripod 
Particle  size 
U,V 

Particle  settling  velocity 

NORWEGIAN  SEA 

Mooring  1,  2,  3 

U,  V,  temperature,  conductivity 
Suspended  sediment  concentration  20  mab 


•  HOPEN  ISLAND 

PMEL  (APL)  mooring 

CTD  156  m 
UV  155  m 


Storfjorden:  UW  Aanderaa 


Figure  1.  Example- of  15-month  time  series  from  150  m  depth  in  Storfjorden. 
Instrument  elevations  approximately  2  m  above  the  seabed. 
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Olgastretet:  USGS  Geoprobe 
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Figure  5.  11-month  time  series  of  velocity  from  Bear  Island  mooring 
Measurements  were  made  20  m  above  the  seabed  and  have  been  rotated 
to  represent  along  (v)  and  across  isobath  (u)  flow  directions. 
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Figure  6.  11-raonth  time  series  from  the  Bear  Island  mooring  2  showing 
suspended  sediment  concentration  20  m  above  the  seabed.  Record 
numbers  represent  hourly  measurements  beginning  on  Sept  15,  1993. 


sal,  %o  deg  C 


Figure  7.  11-month  time  series  from  the  Bear  Island  mooring  2  showing 
temperature  and  salinity  at  20  m  above  the  seabed.  Record  numbers 
represent  hourly  measurements  beginning  on  Sept  15,  1993. 
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Bear  Island  3,  cross-isobath  sediment  flux 
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Figure  8.  11-month  time  series  from  Bear  Island  mooring  3  showing 
suspended  sediment  flux  15  20  m  above  the  seabed  oriented  in 
along-isobath  and  cross-isobath  directions.  Record  numbers 
represent  hourly  measurements  beginning  on  Sept  15,  1993. 
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An  Overview  of  the  NRL  Large  Scale  Modeling  Effort  Studying 
the  Dispersion  of  Radioactive  Contaminants  in  the  Arctic 

Ruth  H.  Preller 
Presented  at  the  Workshop  on 
Modeling  the  Dispersion  of  Nuclear  Contaminants 
in  the  Arctic  Seas 

In  1993,  the  former  Soviet  Union  released  information  pertaining  to  the  dumping  of 
nuclear  waste  into  the  Arctic  Seas.  A  paper  known  as  the  Yablokov  report  or  the  "White 
Paper"  was  published  containing  data  on  both  source  locations  and  the  amount  and  type  of 
material  dumped  at  these  locations.  This  report  indicated  that  the  sources  were  low-level 
solid  and  liquid  waste,  high-level  waste  in  the  form  of  nuclear  reactors,  fuel  rods  with  spent 
nuclear  fuel  and  even  discarded  vessels  containing  nuclear  waste.  In  addition,  the  Ob,  the 
Yenisei  and  the  Pechora  rivers  were  sited  as  sources  of  nuclear  contamination.  The  first 
figure  included  in  this  abstract  is  a  table  from  this  report  showing  the  radionuclide  budget  for 
the  Kara  and  Barents  seas. 

The  object  of  this  study  is  to  determine  what  might  happen  to  this  radioactive  material 
if  and  when  it  enters  the  water  column.  To  look  at  possible  answers  to  this  question,  we  have 
made  use  of  the  Naval  Research  Laboratory’s  coupled  ice  ocean  model.  This  model  consists 
of  the  Hibler  ice  model  coupled  to  the  Cox  ocean  model  and  uses  a  0.28  degree  grid 
resolution  (approximately  25  km).  The  model  extends  from  the  North  Pole  southward  to 
approximately  30  degrees  North  latitude.  For  the  initial  phase  of  this  study,  we  have 
decoupled  the  ice  from  the  ocean  model  and  will  restrict  our  simulations  to  the  ocean  model. 

Several  ten  year  numerical  model  simulations  have  been  run.  In  each  case,  the  model 
was  driven  by  an  annually  averaged  wind  field  from  1986  atmospheric  data  from  the  U.S. 
Navy's  global  model  and  initialized  from  Levitus  climatology.  The  following  five  test  cases 
are  presented  using  sources  described  in  the  Yablokov  report:  low-level  solid  and  liquid 
waste,  high  level  waste,  rivers  (Ob,  Yenisei  and  Pechora),  Sellafield  (Irish  Sea)  and  finally  a 
combination  of  Sellafield  and  the  rivers.  Data  from  a  recent  1992  joint  Russian-Norwegian 
cruise  indicated  that  levels  of  radiation  in  the  Kara  were  ranging  from  5-25  Bq/m**3. 

Models  results  indicated  that  when  the  low  level  sources  were  used,  the  resulting  levels  of 
radiation  were  too  weak.  When  the  high  level  sources  were  used,  the  resulting  levels  of 
radiation  were  too  strong.  When  river  outflow  was  used  as  a  source,  the  resulting  levels  of 
radiation  were  consistent  with  those  measured  in  the  Kara.  The  second  figure  included  with 
this  abstract  shows  the  model  results  at  the  end  of  the  tenth  year  for  the  case  using  rivers  as  a 
source.  Note  that  5  Bq/m**3  is  approximately  0.1  PCi/1  which  is  represented  by  the  pink 
color  contour  in  this  plot.  When  a  simulation  was  run  using  Sellafield  as  a  source,  the  ten 
year  result  agreed  well  with  observation  in  the  Norwegian,  Greenland  and  Barents  Seas,  but 
appeared  to  give  results  that  were  too  weak  in  the  Kara  Sea  (see  the  abstract  by  Kathy  Crane 
in  this  proceedings).  The  final  test  (see  final  figure)  combining  these  two  sources,  the  rivers 
and  Sellafield,  appeared  to  give  the  best  results  when  compared  to  observations. 


Future  work  will  include  a  more  detailed  comparison  to  data  as  it  become  available 
(see  the  abstract  by  Cheng  in  this  proceedings),  continued  modeling  of  the  rivers  as  a  source 
of  radioactivity  in  conjunction  with  the  river  modeling  project  also  presented  at  this  workshop 
(see  abstracts  by  Paluszkiewicz  and  Hibler),  and  modeling  the  possible  transport  of 
radioactive  material  via  sediment  trapped  in  sea  ice. 


Table  5.  Anthropogenic  Radionuclide  Budget 
of  the  Barents  and  Kara  Seas  Ecosystem,  1961-1990 
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Figure  1.  Table  5  from  the  Yablokov  report,  page  29. 


OCEAN  CURRENT/NUCLIDES  (15M) 


1986RIVER 


The  dispersion  of  the  river  outflow  radionuclides  at  15  meters  deep 
at  year  1 0.  The  colors  which  are  in  a  logrithm  scale  represent  the 
radionuclide  level  in  PCi/1.  The  arrows  are  the  ocean  currents 
plotted  at  every  fourth  point.. 


lides  in  PCi/l 


radionuclides  versus  the  measured  Cs-137  and  Sr-90. 
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The  Effect  of  River  Discharge  in  a  Large  Scale 

Ice/Ocean  Model 

Richard  Allard 
Naval  Research  Laboratory 
Stennis  Space  Center,  MS  39529 


Abstract 

The  effect  of  river  inflow  on  the  growth  and  melt  of  ice  along  the  Siberian  Shelf 
is  examined  with  PrPS2.0,  the  spherical  version  of  the  Hibler  model  with  a  .28  deg 
resolution  in  the  Arctic  Basin  coupled  to  the  Cox  ocean  model.  This  coupled  model 
encompasses  a  hemispheric  domain  from  near  30.0  to  the  North  Pole  and  includes 
most  sea  ice-covered  areas  such  as  the  Gulf  of  St.  Lawrence,  the  Yellow  Sea,  the 
Seas  of  Okhotsk  and  Japein  as  well  as  the  Arctic  and  sub- Arctic  areas.  The  model 
is  forced  by  the  Navy  Ocean  Global  Atmospheric  Prediction  System  (NOGAPS). 
A  50-year  (1937-1986)  data  set  of  volume  fluxes  from  9  major  rivers  emptying  into 
the  Arctic  Basin  (8  Russian  and  the  Mackenzie)  is  used.  Each  discharge  is  specified 
at  several  grid  points.  The  largest  river  discharges  (Lena,  Yenisei)  edways  show 
maxima  in  June.  The  rivers  are  specified  as  source  regions  of  0  ppt  salinity  and 
monthly  varying  runoff  temperatures.  Test  cases  examined  include  a  ”no-river” 
case  ,  mean  (1937-1986)  rivers  and  an  ’’extreme”  river  Ccise  in  which  the  maximum 
monthly  discharge  during  a  50-year  period  was  used. 

The  inclusion  of  river  discharge  produces  initially  northward  flowing  surface  cur¬ 
rents  with  peak  velocities  of  6  cm/sec  near  the  Lena  River  in  Jime.  As  discharge 
decreases  during  the  summer  and  autumn,  these  surface  currents  decrease  in  mag¬ 
nitude.  The  combination  of  warmer,  fresher  river  water  during  summer  promotes 
a  thinner  ice  cover  in  the  vicinity  of  the  rivers.  In  the  Laptev  Sea,  there  is  about 
1  meter  less  ice  by  late  August  when  compztred  to  the  test  case  without  rivers.  In 
autmim,  the  fresher  water  has  a  warmer  freezing  temperature,  cdlowing  thin  ice  to 
grow  quicker  when  aii  temperatures  axe  below  freezing.  The  extreme  river  case 
shows  a  quicker  ice  melt  in  late  spring  and  summer,  and  more  ice  growth  in  the 
Laptev  Sea  in  autumn  as  a  consequence  of  a  warmer  freezing  temperature.  The 
Arctic  interior  shows  ice  differences  from  0  to  20  cm,  when  rivers  are  included. 


RIVER  L  OCA  RONS  IN  PIPS2. 0  MODEL  DOMAIN 
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MONTHLY  ARCTIC  RIVER  DISCHARGE 
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Surface  ocean  currents  in  the  Laptev  and  Kara  Seas  a):  with  rivers 
b):  without  rivers. 
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Sea  Ice  Transport  of  Pollutants  in  the  Laptev  Sea 

About  900,000  km^  of  the  polar  pack  ice  is  transferee!  annually  from  the  Arctic  Basin  to 
the  North  Atlantic.  The  largest  portion  of  this  exported  ice  cover  is  created  by  the  large 
scale  divergence  within  the  ice  pack.  A  smaller  portion  of  the  ice  cover  originates  in  the 
marginal  seas,  either  by  fall  freezing  of  the  seasonally  ice  free  waters  or  by  wintertime 
advection  away  from  the  coast. 

The  main  objective  of  this  study  is  to  estimate  the  annual  production  of  ice  in  the  Laptev 
Sea  and  to  determine  its  ultimate  fate.  The  study  is  motivated  by  the  possibility  that  ice 
formed  in  the  Laptev  Sea  may  be  an  agent  for  the  long  range  transport  of  pollutants  such 
as  radionuclides. 

The  strongest  evidence  for  long  range,  non-diffusive,  pollution  transport  is  the  ubiqui¬ 
tousness  of  “dirty  ice”  throughout  the  eastern  Arctic,  e.g  Pfirman  et  al  (1990)  and  Nuern¬ 
berg  et  al  (1993).  Although  much  of  the  ice  is  discolored  by  biogenic  material,  there  is 
clear  evidence  of  fine  grained  lithogenic  material.  While  the  biogenic  material  grows 
throughout  the  arctic,  the  lithogenic  material  can  only  be  incorporated  into  ice  formed  in 
the  coastal  waters. 

In  this  study  we  attempt  to  characterize  the  mean  and  interannual  variability  of  ice  pro¬ 
duction  by  investigating  the  winter  production  and  subsequent  melt  of  ice  in  the  Laptev 
Sea  from  1979  through  1992.  The  general  approach  is  to  associate  pollution  transport 
with  the  net  exchange  of  ice  area  from  the  Laptev  Sea  to  the  perennial  ice  pack.  The  pri¬ 
mary  data  sets  supporting  the  study  are  ice  charts,  ice  motion,  and  geostrophic  wind. 
Unfortunately  we  have  no  direct  measurement  of  pollutants. 
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An  Oriented  Sea  Ice  Dynamics  Model  with 
Contamination  Incorporation,  Transport,  and  Release 

by  Max  D.  Coon 

Northwest  Research  Associates,  Inc. 
and 

Robert  S.  Pritchard 
IceCasting,  Inc. 

A  model  is  being  developed  that  will  account  for  the  essential  physical  pix>cesses  of  sea  ice 
as  a  pathway  for  the  transport  of  radionuclides;  to  characterize  incorporation  of  radionuclides  into 
a  ice  cover,  advection  of  the  contaminated  ice,  and  dispersal  of  the  material  during  transport  and 
melting.  This  model  will  incorporate  a  new  anisotropic  constitutive  law  for  sea  ice  which 
accounts  for  lead  dynamics  directly  (See  Coon  et  al,  1992). 

Viewgraph  1  shows  the  mechanism  to  be  considered  in  the  model. 

Viewgraph  2  lists  mechanics  of  radioactive  waste  incorporated  into  sea  ice. 

Viewgraph  3  shows  radionuclides  released  into  the  nearshore  through  river  outflow.  The 
fresh  water  from  the  river  is  in  contact  with  the  floating  ice  cover.  Radionuclides  in  the  river 
water  are  incorporated  within  the  ice  cover  as  additional  freezing  takes  place. 

Viewgraph  4  lists  features  of  an  oriented  ice  thickness  distribution.  In  this  ice  model, 
which  accounts  for  the  dynamics  of  a  lead,  ice  redistribution  occurs  in  the  leads  which  results  in 
an  oriented  ice  thickness  distribution. 

Viewgraph  5  shows  a  sea  ice  element  with  a  lead  at  angle  0  to  the  x  axis  and  a  local 
coordinate  system  fixed  to  the  lead. 

Viewgraph  6  shows  a  lead  along  the  x  axis  in  physical  space  and  the  failure  surface  for  the 
lead  in  stress  space.  The  failure  surface  for  this  model  depends  on  direction,  which  is  required  to 
describe  the  ice  dynamics  on  the  physical  scale  of  a  few  to  tens  of  kilometers. 

Viewgraph  7  shows  a  refrozen  lead  containing  three  ice  categories  (including  open  water 
and  nilas).  The  three  failure  surfaces  shown  represent  three  lead  ice  thickness  categories. 

Viewgraph  8  lists  the  assumptions  about  the  ice  thickness  distribution. 

Viewgraph  9  illustrates  how  the  assumptions  of  Viewgraph  7  lead  to  equations  for  the 
growth  and  redistribution  of  ice  in  the  thickness  categories  for  both  lead  opening  and  closing. 

Viewgraph  10  lists  radioactive  waste  release  mechanisms  from  ice  to  water. 

Viewgraph  1 1  shows  that  the  goal  of  NWRA  and  ICI  is  to  provide  Ruth  Preller  of  NRL 
with  our  model.  Our  methodology  will  include  interaction  with  other  P.Ls  from  ANW AP. 


Coon,  M.D.,  G.S.  Knoke,  and  D.C.  Echert  (1992),  "Anisotropic  Pack  Ice  Constitutive  Law"  in  proceedings 
of  lAHR  11th  International  Ice  Symposium,  Banff,  Alberta,  June  15-19. 


Mechanisms  to  consider  in  the  model 
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Incorporation  of  rad.  waste  into  sea  ice 
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a)  Four  Ice-Thickness  Categories 


Ice  thickness  categories  and  their  failure  surfaces. 


Assumption  about  Ice  Thickness  Distribution 

•  The  mathematical  equations  and  calculation  procedures  are  being  developed  for 
the  ice  thickness  distribution  with  four  ice  thickness  categories,  based  on  Reimber, 
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When  tlie  j-lli  lead  is  opening,  i.e.,  for  5Vj  >  0,  llie  following  equations  describe  the  rale  of  change 

of  area  fraction  of  the  four  ice  categories. 

(10;  LS5V;  dT: 
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Release  of  rad.  waste  from  sea  ice  to  water 
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1  A  numerical  study  of  radionuclide  dispersal  in  the 
Kara  Sea 
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Santa  Cruz  CA  95064 

Fax:  408-459-4882 
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email:  dcsiv@cascade.ucsc.edu 

A  number  of  physical  processes  may  be  important  in  the  transport  of  radionuclides  from 
Russian  rivers  into  the  Arctic  Ocean.  These  include  transport  by  ocean  currents,  sea  ice 
and  geological  pathways.  The  amount  of  radionuclide  transported  by  each  of  these  paths  is 
dependent  on  time  of  year,  as  river  flow  rates,  nearshore  ice  fields  and  atmospheric  mixing 
of  the  ocean  change  seasonally.  The  long  term  goals  of  this  research  are  to  understand  how 
river  borne  contaminants  are  partioned  between  sea-ice,  sea  water  at  various  depths  and 
sediment  on  the  shelves?  The  emphasis  of  this  study  is  on  the  Kara  Sea  and  in  the  inflow 
from  the  adjacent  Ob  and  Yenesei  river  systems. 

A  regional  ocean  circulation  model,  presently  in  use  in  numerous  ocean  circulation  studies, 
is  being  adapted  for  application  to  the  Kara  Sea.  The  model  will  be  used  in  a  study  of  the 
various  processes  that  may  contribute  to  the  redistribution  of  radionuclides  in  the  Kara  Sea. 
The  important  modifications  of  the  coastal  ocean  model  for  Kara  Sea  application  involve 
incorporating  ice  and  radionuclide  loss  terms. 

Presently  we  are  forcing  the  model  with  a  coastal  buoyancy  source  over  simple  shelf 
topography.  The  magnitude  of  the  buoyancy  forcing  will  be  adjusted  to  correspond  to  Ob 
and  Yenesei  flow  rates  for  different  seasons.  In  subsequent  simulations  the  added  effects 
of  winds  will  then  be  considered.  Initial  simulations  will  not  involve  existing  ocean  flow 
fields.  It  is  possible  to  incorporate  these  in  subsequent  simulations  by  relaxing  offshore 
density  fields  back  to  a  prescribed  ocean  climatology.  After  the  sensitivity  of  river  plume 
distributions  to  the  above  factors  is  understood,  we  will  repeat  simulations  with  bathymetry 
fields  appropriate  to  the  Arctic  shelf  adjacent  to  the  Ob  and  Yenesei  rivers. 

It  will  not  be  feasible  to  examine  the  transport  of  radionuclides  in  all  the  possible  paths 
discussed  above.  We  will  however  consider  some  aspects  of  these.  It  will  be  possible  to 
include  the  dynamical  (but  not  thermodynamical)  effects  of  ice.  In  nearshore  regions  for 
example,  the  inclusion  of  free  drift  ice  cover  can  be  readily  implimented  to  examine  ice  edge 
ocean  transport  phenomena.  The  effects  of  pack  ice  may  be  incorporated  by  the  presense 
of  lateral  boundaries  which  simply  restrict  river  flow.  Incorporation  of  convective  forcing, 
associated  with  ice  formation  is  illustrated  below. 

As  a  first  step,  we  have  been  working  to  identify  space  and  times  scales  associated  with 
river  forced  flow  in  the  Kara  Sea.  These  determine  model  resolution  and  duration  of  simula¬ 
tions.  At  present,  it  appears  that  model  resolution  of  order  several  kilometers  and  seasonal 
time  scales  are  important.  At  these  resolutions  we  have  been  able  to  reproduce  the  results 
of  other  investigators  using  the  same  model  to  study  near  shore  buoyancy  forced  flows. 
Because  this  effort  is  new  this  year,  no  major  results  are  yet  available.  Several  preliminary 


simulations  are  shown  in  figure  1  and  contrast  the  different  responses  of  the  coastal  ocean 
to  dense  versus  fresh  coastal  ocean  buoyancy  forcing.  Bottom  topography  is  represented  by 
a  simple  linear  sloping  bottom  which  slopes  away  from  the  coast,  ranging  in  depth  from  20 
to  100  m.  Experiment  1  is  forced  with  dense  water  source  at  the  surface  adjacent  to  the 
coast  and  is  representative  of  the  flow  forced  in  a  nearshore  polynya,  where  ice  formation  and 
brine  rejection  are  occurring.  Experiment  2  is  forced  with  a  fresh  water  source  at  the  surface 
representative  of  river  inflow  into  the  coastal  ocean  in  the  absence  of  ice.  Both  processes 
may  occur  in  the  Kara  Sea  depending  on  season.  The  physical  processes  occurring  in  the  two 
experiments  are  quite  different.  In  experiment  1,  convective  instability  leads  to  baroclinically 
unstable  wave  growth  and  breakup  of  the  flow  into  eddies.  The  sloping  bottom  enhances 
the  unstable  wave  growth.  The  river  plume  in  experiment  2  is  surface  trapped  and  evolves 
stably  as  a  coastal  Kelvin  wave. 

The  Kara  Sea  model  effort  interfaces  with  the  research  of  a  number  of  other  investigators. 
Researchers  at  the  University  of  Washington  are  in  the  process  of  obtaining  Russian  data  sets 
for  the  Kara  Sea,  which  will  provide  valuable  input  for  model  initialization  and  verification. 
Data  sets  taken  by  Jim  Brooks  at  TAMU  will  also  provide  information  about  the  distribution 
of  river  plumes  near  the  mouth  of  the  Ob  and  Yenesei  rivers.  Elow  rates  and  geological 
transport  processes  in  the  Ob  and  Yenesei  river  -estuary  system  are  being  studied  by  the 
group  at  the  Pacific  Northwest  Laboratory.  Input  to  the  Kara  Sea  modeling  will  be  taken 
from  their  findings.  The  results  of  regional  Kara  Sea  modeling  should  also  benefit  the 
larger  scale  Arctic  basin  modeling  being  done  at  NRL,  in  the  form  of  improved  assessment  of 
radionuclide  residence  time  in  the  Kara  Sea  and  improved  boundary  conditions  for  the  larger 
domain  model.  Finally  it  will  be  possible  to  assess  the  importance  of  tides  in  redistribution 
of  river  borne  contaminants.  This  will  be  done  in  collaboration  with  David  Brooks  at  TAMU 
who  will  perform  compatible  experiments  using  a  tidally  driven  model. 
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Experiment  1.  Dense  water  source  at  the  surf ace. Left  column  is 

salinity  distribution  at  the  surface,  right  column  is 
salinity  distribution  at  z=-20  m. 
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Experiment  2.  Fresh  water  source  at  the  surface.  Left  column  is 

salinity  distribution  at  the  surface,  right  column  is 
salinity  distribution  at  z=-20  m. 
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Multiple-River  Interactions  in  Buoyancy-Driven  Coastal  Currents 

David  A.  Brooks 
Department  of  Oceanography 
Texas  A&M  University 
College  Station,  TX  77843 

Coastd  and  shelf  circulation  patterns  are  often  strongly  influenced  by  fresh  water  outflows 
from  rivers.  The  fresh-water  buoyancy  injection  provides  a  potential  energy  source  that, 
due  to  the  earth's  rotation,  eventually  sets  up  a  nearshore  current  that  is  trapped  against  the 
coast  in  the  Kelvin-wave  sense  (right  bounded  in  the  Northern  Hemisphere).  Numerical 
models  and  observations  indicate  that  near  the  river  mouth  the  freshened  water  first  forms  a 
bolus  or  bulge  of  low  salinity  water  that  spreads  in  both  directions  along  the  coast.  The 
buoyancy  driven  coastal  current  eventually  escapes  the  bulge  and  spreads  along  the  coast. 

In  a  coastal  region  influenced  by  multiple  rivers,  complex  interactions  can  occur.  When  the 
river  mouths  are  nearby,  the  bulges  from  the  adjacent  rivers  can  merge  and  act  essentially 
as  a  single  source  for  a  coastal  current.  On  the  other  hand,  when  the  mouths  are  separated 
sufficiently  that,  for  prescibed  flow  rates,  the  bulges  remain  distinct,  the  'upstream'  coastal 
current  can  be  redirected  offshore  as  it  encounters  the  strong  thermohaline  fronts  associated 
with  the  'downstream'  river  bulges.  Such  a  case  occurs,  for  example,  in  the  Gulf  of  Maine, 
where  the  combined  plume  from  the  Kennebec  and  Androscoggin  Rivers  intermpts  the 
Maine  coastal  current  that  is  set  up  in  part  by  outflow  from  the  Penobscot  River  farther  to 
the  east  (see  Figure).  Model  results  show  that  the  Kennebec- Androscoggin  plume  partially 
deflects  the  coastal  current  offshore  into  the  central  Gulf  of  Maine,  where  the  nutrient  flux 
influences  the  annual  cycle  of  productivity. 

Multiple  rivers  are  common  on  continental  shelves,  so  a  similar  effect  can  be  expected  in 
many  regions.  For  example,  the  Ob  and  Yenisee  Rivers  in  the  Kara  Sea  of  the  Arctic  Ocean 
have  very  large  run-off  in  spring  months  and  no  doubt  interact  in  complex  ways.  Pollutant 
pathways  in  the  Kara  Sea  and  toward  the  Alaska  shelf  in  the  associated  coastal  current 
therefore  depend  to  some  extent  on  the  way  in  which  the  individual  river  plumes  interact. 
Model  studies  are  planned  to  begin  examining  the  ways  in  which  these  interactions  occur 
and  how  they  depend  on  river  flow  rates  and  timing  of  maximum  volume  fluxes. 


Surface  velocity  and  salinity,  t=20  days,  June,  1982 


FIGURE  CAPTION 

Model  tidal-residual  surface  current  and  salinity,  showing  the  interaction  of  river  plumes  20 
days  after  the  start  of  a  run  intialized  with  hydrographic  data  from  June,  1982.  The  region 
shown  is  a  coastal  subregion  of  the  full  modeled  domain,  which  included  the  Gulf  of 
Maine,  Bay  of  Fundy  and  Georges  Bank  on  the  continental  shelf.  The  mouths  of  the 
Penobscot  (P)  and  Kennebec- Androscoggin  (KA)  Rivers  are  shown  by  the  filled  circles. 
The  contour  interval  is  5  ppt  for  salinity  <30  ppt  and  0.5  ppt  for  salinity  >30  ppt.  Dotted 
contours  show  detail  between  25  and  30  ppt  with  1  ppt  interval.  Vector  speed  scale  (cm  s‘ 
1)  is  indicated.  The  horizontal  arrowheads  identify  vertical  sections  not  shown  here  (see 
Brooks,  1994,  for  additional  details). 

Reference: 

Brooks,  D.A.,  1994:  A  Model  Study  of  the  Buoyancy-Driven  Circulation  in  the  Gulf  of 
Maine.  J.  Phys.  Oceanogr.,  24(11),  2387-2412. 
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S.J.  Timms’,  N.M.  Lynn’,  M.E.  Mount^,  Y.  Sivintsev^ 
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ABSTRACT 


In  progressing  its  work  for  the  International  Arctic  Seas  Assessment  Project  (lASAP),  under  the 
auspices  of  the  International  Atomic  Energy  Agency  (IAEA),  the  Source  Term  Working  Group 
developed  a  simple  spreadsheet  model  to  predict  a  radiation  release  profile  into  the  Kara  Sea 
from  reactor  fuel  and  activated  components  dumped  within  the  reactor  compartment  of  the 
icebreaker  Lenin. 

The  model  accounts  for  the  degradation  of  containment  materials  through  corrosion  and  other 
mechanisms,  and  predicts  annual  release  rates  to  four  thousand  five  hundred  years  into  the 
future.  Five  scenarios  are  modelled,  depending  on  the  integrity  of  the  levels  of  containment, 
a  major  factor  being  the  very  slow  rate  of  corrosion  of  the  UOj  fuel.  At  the  turn  of  the  century! 
total  release  rates  are  shown  to  vary  between  <0.01  to  50  Ci/yr  (0.37  GBq  to  1.85  TBq/yr) 
depending  on  the  model.  By  2050,  the  best  of  the  scenarios  are  predicting  release  rates  on 
the  order  of  0.6  Ci/yr  (22  GBq/yr)  and  at  2500,  the  rate  is  down  to  0.2  Ci/yr  (7.4  GBq/yr), 
dropping  by  a  decade  until  the  year  4000;  after  this,  the  radioactive  components  will  all  have 
been  released  from  their  containment  or  decayed  away. 

The  model  is  being  applied  to  other  forms  of  solid  radioactive  waste  dumped  in  the  Kara  and 
Barents  Seas,  in  a  programme  of  work  for  the  lASAP,  aimed  at  assessing  the  collective  release 
profile  from  all  significant  dump  sites  in  this  region. 

INTRODUCTION 

Background 


1.  In  the  Spring  of  1993,  an  English  translation  of  the  Russian  report.  Facts  and  Problems 
Related  to  Radioactive  Waste  Disposal  in  Seas  Adjacent  to  the  Territory  of  the  Russian 
Federation was  released.  The  findings  presented  in  this  report  were  the  result  of  a  scientific 
study  commissioned  in  October  1992  by  the  Office  of  the  President  of  the  Russian  Federation. 
The  White  Book,  as  the  report  was  later  called,  reported  that:  (1)  between  1965  and  1988,  16 
naval  reactors  from  seven  former  Soviet  Union  submarines  and  the  icebreaker  Lenin,  each  of 
which  suffered  some  form  of  reactor  accident,  were  dumped  at  five  sites  in  the  Kara  Sea; 


radioactive  waste  was  discharged  at  sites  in  the 

leJr^  between  1964  and  1991,  low  and  intermediate 

level  solid  radioactive  waste  was  dumped  at  sites  in  the  Barents  and  Kara  Seas. 

fn  discarded  naval  reactors,  six  of  the  16  contained  their  spent  nuclear  fuel  fSNFI 

I? t''®  "1^®®  '■e™  naval  readore  wS 
disposed  of  in  a  reinforced  concrete  and  stainless  steel  shell  container.  The  vast  maioritv  of 
the  low  and  intermediate  level  solid  radioactive  waste  was  disposed  of  in  containprc  nT.  u  ^ 

navairaSs  and““^^ 

level  solid  radioactive  waste  varied  in  depth  from  12  to  380  m. 

T^f  reported  estimates  of  total  radioactivity  at  the  time  of  disoosal  th^« 

were  (1)  2,300  kCi  (85  PBq)  of  fission  products  in  the  SNF;  (2)  100  kCi  (O.fpBq)  of  ^Co  in 
he  reactor  cornponents.  (3)  24  kCi  (0.89  PBq)  of  unspecified  origin  in  the  low  level  iTquid 
pfin"?  n?'"®  waste,  over  50%  of  which  was  discharged  in  the  Barents  Sea;  and  (4)  16  kCi  0  59 
PBq)  of  unspecified  origin  in  the  low  and  intermediate  level  solid  radioactive  waste  over  95% 
Of  Which  was  discarded  ih  the  Kara  Sea,’  With  rare  exception,  no  radlonucllderw?re  M^d 
and  no  estimate  provided  for  the  current  levels  of  radioactivity. 

International  Arctic  Seas  Assessment  Project 

thP<5 concern  over  the  possible  health  and  environmental  effects  from  disposal  of 
Inttrn5  wastes  in  the  shallow  waters  of  the  Arctic  Seas  prompted  the 

responsibilities  to^he  L^don 

lAciAP  il  ^  International  Arctic  Seas  Assessment  Project  (lASAP)  ^  The 

or  Actual  and  Potential  Consequences  of  Dumping  Radioactive  Waste  into  Arctic  Seas. 

5^  The  stated  objectives  of  the  lASAP  are  to:  (1)  assess  the  risks  to  human  health  and  to  the 
nvironment  associated  with  the  radioactive  waste  dumped  in  the  Kara  and  Barents  Seas  and 
thev  possible  remedial  actions  related  to  the  dumped  wastes  and  to  advise  on  whether 

es®abnsLd^7i?T  1°  objectives,  four  working  groups  were 

established.  (1)  Source  Term,  (2)  Existing  Environmental  Concentrations^  f3)  Transfer 

echanisms  and  Models,  and  (4)  Impact  Assessment  and  Remedial  Measures. 

Source  Term  Working  Group 

the  Graw  was°he1d'i"lL^il1ir®-°^  ''®^®after  identified  as 

tne  uroup,  was  held  at  the  IAEA,  in  Vienna,  Austria  on  January  10-14,  1994  The  purposes 

aliessmenr  f  progress  and  needs  on  the  source  term  evaluation  for  the 

dnmnf^  1  environment  associated  with  nuclear  reactors 

Trl^  'f  prepare  a  preliminary  source  term  description  to  be  used  by 

the  Transfer  Mechanisms  and  Models  Working  Group.  ^ 

7.  During  the  meeting,  the  Group  reviewed  and  discussed  source  term  reports  prepared  bv 
consulting  members,  one  from  the  United  States  and  one  from  Russia  Th^es^ate  p^^^^^^^^ 
by  the  consulting  member  from  the  United  States  was  for  the  inventory  of  selected  acSes 
and  fission  products  in  the  SNF  and  activation  products  in  the  reactor  components  fnd  p  ^ 


system  corrosion  products  of  the  16  naval  reactors.^  Estimates  for  the  Lenin  radionuclide 
inventory  were  based  upon  core  history  information  obtained  from  translations  of  Russian 
articles  available  in  the  open  literature,  while  the  estimates  of  the  nuclear  submarine 
radionuclide  inventories  were  based  upon  core  history  information  developed  from  open 
literature  Western  estimates  of  their  operating  characteristics.  The  estimate  prepared  by  the 
consulting  member  from  Russia  was  for  the  inventory  of  selected  actinides  and  fission  products 
in  the  SNF  and  activation  products  in  the  reactor  components  of  the  Lenin."*  In  addition  to  the 
inventory  estimate,  this  report  provided  detailed  information  on  the  operating  history  of  the 
Lenin:  condition  of  the  reactor  compartment,  reactor  pressure  vessels  (RPV),  and  the  SNF  prior 
to  disposal,  and  the  disposal  operations.  From  the  standpoint  of  the  Group,  the  later  of  these 
two  source  term  estimates  was  definitive  for  the  Lenin. 

8.  The  Group  was  also  provided  with  preliminary  information  from  the  September  -  October 
1993  joint  Norwegian/Russian  expedition  to  the  Tsivolka  Inlet,  Stepovoy  Inlet,  and  Novaya 
Zemlya  Depression  disposal  sites  in  the  Kara  Sea.^  Included  among  this  information  were 
preliminary  results  from  underwater  camera  investigations  of  specific  objects  and  in  situ  gamma 
spectroscopy  of  bottom  sediments. 

9.  As  the  information  from  the  joint  Norwegian/Russian  expedition  was  preliminary  at  best, 
and  because  Russian  estimates  for  the  radionuclide  inventories  in  the  nuclear  submarine  naval 
reactors  would  not  be  available  for  evaluation  before  the  Fall  of  1 994,  the  Group  focused  its 
efforts  on  the  Lenin.  The  following  section  summarizes  the  information  pertinent  to 
development  of  the  Lenin  source  term. 

LENIN  SOURCE  TERM* 


Background 

10.  Launched  in  Leningrad  in  1959,  the  icebreaker  Lenin  was  the  first  nuclear  merchant  ship 
in  the  world.  During  her  31  years  in  commission,  the  Lenin  had  two  separate  steam  generating 
plants  (SGPs).  The  first  SGP  contained  three  pressurized  water  reactors  (PWRs)  of  90  MW 
maximum  thermal  power  each  and  operated  from  1959  to  February  1965,  when  during  routine 
repair  of  the  SGP,  an  operator  error  allowed  the  core  of  the  N2  PWR  to  be  left  without  water 
for  some  period  of  time.  As  a  consequence,  a  part  of  the  reactor  core  was  damaged  due  to 
residual  heat.  It  is  the  characteristics  of  this  first  SGP  that  form  the  basis  for  the  Lenin  source 
term. 

Characteristics  of  the  First  SGP 

11.  The  entire  SGP,  including  the  steam  generators  and  pumps,  was  located  in  the  centre  of 
the  ship  in  an  isolated  reactor  compartment.  The  three  PWRs  were  aligned  vertically  in  a  plane 
perpendicular  to  the  keel  and  were  surrounded  by  a  large  iron-water  tank  shield.  The  biological 
shield  located  above  the  three  PWRs  was  made  of  limonite  concrete  and  a  heat  resistant 
composition  of  graphite  and  boron. 

12.  Each  cylindrical  RPV  was  made  from  low  alloy  steel,  with  approximate  dimensions:  2-m 
diameter,  5-m  height,  and  40-mm  thick  wails  with  a  5-mm  internal  stainless  steel  (SS)  layer  of 
cladding.  The  height  and  diameter  of  each  PWR  core  was  approximately  1.6  m  and  1  m, 


respectively,  with  each  core  containing  219  cylindrical  technical  fuel  channels  (TFCs)  with  36 
pin-type  fuel  elernents.  The  fuel  elements  contained  fuel  pellets  made  from  sintered  UO 
enriched  to  5.0%  and  clad  with  a  Zr-Nb  alloy  or  SS.  ^ 


1 3.  To  reduce  heat  and  radiation  effects  on  each  RPV  and  subsequently  extend  their  operatino 
lives,  thermal  shields  consisting  of  five  concentric  SS  cylinders  of  20-  to  30-mm  thickness  each 
(185  mm  total)  were  positioned  around  each  reactor  core.  Two  of  these  five  concentric 
cylinders  were  welded  to  the  upper  and  lower  rings,  which  connected  with  the  upper  and  lower 
plates  used  for  fixation  of  the  TFCs.  ^ 


14.  The  first  fuel  load  lasted  from  1959  to  1962  and  consisted  of  80  kg  of  in  each  reactor 
core.  The  integrated  power  productions  with  this  first  fuel  load  were  equal  to  17.8  GWd  for  the 
N1  PWR  and  18.0  GWd  for  the  N2  and  N3  PWRs.  Refuelling  occurred  in  1963  The  second 
fuel  load  lasted  from  1963  to  1965  and  consisted  of  129  kg  clad  in  SS  in  the  N1  PWR  and 
75  kg  U  clad  in  Zr-Nb  in  the  N2  and  N3  PWRs.  The  N2  PWR  operated  from  July  19  1963 
to  November  17,  1963  and  from  June  22,  1964  to  November  13,  1964  for  a  total  period  of  267 
days.  The  mean  power  of  the  N2  PWR  during  the  second  fuel  load  was  53  MW.  At  the  time 
of  the  reactor  accident,  the  integrated  power  production  for  the  N2  PWR  was  approximately 
14.2  GWd  thermal  and  the  burnup  was  equal  to  18.9  MWd/metric  tonne  initial  heavy  metal. 
Integrated  power  productions  for  the  second  fuel  load  were  equal  to  22  5  GWd  for  the  N1  PWR 
and  17.5  GWd  for  the  N3  PWR. 


15.  Initially,  all  undamaged  and  damaged  SNF  was  removed  from  the  three  RPVs  However 
as  a  consequence  of  the  accident,  only  94  of  the  21 9  TFCs  from  the  N2  RPV  could  be  removed 
for  normal  disposal.  The  remaining  125  TFCs  were  removed  with  three  of  the  concentric  SS 
cylinders  of  the  thermal  shield  and  placed  into  a  cylindrical  reinforced  concrete  and  SS  shell 
container  of  inner  and  outer  diameter  of  approximately  1  m  and  2  m,  respectively.  This 
reinforced  concrete  and  SS  shell  container,  hereafter  identified  as  the  box,  was  filled  with  a 
hardening  compound  based  on  furfural  resin  and  then  stored  on  land  for  approximately  two 


16.  Before  disposal,  the  primary  circuit  loops  and  equipment  were  washed,  dried,  and 
hermetically  sealed;  the  most  radioactive  inner  parts  of  each  RPV  were  covered  with  the  furfural 
based  hardening  compound;  the  ceiling  of  the  reactor  compartment  was  equipped  with  special 
pressure  relief  valves;  the  box  with  the  SNF  was  transferred  from  its  land  storage  location  to 
e  top  of  the  SGP  steel  shielding;  and  the  hole  in  the  reactor  compartment  ceiling  where  the 
box  was  inserted  was  covered  with  a  10-mm  thick  carbon  steel  plate  and  sealed  by  weld. 

1 7  In  August  1 967,  the  Lenin  reactor  compartment  with  her  three  RPVs  and  SNF  was  dumped 
in  the  shallow  water  of  Tsivolka  Inlet,  as  shown  in  Figure  1,®  at  an  estimated  depth  of  50  m 
irectly  from  the  Lenin  through  the  bottom  of  her  hull.  A  diagrammatic  representation  of  the 
entire  dumped  assembly  appears  in  Figure  2.* 


Radionuclide  Inventories 


Fission  Products  and  Actinides 

1 8.  Calculations  of  core  cell  burnup  in  the  N2  PWR  were  performed  by  the  spectral  code 
CETERA.®  The  fission  product  and  actinide  activities  were  estimated  using  the  RECOL^  library 
data  base,  which  was  generated  on  the  bases  of  the  latest  versions  of  the  evaluated  nuclear 
data  files,  ENDF/B-V,  with  corrections  based  on  the  results  of  critical  experiments.®  The 
criticality  problem  was  solved  for  a  realistic  3-D  geometry  model  of  a  TFC  by  Monte-Carlo  with 
RECOL  and  checked  with  MCNP  for  fresh  fuel  load.  One-group  cross  sections  were  prepared 
for  burnup  calculation  of  critical  loads  of  both  fresh  and  spent  fuel  and  input  to  ORIGEN-2’°  for 
detailed  radionuclide  content  calculations. 

Activation  Products 

19.  All  internal  reactor  constructions  were  made  from  SS;  the  RPVs  were  made  from  low  alloy 
steel.  Almost  all  activation  products  associated  with  the  Lenin  PWRs  originated  in  the  internal 
constructions  made  from  SS. 

20.  Analysis  of  neutron  activation  shows  that  only  four  radionuclides  are  of  consequence  at 
10  or  more  years  after  reactor  shutdown.  These  are  ®®Ni,  ®°Co,  ®®Ni,  and  Activities  of  ®®Ni 
and  ®°Co  were  estimated  from  the  following  equation; 

A  =  3.1x10^®  (P„„J  [p,  +  p,  +AK  (Z,  /  Z,.J]  TAH>, 

where  3.1x10’®  is  the  number  of  fissions/second-MW;  P,„„„  is  the  mean  thermal  power  of  the 
reactor,  MW;  p^  is  the  fraction  of  the  neutron  flux  captured  in  the  main  reactivity  compensating 
lattices  per  unit  fission;  p,  is  the  fraction  of  the  neutron  flux  captured  in  the  additional  reactivity 
compensating  lattices  per  unit  fission;  AK  is  the  fraction  of  the  neutron  flux  captured  in  the 
shielding  assemblage  and  lower  plate  per  unit  fission;  Z,  is  the  macroscopic  capture  cross 
section  for  thermal  and  intermediate  neutrons  in  the  steel  of  the  shielding  assemblage,  in  cm’’; 
Zj*,„  is  the  macroscopic  capture  cross  section  for  thermal  and  intermediate  neutrons  in  the  steel 
and  water  of  the  shielding  assemblage,  in  cm  ’;  T  is  the  irradiation  time,  in  days;  A  is  the 
radioactive  decay  constant  for  ®®Ni  or  ®°Co,  in  days'’;  and  'P  is  the  ratio  of  the  activation  cross 
section  of  Ni  or  Co  to  the  total  activation  cross  section  of  the  steel  under  irradiation. 

Activities  for  Ni  and  C  were  estimated  from  a  British  calculation  of  activation  products  in  the 
reactor  components  of  a  generic  nuclear  powered  submarine  one  year  after  shutdown.” 

^  presents  a  summary  of  the  estimated  radioactivity  in  the  box  that  is  sealed  within  the 
reactor  compartment  of  the  icebreaker  Lenin  for  1993.®  “” 

Table  2  presents  a  summary  of  the  estimated  radioactivity  for  activation  products  in  the  three 
reactor  pressure  vessels  that  are  sealed  within  the  reactor  compartment  of  the  icebreaker  Lenin 
for  1993.“” 


LENIN  RELEASE  MODEL 


'  V  '^3s  developed  to  provide  a  workable  estimate  of  the  rate  of  release 

of  activity  into  the  environment.  With  so  many  unknowns,  especially  with  respect  t^S 
behaviour  of  oxide  fuel  in  sea  water,  corrosion  rates,  and  the  effectiveness  of  the  various 

in  ^  assumptions  were  necessary  as  described  below 

in  al  cases,  such  assumptions  are  pessimistic  so  as  to  over  rather  than  under-estimate  the  rate 
of  release  of  activity.  The  model  has  the  ability  to  vary  most  of  the  operating  parameters  for 

22.  All  the  fission  products  and  actinides  are  in  the  box  placed  on  top  of  the  SGP  This  box 

also  contains  three  thermal  shields  from  the  N2  RPV,  representing  an  estimated  27%  of  the 

total  inventory  of  activation  products.  The  remaining  73%  of  activation  product  are  In  the 

O fToy  u  ^  further  assumed  that  the  activation  products  in  the  RPVs  are  distributed  in  the  ratio 
80%  thermal  shields  10%  RPV  cladding,  10%  RPV  walls.  a'siriouted  in  the  ratio 

23.  To  cope  with  the  different  forms  of  radioactive  material  and  their  containment 
arrangements,  the  model  is  based  upon  the  four  cases  depicted  in  Figure  3  with  each  case 

^Xh^rTH-  sof  containment.  Case  I  assumes  no  containment  whatsoever 

with  the  radioactive  material  directly  in  contact  with  the  sea.  Cases  II  to  IV  progressively  add 
containment  barriers  that  have  an  inhibiting  effect  on  the  expected  release  profile  from  the 

'■  ^  ^  assumed  that  adds  a  degree  of  realism  by 

considering  the  collapse  of  support  for  the  box  and  consequent  displacement  of  its  lid. 

process  of  release  of  radionuclides  from  the  activated  RPV  walls  and  SS  thermal 
shields  IS  deemed  to  be  by  corrosion.  For  fission  products  and  actinides  in  the  oxide  fuel,  once 
ea  water  penetrates  through  to  the  fuel,  the  release  process  is  assumed  to  be  immediate  from 
mfi.?  followed  by  a  very  slow  dissolution  of  fuel  grains.  Reasonably  accurate 

hppn  nLw release  into  the  sea  are  clearly  crucial  for  realistic  modelling,  and  have 
ntnarfrlfa  t  ^  detailed  survey  of  current  literature.  The  result  of  this  work  by  the 

rennrtpm  Technology,  Royal  Naval  College,  Greenwich,  has  been 

reported  separately  and  is  the  source  document  for  the  majority  of  rates  of  material 
decomposition  used  in  this  paper.  maienai 

The  constituent  elements  of  the  model  are  discussed  in  more  detail  under  case  headings: 

Case  I  -  No  Containment  Barriers 

SS  internal  clad  and  concentric  SS  thermal 

u®  pins:  and  remaining 

three  SS  thermal  shields  from  the  N2  RPV  are  directly  in  contact  with  the  sea. 


Activation  Release  from  the  RPV  Walls  and  Thermal  Shields 

26.  Based  on  published  data/  vessel  geometry  is  assumed  to  be  as  shown  in  Figure;  4. 

ss  ss 

RPV  CLAD  SHIELDS  -> 


Outside  RPV 

Reactor 

Compartment 


Inside  RPV 


Thicknesses  in  mm.  40  5  14.4  14.4 

Figure  4:  RPV  geometry  showing  RPV  wall , cladding  and 

the  internal  stainless  steel  thermal  shields 

0  immediately  from  the  outside  at  an  assumed  rate  of 

.  5  nirn/yr,  releasing  activation  products  into  the  surrounding  water.  Similarly  the  inner 

ayers  of  SS  cladding  and  SS  thermal  shields  corrode  into  the  sea,  but  at  a  rate  of  0  02  mm/yr” 
from  both  sides  simultaneously.  ^ 

Activation  Release  from  the  N2  PWR  Fuel  Rods  and  Thermal  Shields 

27.  As  for  the  SS  thermal  shields  in  the  RPVs,  the  three  from  the  N2  RPV  are  assumed  to 
corrode  away  on  both  sides  at  a  rate  of  0.02  mm/yr.”  The  damaged  fuel  from  the  N2  PWR  is 

SXallov  Ha!!!?"®®  4-5-mm  diameter.  10-g/cm^  density,  carried  in 

2r-Nb  alloy  cladding.  In  the  absence  of  details  of  the  extent  of  damage  to  the  125  fuel  pins 

does  cladding  material  has  been  fully  breached.^and 

does  not  constitute  a  barrier  to  fission  product  release  from  the  fuel. 

contact  with  sea  water,  an  immediate  release  of  20%  of  the 
inventory  is  assumed  from  fuel  grain  boundaries.  The  remainina  fission 

Sf  30x?0-4"/cr^4a^^^^^  are  released  through  dissolution  of  fuel  grains  at  a  pessimistic  rate 
of  30x10  g/cm  /day.  For  fuel  pins  with  a  4.5-mm  diameter  and  10-g/cm^  density  this  eouates 

inear  surface  for  the  circular  fuel  pin  and  an  evenly  spread  rate  of  dissolution.  Allowing  for 


the  geometry  of  the  pins,  activity  release  from  dissolution  will  commence  at  its  maximum  rate 
reducing  until  the  pins  finally  disappear. 

Model  results 


29.  Figure  6:  This  shows  the  initial  high  release  [36  kCi/yr  (1.3  PBQ/yr)]  due  to  the  20%  fission 
product  inventory  released  immediately  to  the  sea.  This  drops  to  56  Ci/yr  (2.1  TBq/yr)  by  the 
year  2000,  and  drops  further  after  2300  when  the  steels  have  largely  corroded  away.  From 
then  on,  the  only  release  is  due  to  the  corrosion  of  the  remaining  SNF.  By  the  year  3000  the 
release  rate  is  0.1 1  Ci/yr  (4.1  GBq/yr)  and  at  4000,  it  is  0.02  Ci/yr  (0.74  GBq/yr).  After  the  year 
4000,  the  SNF  has  disappeared. 


Case  II  -  Reactor  Compartment  Containment  Barrier 

30.  The  reactor  compartment  walls  are  assumed  to  form  a  boundary  around  the  radioactive 
waste  material,  and  to  thus  modify  the  release  into  the  surrounding  sea  by  the  containment 

effectiveness  factor  K^,  with  the  release  at  time  t,  in  years,  for  fission  and  activation  products 
given  by: 

Aj(t)  =  Aj[Case  l(t)]  x 

31.  As  shown  in  Figure  2,  pressure  release  valves  allowed  sea  water  into  the  reactor 
compartment  on  sinking.  Corrosion  is  therefore  assumed  to  occur  on  either  side  of  the  hull,  but 
at  a  reduced  rate  internally  because  of  low  oxygen  levels.  The  weakest  area  of  containment 
is  in  the  vicinity  of  the  weld  securing  the  mild  steel  cover  plate.  At  a  mean  corrosion  rate  of 
0.166  mm/yr,’^  on  both  sides,  the  estimated  life  of  an  assumed  10-mm  thick  weld  is  30  years. 
To  this  point,  the  effectiveness  of  the  containment  boundary  is  considered  to  be  100%  (Kc=0): 
thereafter,  it  is  assumed  to  be  90%  effective  (Kp=0.1).  General  corrosion  of  the  10-mm  thick 
mild  steel  hull  in  this  area  continues  at  a  rate  of  0.075  mm/yr^^  on  both  sides  of  the  plating,  until 
complete  failure  eliminates  the  effectiveness  of  this  boundary  (Kc=1)  67  years  from’  first 
immersion  in  sea  water  in  1967.  These  values  of  are  injected  into  the  model  as  step 
changes. 

Model  Results 

32.  Figure  7.  Adding  the  reactor  containment  barrier  to  Case  I  reduces  the  initial  release  to 
the  sea.  It  is  assumed  that  the  20%  soluble  fraction  escapes  in  1 996  [950  Ci/yr  (35  TBq/yr)] 
when  the  reactor  compartment  closure  weld  fails  and  K,.  =  0.1.  Once  the  weld  fails,  release 
rates  drop  to  below  10  Ci/yr  (370  GBq/yr)  and  then  rise  to  15  Ci/yr  (560  GBq/yr)  as  the  full 
containment  barrier  corrodes  away,  then  follow  the  same  pattern  as  Case  I. 


Case  III  -  Furfural  Containment  Barrier 

33.  Information  is  lacking  on  the  precise  composition  of  the  hardening  mixture  furfural,’  in 
which  the  radioactive  material  is  encased.  It  is  known  that  it  is  based  on  a  chemical  obtained 
from  the  acid  hydrolysis  of  the  polysaccharides  of  oat  husks,  corn  cobs  or  straw,  and  that  it  can 


be  polymerized  with  compounds  such  as  acetone  or  urea  to  form  a  hard  resin.  An  effective 
lifetime  of  500  years  is  quoted’  for  this  material.  However,  in  the  absence  of  reliable  data  on 
the  performance  of  furfural  in  such  environmental  conditions,  a  conservative  lifespan  of  100 
years  in  the  radiation  environment  is  assumed  in  the  model.  It  has  the  effect  of  inhibiting  the 
activity  released  from  ail  components  encased  in  furfural  by  the  factor  K^.  The  multiplying  factor 
K,  is  zero  until  the  Lenin  reactor  compartment  is  dumped,  and  then  ramps  to  unity  over  the 
ensuing  100  years.  The  breakdown  of  the  material  is  assumed  to  be  due  to  radiation  damage, 
cracking,  and  decomposition  in  sea  water. 

Model  Results 


34.  Figure  8;  Adding  the  protection  of  furfural  to  the  box  and  interior  of  the  RPVs  cuts  down 
the  release  rates  in  the  early  years.  The  activation  products  corroding  from  the  outside  of  the 
RPV  walls,  as  they  are  unprotected  by  the  resin,  are  still  released,  but  the  remainder  is 
encapsulated  in  the  resin.  The  form  of  the  release  rate  follows  Case  II  after  the  resin  loses  its 
complete  effectiveness  in  100  years,  but  the  initial  soluble  pulse  has  been  reduced  by  a  third 
to  290  Ci/yr  (IITBq/yr).  Choosing  a  lifetime  of  100  years  for  the  resin  means  that  its 
effectiveness  is  comparable  to  the  reactor  compartment  itself  (Case  II).  The  spreadsheet  model 
can  be  run  with  a  range  of  furfural  lifetimes  and  then  its  barrier  properties  can  be  shown  to 
extend  beyond  the  Case  II  protection. 

Case  IV  -  Box  Containment  Barrier 

35.  The  fuel  pins  and  three  SS  thermal  shields  from  the  N2  PWR  are  sealed  within  the  box 
with  welded  closure  lid,  and  the  RPV  heads  are  in  place.  For  the  latter,  as  bottom  and  top  RPV 
enclosures  are  only  bolted  into  position,  sea  water  is  assumed  to  have  penetrated  at  the  time 
of  disposal  and  thus  interfaced  with  the  furfural  to  commence  its  degradation. 

Corrosion  of  the  SS  Shell  Around  the  Box 

36.  The  box  geometry  is  shown  in  Figure  5.  Assuming  a  10-mm  thick  closure  weld  and  a 
mean  pitting  corrosion  penetration  rate  of  0.5  mm/yr,’^  initial  penetration  will  occur  after  20 
years,  admitting  sea  water  into  the  box  to  interface  with  the  furfural  and  start  its  degradation. 
With  a  general  bulk  corrosion  rate  of  0.02  mm/yr,”  complete  penetration  of  the  SS  shell  around 
the  concrete  can  be  expected  500  years  later. 
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encapsulated  in  furfural . 


Deterioration  of  the  Concrete  in  Sea  Water 


37.  There  is  little  data  on  the  performance  of  concrete  in  sea  water.  Available  literature’^ 
indicates  that  a  concrete  structure  such  as  this  should  have  a  long  life  if  it  does  not  suffer 
cracking  from  radioactive  heating  or  mechanical  damage.  However,  changes  in  the  concrete 
structure  over  time  are  likely,  and  will  alter  its  effectiveness  as  a  barrier  to  the  diffusion  of 
radioactive  species  and  to  the  inward  diffusion  of  salt  water  and  dissolved  oxygen.  The  model 
assumes  such  diffusion  commences  after  general  sea  water  penetration  of  the  clad  and 
increases  at  a  linear  rate  for  a  further  200  years  until  the  concrete  is  assumed  completely 
porous.  The  box  containment  factor ,  K^,,  therefore  has  a  value  of  zero  until  failure  of  the  box 
closure  weld  and  then  ramps  to  unity  over  the  subsequent  200  years.  This  is  less  than  the  time 
for  the  complete  degradation  of  the  SS  of  the  box  itself. 

Model  Results 

38.  Figure  9:  Although  the  outside  RPV  walls  always  contribute  to  the  release  rate  after  the 
reactor  compartment  walls  have  been  breached  (Case  II  and  Case  III),  the  box  is  shown  to 
reduce  the  overall  rate  to  below  1  Ci/yr  (37  GBq/yr).  The  early  activity  release  comes  from  the 
RPVs  as  the  furfural  degrades  (Case  III).  When  the  box  weld  fails,  additional  activity  will  be 
released.  The  concrete,  with  its  lifetime  of  200  years,  and  the  outer  SS  shell  will  inhibit  the 
release.  Overall  rates  do  not  exceed  0.7  Ci/yr  (26  GBq/yr)  by  the  year  2100.  The  concrete 

ceases  to  provide  protection  long  before  the  outer  SS  shell  corrodes  completely  just  before  the 
year  2500. 


Case  V  -  Box  Support  Collapse 

39.  Assuming  that  the  Lenin  reactor  compartment  has  remained  upright  on  the  sea  bed,  with 
the  box  securely  located,  this  final  case  considers  the  consequences  of  the  box  losing  its  lid 
through  displacement  following  failure  of  its  supporting  deck  and  hull  structure.  If  the  mild  steel 
IS  assumed  to  corrode  away  at  a  rate  of  0.075  mm/yr,’^  a  15-mm  support  framework,  corroding 
on  both  sides,  would  collapse  in  100  years  from  first  immersion  in  sea  water,  nearly  the  same 
time  of  the  complete  failure  of  the  surrounding  hull  structure  and  its  containment  (Case  I).  At 
this  point  the  massive  box  is  assumed  to  topple  and  lose  the  lid,  eliminating  its  effectiveness 
as  a  containment  barrier  (Kt,-1).  All  that  remains  to  inhibit  the  release  of  fission  products, 

X  activation  products  will  be  the  furfural,  and  with  the  lifetime  chosen  for  this  model 

of  100  years,  this  provides  no  barrier. 

40.  It  was  also  assumed  in  this  scenario,  that  the  20%  soluble  fraction  of  fission  products  will 
be  released  when  the  box  support  collapses. 

Model  results 


41.  Figure  10:  Before  the  box  collapse,  release  rates  are  identical  to  Case  IV.  When  the  hull 

frames  corrode  away  in  100  years,  the  box  lid  breaks  off  and  the  SNF  is  exposed  to  the  sea 
water,  by  now  without  the  protection  of  furfural  (and  the  box).  Figure  10  shows  the  initial  pulse 


of  the  20%  soluble  fraction  of  the  fission  product  inventory,  released  when  the  box  breaks  open 
in  2066,  at  a  rate  of  1900  Ci/yr  (70  TBq/yr)  to  the  sea.  Afterwards,  the  rate  drops  down  to  1.9 
Ci/yr  (70  GBq/yr)  at  2100,  0.5  Ci/yr  (19  GBq/yr)  by  2200,  and  then  mirrors  Case  I,  with  no 
containment. 


CONCLUSIONS  AND  FURTHER  WORK 

42.  All  the  scenarios  are  brought  together  in  Figure  1 1 ,  which  shows  that  Case  IV  provides 
the  best  level  of  protection,  as  might  be  expected.  With  all  the  pessimistic  assumptions 
(paragraph  21)  on  the  furfural  and  concrete  lifetimes,  release  rates  never  exceed  0.7  Ci/yr  (26 
GBq/yr),  and  stay  within  the  0.1  to  0.2  Ci/yr  (3.7  to  7.4  GBq/yr)  band  for  the  next  500  years. 

43.  The  box  collapse  scenario.  Case  V,  represents  the  greatest  risk,  allowing  the  soluble 
fraction  of  fission  products  to  be  washed  out  in  the  year  the  structure  fails  and  the  lid  breaks 
off.  This  case  is  most  sensitive  to  the  lifetime  assumed  for  furfural  protection. 

44.  Recommendations: 

a.  As  the  box  collapse  scenario  is  deemed  to  be  the  "worst"  of  the  scenarios  that  could 
be  envisaged  realistically,  it  is  suggested  that  in  the  next  exploratory  cruise  to  the  Kara  Sea, 
the  Lenin  reactor  compartment  be  located  in  Tsivolka  Inlet  and  the  state  of  the  structure  be 
photographed.  In  particular,  the  angle  of  repose  and  the  state  of  the  closure  weld  on  the  top 
deck  should  be  checked,  and  any  indication  of  general  structural  failure  noted. 

b.  More  data  must  be  found  on  the  lifetime  of  furfural  in  the  Arctic  Ocean  environment  and 
a  research  project  should  be  initiated  using  the  original  formulation  of  the  resin. 

c.  Similarly,  the  effective  lifetime  or  barrier  efficiency  presented  by  the  concrete  in  the  box 
needs  to  be  evaluated. 

44.  Further  work: 

a.  The  concepts  of  this  spreadsheet  model  developed  for  the  Lenin  could  equally  well  be 
applied  to  the  other  nuclear  material  that  still  remains  on  the  sea  floor.  It  is  the  intention  of  the 
Group  to  extend  the  study,  using  information  on  nuclide  inventories,  activation,  composition, 
structures,  and  the  containment  boundaries  of  all  the  disposed  material,  when  the  data 
becomes  available.  This  way,  a  realistic  estimate  can  be  made  for  the  radionuclide  release 
rates  into  the  Barents  and  Kara  Seas. 
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Table  1. 


Estimated  radioactivity  in  the  reinforced  concrete  and  stainless  steel  shell 

reactor  compartment  of  the  icebreaker  Lenin  for 


Radionuclide 

Fission  products 
'"'Cs 

137mB3 

“Sr 

soy 

“Kr 

’“"Pm 

®®Tc 

_ 

Subtotal  _ 

Activation  products 

“Ni 


Subtotal 

Actinides 

"“’Pu 

239pu 

^“°Pu 

^"Pu 

Subtotal 


Curies  (Ci) 

1.331x10“ 

1.259x10“ 

1.208x10“ 

1.208x10“ 

4. 768x1  CP 

2.969x10^ 

4.827x10’ 

2.371x10’ 

3.629x10° 

2.034x10° 

1.000x10’^ 


Activity 

Becquerels  (Bq) 

4.925x10’“ 

4.658x10’“ 

4.470x10’“ 

4.470x10’“ 

1.764x10’° 

1.099x10’° 

1.786x10’° 

8.773x10” 

1.343x10” 

7.526x10’° 

3.700x10* 


Percent  (%) 

26.14 

24.73 

23.73 
23.73 

0.937 

0.583 

0.095 

0.047 

0.007 

0.004 

0.00002 


5.091x10“ 

1.884x10’® 

100.00 

1.350x10° 

4.995x10’° 

77.23 

3.510x10° 

1.299x10’° 

20.08 

3.780x10’ 

1.399x10’° 

2.16 

9.180x10° 

3.397x10” 

0.53 

1.748x10° 

6.468x10’° 

100.00 

1.908x10° 

7.060x10’° 

82.13 

1.891x10° 

7.000x10’° 

8.14 

1.363x10° 

5.040x10’° 

5.87 

6.105x10’ 

2.260x10’° 

2.63 

2.856x10’ 

1.060x10’° 

1.23 

2.325x10° 

8.601x10’° 

100.00 

Table  2.  Estimated  radioactivity  in  the  three  reactor  pressure  vessels  within  the  reactor 
compartment  of  the  icebreaker  Lenin  for  1993/  ^^ 


Activity 

Radionuclide 

Curies  (Ci) 

Becquerels  (Bq) 

Percent  (%) 

Activation  products 

“Ni 

3.650x10^ 

1.351x10’“ 

77.23 

“Co 

9.490x10^ 

3.511x10’^ 

20.08 

“Ni 

1.022x10^ 

3.781x10’^ 

2.16 

2.482x10’ 

9.183x10” 

0.53 

Subtotal 

4.726x10^ 

1.749x10’“ 

100.00 

Figure  1:  The  LENIN  dump  site  in  Tsivolki  Inlet,  Novaya  Zemlya 


Mild  Steel  -  10mm 


Figure  2:  Diagramatic  representation  of  the  LENIN  reactor  compartment  (Based  on 
FR0f@r@nc@ 

SW  =  seawater,  SS  =  stainless  steel,  Furf  =  furfural,  TFCs=  Technological 
Fuel  Channels  (containing  the  damaged  fuel) 
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Figure  3:  The  four  scenarios  in  the  LENIN  Release  Model  (See 


para  23  for  description) 
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Figure  9:  Case  IV  Box  Containment  Barrier 
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Figure  10:  Case  V  Box  Support  Collapse 
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Rgure  1 1 :  Cases  I  to  V  All  Scenarios  Compared 
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Overview  of  Modeling  from  Land  Based  Sources  Through  the  Ob  and  Yenisei  Rivers 

and  Estuaries  to  the  Kara  Sea 

Terri  Paluszkiewicz 
Pacific  Northwest  Laboratories 
Marine  Sciences  Laboratory 
1529  W.  Sequim  Bay  Road 
Sequim,  Washington  98382 

The  objective  of  this  research  is  to  couple  models  with  existing  and  forthcoming  data 
to  estimate  radionuclide  fluxes  to  the  Kara  Sea  from  the  Ob  and  Yenisei  River  systems.  Our 
approach  is  to  use  a  hierarchy  of  simple  models  to  quantify  the  transport  and  fate  of 
radionuclide  contaminants  via  various  pathways  through  the  system.  The  system  consists  of 
the  sources  of  past  and  (possible)  future  contaminant  releases  from  Russian  nuclear  defense 
sites  on  land,  the  environment  (terrestrial,  ground  and  surface  water),  between  the  location  of 
the  sources  and  the  Ob  and  Yenisei  Rivers,  the  Ob  and  Yenisei  Rivers  and  their  tributaries, 
the  Ob  and  Yenisei  Estuaries,  and  the  confluence  region  where  the  estuaries  flow  into  the 
Kara  Sea.  There  are  five  integrated  tasks  that  together  provide  data  and  numerical  modeling 
to  quantify  the  fluxes  through  the  entire  system  and  each  component  (riverine,  estuarine,  or 
inner  shelf)  within  the  system. 

The  contaminants  that  we  are  following  through  the  system  are  ®°Sr,  ^^’Cs,  and  ^^^Pu. 
The  basic  algorithms  governing  the  adsorption  and  desorption  of  a  dissolved  contaminant  with 
sediments  and  solids  are  used  in  the  budget  models.  We  are  evaluating  the  following 
scenarios  for  both  river  systems:  a  steady  release  of  contaminants  and  a  catastrophic  release, 
corresponding  to  a  flood  and  dam  failure  and  to  the  fast  (in  geological  time)  movement  of 
groundwater  from  injection  wells.  Both  scenarios  will  be  evaluated  over  all  seasons.  The 
steady  release  simulations  will  yield  a  baseline  for  comparison  with  samples  being  taken  in 
the  estuary  and  on  the  shelf  by  other  program  participants. 

An  initial  quantification  has  been  made  of  the  inventory  at  nuclear  facilities  that 
include  research  institutes,  fuel  fabrication  facilities,  uranium  mining  and  materials, 
production  reactors,  and  spent  nuclear  fuel  reprocessing  sites.  The  latter  sites,  Mayak, 
Tomsk-7,  and  Krasnoyarsk-26  are  believed  to  be  responsible  for  the  majority  of  the 
radioactive  contamination  that  is  found  in  the  major  river  systems  that  feeds  into  the  Arctic 
Ocean  through  the  Kara  Sea.  An  analysis  of  basic  residence  times  under  simplified 
conditions  has  been  used  to  establish  that  the  models  were  configured  properly  for  each  river 
system.  Multiple-parameter  range  studies  were  run  to  determine  the  sensitivity  to  ranges  in 
the  empirical  coefficients  used  in  the  sediment  transport,  sediment-contaminant  transport 
algorithms.  Initial  results  indicate  the  need  for  further  data  on  the  composition  of  bed-load 
sediments  and  the  need  to  use  a  calibration  data  set  to  establish  the  mass  transfer  coefficient. 


Assessment  of  the  Contamination  of  the 
Kara  Sea  from  Radionuclide  Transport 
through  the  Ob  and  Yenisei  River  Systems 

and  amended  scope  of  work 

Refinement  and  Data  Calibration  of  the 
Radionuclide  Transport  Model  for  the 
Ob  and  Yenisei  River  Systems 


Participants: 

Terri  Paluszkiewicz 
Peter  Becker 
Don  Bradley 
Don  Denbo 
Mike  Foley 
Lyle  Hibler 
Chris  Sherwood 
Yasuo  Onishi 
Gene  Whelan 

Collaborators 

Ted  Callender 
Ed  Landa 
Robert  Meade 
TYPHOON 

Soviet  Hydrological  Institute 


The  study  will  address  the  following  specific  questions: 


•  Under  conditions  of  normal  river  flow,  with  existing 
dam  structures,  and  a  steady  continuous  release  of 
contaminants,  what  are  the  concentration  levels  of 
^Sr,  ^^^Cs  and  ^^^Pu  in  the  Ob  and  Yenisei  Rivers, 
their  estuaries,  and  the  inner  Kara  Sea? 


•  Under  catastrophic  conditions,  such  as  a  dam-break, 
tank-failure,  or  flood,  (i.e.  a  large  pulse-like  release  of 
contaminants),  what  are  the  concentration  levels  of 
^Sr,  ^^^Cs  and  ^^^Pu  in  the  Ob  and  Yenisei  Rivers, 
their  estuaries,  and  the  inner  Kara  Sea? 


Project  Objectives 


Couple  simple  models  with  existing  and  forthcoming  data 
to  estimate  radionuclide  fluxes  to  the  Kara  Sea  from  the  Ob 
and  Yenisei  River  systems. 

Answer  the  following  questions  using  existing  data, 
observations  that  are  presently  being  collected,  and  simple 
but  robust  budget  models. 

-  What  are  the  basic  pathways  and  partitioning  of  the 
contaminants  (dissolved,  or  sorbed)  in  the  rivers,  estuary, 
and  inner  shelf? 

-  What  are  the  dominant  processes  that  control  these 
pathways? 

-  Which  processes  require  additional  study  and  what 
additional  data  is  most  critical? 


Task  Integration 


Planned  Approach 

Quantitative  analysis  of  existing  contaminant  sources. 

Collaboration  with  Lawrence  Berkeley  Lab.  (studies  of 
Mayak). 

Collaboration  with  FSU  Typhoon  scientists  on  their 
ongoing  riverine  and  estuarine  studies. 

Collaboration  with  USGS  Mobilization  and  Transport 
Study. 

Assessment  of  transport  risk  from  terrestrial  sources  to  the 
river  using  the  MEPAS  methodology. 

Refinement  of  contaminant  transport  parameterization  if 
data  are  sufficient. 


SPOT  Image  of  the  Mayak  Site 
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Riverine  Modeling 


Some  key  processes  in  the  conceptual  model: 

Time-dependent  transport  of  water,  dissolved 
contaminants,  sediments,  contaminants  sorbed  onto 
sediments 

Deposition  and  Resuspension  of  contaminants  and 
sediments,  contaminants  sorbed  onto  sediments 

Ice  formation  and  breakup,  ice  damming  of  flow,  ice- 
incorporation  of  sediment  and  contaminant 


Incremental  improvements  in  model  complexity: 

•  Phase  1: 

Evaluate  sediment  and  contaminant  transport 
parameterizations  under  simplified  (steady)  flow 
conditions; 

'Tune"  model  using  existing  data  and  understanding; 

Evaluate  data  and  parameterization  weaknesses; 

Determine  residence  time  and  basic  behavior  of 
system. 

•  Phase  2: 

Evaluate  the  system  response  for  time-dependent  flow 
conditions; 

Examine  the  system  response  for  key  transport 
scenarios  (flooding,  ice  damming,  freshet,  freeze-up); 

Incorporate  additional  data,  calibrate  and  validate 
when  possible. 

•  Phase  3: 

Evaluate  system  response  to  risk-related  scenarios 
(container-rupture,  dam-break,  flood)  using  full  set  of 
data  and  algorithm  improvements. 


3-D  Primitive  Equation  Model:  Mellor  and  Blumberg, 
SPEM,  Davies,  Semtner  and  Chervin 


Attributes: 


Detractions: 


-Reveals  more  complex  circulation  and 
interactions 

-Could  provide  circulation  estimates  to  box 
models  of  sediments  and  contaminants 
-Easily  accepted  by  the  unsuspecting 
viewer 

-Our  level  of  observation-based  knowledge 
would  not  support  a  consistency  check 
with  model  results. 

-High  level  of  complexity  in  adding  ice 
-Requires  considerable  data 
-Parameterizations  are  inappropriate 


Hierarchy  of  Estuarine  Models 

Box  Model:  Officer  (1980a, b)  Officer  and  Kester  (1991) 
Attributes:  -Data  requirement  matches  availability 

-Mcludes  sediment  transport 
-Inclusion  of  ice  possible  (Hamblin,  1989) 
-Will  reveal  basic  patterns  of  distribution 

Detractions:  -Assumes  steady  state 

-Assumes  tidal  mixing  <  horiz.  net  circ. 
ie,  gravitational  processes  are  dominant 
-Will  not  reveal  new  information  about  the 
estuarine  processes 

2-D  Numerical  Models:  Bowden  and  Hamilton  (1975),  Oey 
et  al.  (1985),  Festa  and  Hansen  (1976), 
RMA2-RMA4(ACOE) 

Attributes:  -Reveals  more  complex  interactions  in 

estuarine  circulation 
-Will  reveal  more  information  about 
circulation  processes 
-Computational  efficiency 
-More  responsive  to  complex  physical 
forcing 

Detractions:  -Requires  more  data  (N^,  Rj,  U,V) 

-More  complexity  in  including  sediment 
dynamics  (except  RMA4) 

-More  complexity  in  including  ice 
*Box  model  solution  of  Festa  and  Hansen  (1978)  has  been 
compared  with  box  model  of  Officer(1980) 


Estuarine  Modeling 

Identify  the  characteristics  of  the  estuary: 

Salt-wedge  (river-flow  dominant) 

Two-layer,  entraining  (river-flow  dominant  modified 
by  tidal  currents 

Two-layer  with  vertical  mixing  (river  flow  and  tidal 
mixing) 

Vertically  homogeneous 
with  lateral  variation 
laterally  homogeneous 

Identify  key  processes  in  conceptual  model: 

Time-dependent  transport  of  fresh  and  salt  water, 
dissolved  contaminants,  sediments,  contaminants 
sorbed  onto  sediments 

Advection,  dispersion,  and  diffusion  due  to  tidal 
forcing,  wind-forcing,  river  inflow 

Ice  formation  and  breakup,  ice  damming  of  flow,  ice- 
incorporation  of  sediment  and  contaminant 

Scenarios  to  evaluate: 

Wind  mixing  and  denevilation 

Brine  rejection  effects  on  stratification  during  ice 
season 

Degree  of  homogeneity  during  freshet 


SEDIMENT  TRANSPORT  PROCESSES 


Description  of  the  Yenisei  Estiiary 

Primary  dynamic  influences:  river  runoff^  wind^  tides,  ice 

360  km  long,  35-100  km  wide 

Inflow:  20,000  m^  s'^  average 
8,000  m^  s"^  low 
80,000  m^  s'^  high 

Primarily  semidiurnal  tides 

Flow:  is  dominated  by  river  discharge  with  local  currents 
up  to  50  cm  s’^  northward; 

Wind-driven  ”pile-ups”  influence  the  head  of  the  estuary 
with  sea-level  elevations  of  up  to  6.0  m  during  minimum 
discharge  and  strong  fall  winds  15-20  m  s’^;  related 
currents  up  to  8  cm  s’^;  long-wave  propagation  noted. 

Mixing  Regime:  the  Yenisei  estuary  is  stratified  near  the 
mouth.  The  extent  of  the  salt  wedge  is  not  well 
documented 

Longitudinal/Latitudinal  Variations:  No  indications  of 
latitudinal  variations;  indications  are  that  in  the  shallower 
portions  of  the  estuary  the  waters  are  more  well-mixed  and 
homogenous. 


Bottom  Salinity,  August  1993 


Data  Courtesy  of  J.  Brooks,  GERG 


Bottom  Salinity,  August  1993 


70*  7Z  74*  76*  78*  86° 


Data  Courtesy  of  J.  Brooks,.  GERG 


Description  of  the  Ob  Estuary 

Primaiy  dynamic  influences:  river  runoff,  wind,  tides,  ice 
800  km  long,  60-100  km  wide 

Inflow.  5j0.5  km^  yr  ^  (16,800  m^  s'^)  average 
404.0  km^  yr"^  (12,800  m^  s‘^)  low 
662.0  km^yr'X21,000  m^  s’^)  high 

Volume  of  Ob  is  450  km^;  waters  renewed  completely 
during  the  year 

Primarily  semidiurnal  tides;  max  amplitude  1.85  m  at 
mouth  decreasing  to  0.30  m  at  head 

Flow:  is  dominated  by  river  discharge  with  local  currents 
up  to  50  cm  s  northward  during  freshet  and  <  2  cm  s'^ 
during  minimum  discharge; 


Wind-driven  "pile-ups”  influence  the  head  of  the  estuaiy 
with  sea-level  elevations  of  up  to  4.0  m  during  minimum 
discharge  and  strong  fall  winds 

Mixing  Regime:  the  Ob  estuaiy  is  well  mixed 
(homogenous  from  the  head  to  71  N)  and  somewhat 
stratified  near  the  mouth.  The  extent  of  the  salt  wedge  is 
not  well  documented 

L^ngitudinal/Latitudinal  Variations:  Indications  are  that 
the  freshwater  flow  is  right-bounded  along  the  estuary  and 
tends  right  upon  reaching  the  confluence  region. 


o 


Hierarchy  of  Estuarine  Models 

Box  Model:  Officer  (1980a, b)  Officer  and  Kester  (1991) 
Attributes:  -Data  requirement  matches  availability 

-Includes  sediment  transport 
-Inclusion  of  ice  possible  (Hamblin,  1989) 
-Will  reveal  basic  patterns  of  distribution 

Detractions:  -Assumes  steady  state 

-Assumes  tidal  mixing  <  horiz.  net  circ. 
ie.  gravitational  processes  are  ddminant ' 
-Will  not  reveal  new  information  about  the 
estuarine  processes 

2-D  Numerical  Models:  Bowden  and  Hamilton  (1975),  Oey 
et  al.  (1985),  Festa  and  Hansen  (1976), 
RMA2-RMA4(ACOE) 

Attributes:  -Reveals  more  complex  interactions  in 

estuarine  circulation 
-Will  reveal  more  information  about 
circulation  processes 
-Computational  efficiency 
-More  responsive  to  complex  physical 
forcing 

Detractions:  -Requires  more  data  (N^,  Rj,  U,V) 

-More  complexity  in  including  sediment 
d3mamics  (except  RMA4) 

-More  complexity  in  including  ice 
*Box  model  solution  of  Festa  and  Hansen  (1978)  has  been 
compared  with  box  model  of  Officer(1980) 
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Preliminary  Assessment  of  Radioactive  Contaminant  Transport  from  Catastrophic 
Release  Scenarios  in  the  Ob  and  Yenisei  River  Systems 

L.  F.  Hibler 
Battelie/MSL 

1529  W.  Sequim  Bay  Road 
Sequim,  Wa  hington  98382 
(206)081-3642 

e-mail:  lf_hibler@ccmaiLpnl.gov 

M.  C.  Richmond 
Battelle 

Richlar  d,  WA  99352 
(20u)376-8319 

There  are  significant  inventories  of  radioactive  waste  stored  in  the  Ob  and  Yenisei 
watersheds  of  the  Siberian  Arctic.  The  purpose  of  this  study  was  to  provide  a  defensible 
means  of  estimating  the  contamination  that  would  reach  the  Kara  Sea  under  a  catastrophic 
release  scenario  from  three  sites  (Krasnoyarsk,  Tomsk,  and  Mayak).  A  dynamic,  one¬ 
dimensional  sediment/contaminant  transport  model  was  used  to  study  the  transport  of  Cs-137 
under  a  variety  of  flow  conditions.  Analysis  of  the  sediment/contaminant  transport  algorithms 
showed  that  the  interaction  of  key  parameters  (critical  shear  stress,  erodibility,  and 
partitioning  and  rate  constants)  had  a  pronounced  effect  on  the  calculated  contaminant 
transport.  A  suite  of  sensitivity  studies  was  performed  to  determine  acceptable  ranges  for 
these  parameters,  and  these  ranges  were  compared  with  those  from  similar  analyses  of  other 
river  systems.  This  transport  analysis  provided  contaminant  boundary  conditions  for  the 
complementary  studies  of  the  estuarine  and  shelf  contaminant  transport. 

This  research  was  supported  by  the  Office  of  Naval  Research  under  Contract  25110. 


Oh  and  Yenisei  River  Modeling  Study 


Study  Goals: 

Estimate  the  Fate  and  Behavior  of  Radioactive  Contamination 
Released  from  Krasnoyarsk( Yenisei),  Mayak  and  Tomsk  (Ob) 

Provide  Estimated  Boundary  Conditions  for  an  Estuarine  Model 


Approach: 

Use  One-Dimensional  Dynamic  Models  as  Tools  for 
Objective  Characterization  of  Sediment  and  Contaminant 
Transport  within  the  Ob  and  Yenisei  Rivers 

Status: 

Performing  Sensitivity  Studies  on  the  Yenisei  River  Application 

Configuring  similar  study  for  the  Techa-Iset-Tobol-Irtysh-Tom- 
Ob  Application 


Future  Work: 

Use  of  CHARIMA  Model  for  linked  Hydrodynamic  and 
Transport  Modeling 


Ob  River  Discharge  at  Salekard 


Yenisei  River  Discharge  at  Igarka 


Annual  Suspended  Sediment  Load  (10®  tons) 


Annual  Runoff  (km®) 
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Contaminant  Inventory  Distribution 


Percent  of  Inventory 


Contaminant  Inventory  Distribution 


Modeled  Contaminant  and  Sediment  Process: 


Sediment  Transport: 

[Time  rate  of  change  in  Suspended  Sediment  Concentration]  + 

[Advection]  =[Dispersion]  +[Sources]  +[Erosion]  +[Deposition] 


Dissolved  Contaminat  Transport: 

[Time  Rate  of  change  in  Dissolved  Contaminant  Concentration] 
[Advection]  =[Dispersion]  -Pecay]  +  [Sources]  - 
[Sorbed  Phase  Exchange  From  Suspended  Sediment]  - 
[Sorbed  Phase  Exchange  From  Bed  Sediment] 


+  UA^  =  ^  eA^  -  XAC 
dt  dx  dxV  dxy 

-  Qf  -  AK  {KdSC  -  G)  -  5T  ( 1  -  Ti)  dKi,  (K^C  -  G^) 
Sediment-Sorbed  Contaminant  Transport: 

[Time  Rate  of  Change  in  Sediment  Sorbed  Contaminant  Concentration]  = 
[Advection]  = 

[Dispersion]  -[Decay]  +[Sources]  -{-[Erosion]  -Peposition]-{- 
[Dissolved  Phase  Exchange] 


dG  TT  d( 

+  UAr^  =  ^  £A-^  -  aAG 
dt  dx  dxy  dxJ 


QIG  +  B  GbSr 


SdG 
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Summer  (High  Flow)  Hydrology: 


(a) 

(b) 

(c) 

Length  (km) 

270 

460 

580 

620 

Slope 

0.00016 

0.00006 

0.00009 

0.00002 

Area  (m^) 

64977 

93733 

106822 

130851 

Depth  (m) 

19.3 

22.7 

24.0 

26.3 

Vol  (km^) 

18 

43 

62 

81 

Q  (m^/s) 

39353 

57050 

65131 

80000 

0.61 

0.61 

0.61 

0.61 

T  (days) 

5 

9 

lO 

12 

(dynes/cm  ) 

303 

133 

212 

53 
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SPOT  Image  of  the  Mayak  Site 
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Transport  Model 
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The  NRL  NORAPS  Model 


William  Thompson 
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ABSTRACT 


A  version  of  the  Navy  Operational  Regional  Atmospheric 
Prediction  System  (NORAPS)  will  be  used  to  provide  atmospheric 
forcing  to  PIPS  over  the  Kara  Sea.  NORAPS  is  a  hydrostatic 
primitive  equation  model  with  split  explicit  time  differencing 
and  time  dependent  lateral  boundary  conditions  from  NOGAPS.  In 
the  analysis  and  initialization  phase,  an  optimum  interpolation 
scheme  is  used  to  blend  current  observational  data  with  a 
previous  12  h  forecast  and  an  incremental  update  is  used.  The 
physical  parameterizations  used  in  this  version  of  the  model 
include  a  second-order  closure  turbulence  scheme,  fully 
interactive  clouds  and  radiation,  and  a  surface  energy  budget  for 
snow  and  ice  covered  surfaces.  The  physical  parameterizations 
use  implicit  time  differencing. 

For  the  Kara  Sea  runs,  20  km  horizontal  resolution  will  be 
used  with  a  domain  extending  2100  km  E-W  by  1600  km  N-S.  The 
domain  thus  extends  from  the  Eastern  Barents  Sea  to  the  western 
Laptev  Sea  and  from  approximately  75N  to  83N. 


William  T.  Thompson 
Naval  Research  Laboratory 
Marine  Meteorology  Division 
7  Grace  Hopper  Ave. 
Monterey,  CA  94943-5502 
(408)  656-4733 
(408)  565-4769  (fax) 
thompson@nr Imry . navy .mil 


COASTAL  MESOSCALE  DATA  ASSIMILATION  SYSTEM 

ANALYSIS /INITIALIZATION 

•  Optimum  Interpolation  Analysis 

•  Incremental  Update 

•  Nonlinear  Vertical  Mode  Initialization 
FORECAST  MODEL  FEATURES 

•  Hio-h  Vertical  and  Horizontal  Resolution 

20  km  horizontal  resolution 

30  levels  in  the  vertical;  14  in  the  lowest  1  km 

•  Hydrostatic 

•  Primitive  Equation 

.  Time  Dependent  Lateral  Boundary  Conditions  from  the  Global  Model  (NO¬ 
GAPS) 


•  Split  Explicit  Time  Differencing 


MODEL  PHYSICS 

.  Second-Order  Closure  Turbulence  Parameterization 
Mellor-Yamada  (1974)  level  3 
•  Fully  Interactive  Cloud- Radiation 
Two-Stream  Radiative  Transfer 

Cloud  Parameterization  follows  Sommeria  and  Deardorff  (19 1 9) 

•  Surface  Layer  Parameterization 

Loius  (1979)  over  land  and  sea  ice 
Liu  (1979)  over  open  water 

•  Implicit  Time  Differencing  for  Physics 


Izk  U?  hjt 


NORAPS  grid  1,  109  x  82  x  36  20.00  km 

500.0  mb  geop  ht 
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ABSTRACT 


Use  of  Multi-channel  NOAA  AVHRR  Data  in  Studying  River  Outflow 
Effects  and  Ice  Thickness  Patterns 

Robert  W.  Fett 
SAIC/Consultant 

Characteristic  albedo  and  temperature  changes  have  been  used 
to  evaluate  ice  thickness  in  NOAA  AVHRR  data  in  the  Banks  Island  - 
Beaufort  Sea  region.  Similar  attempts  have  been  made  in  the  past 
using  research  aircraft.  A  "spin-off"  from  this  research  using  NOAA 
data  in  the  Mackenzie  Bay  region  revealed  an  anomalous  warm  plume 
extending  several  hundred  kilometers  westward  along  the  North  Slope 
of  Alaska.  This  plume  is  believed  to  be  an  outflow  effect  of  the 
Mackenzie  River.  Reduced  albedo  in  the  plume  may  be  a  result  of 
river  sediment  which  has  become  embedded  in  the  ice  and  carried 
westward  by  the  Beaufort  anticyclonic  gyre.  The  example  suggests 
that  river  outflow  effects  and  the  effect  of  anomalous  warm  spots 
in  the  Laptev  Sea  may  also  be  monitored  using  these  data. 


Table  1 .  Sea  ice  classification  scheme  (World  Meteorological  Organization  1970). 
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The  Sea  of  Okhotsk,  located  in  the  northwestern  corner  of  the  north  Pacific  Oceaii 
between  45°N  and  65°N,  is  covered  by  low-salinity  water  in  the  surface.  This  large  margina 
sea  possesses  a  quite  unique  oceanographic  setting  that  is  not  found  elsewhere  on  a  similar 
scale.  First,  it  is  mostly  covered  by  sea  ice  during  winter  despite  its  relatively  southern 
location.  Second,  it  exhibits  a  very  strong  seasonal  cycle.  The  sea  surface  temperature 
fSST)  varies  from  over  15°C  at  the  summer  to  near  freezing  point  at  winter.  Unlike  tne 
Arctic  Ocean,  Okhotsk  sea  ice  completely  melts  in  the  summer  despite  its  extensive  coverage 
in  winter.  Third,  its  deep  basin  is  covered  by  a  dicothermal  layer  between  30  and  150  meters 
in  depth.  The  dicothermal  layer  water  usually  remains  near  freezing  through  the  entire  year. 
The  Sea  of  Okhotsk  also  shows  a  very  strong  interannual  variability  as  revealed  from  satellite 
passive-microwave  observation  of  sea  ice  and  from  a  few  available  in  situ  measurements. 

A  simple  coupled  sea  ice  and  ocean  model  is  used  to  examine  some  key  physical  pro¬ 
cesses  that  maintain  this  unique  oceanographic  state  and  to  explore  some  plausible  mec  a- 
nisms  that  contribute  to  the  observed  interannual  variability.  The  o^an  model  is  a  modified 
version  of  the  1.5-dimensional  model  of  Killworth  and  Snaith  (1984).  This  model  only  resolves 
laterally-integrated  vertical  profiles  of  temperature,  salinity  and  upwelling.  It  is  forced  by 
a  prescribed  inflow  of  the  Pacific  Ocean  water,  precipitation-evaporation,  runoff,  freshwater 
flux  associated  with  sea-ice  melting  and  freezing,  sinking  of  shelf  water,  and  inflow  from  the 
Sea  of  Japan.  An  outflow  to  the  Pacific  Ocean  is  also  specified.  A  thermodynamic  sea-ice 
model  (Welander,  1973)  is  used  to  calculate  the  seasonal  vmation  of  sea-ice  thickness  The 
forcings  are  specified  by  fitting  them  to  some  scarce  observational  data  mainly  from  published 
Japanese  and  russian  literatures.  Despite  its  simplicity,  this  model  successfully  captures  the 
main  features  of  the  Sea  of  Okhotsk.  For  example,  the  model  produces  a  dicothermal  layer 
between  50  and  150  meters,  a  transitional  layer  beneath  the  dicothermal  layer,  and  a  highly 
seasonal  surface  mixed  layer.  In  our  model,  the  dicothermal  layer  is  generated  by  the  winter 
convection,  and  the  transitional  layer  is  due  to  a  mixing  between  shelf  watCT  and  Pacific 
Ocean  intermediate  water.  These  model  results  agree  with  the  hypothesis  of  Kitani  (1973). 
More  interestingly,  the  model  also  exhibits  a  strong  interannual  variability  characterized  by 
alternating  between  two  deep  convections  to  the  base  of  the  dicothermal  layer  and  two  sh^- 
low  convections  to  the  base  of  the  mixed  layer.  After  two  successive  deep  convections,  the 
dicothermal  layer  is  cold  and  salty.  Sea  ice  is  also  thicker  during  these  two  deep  convections. 
In  the  following  two  years,  the  increased  melting  water  to  the  imxed  layer  together  with  a 
denser  dicothermal  layer  restrict  a  winter  convection  only  within  the  surface  layer.  Deep 
convection  resumes  after  two  years  of  shallower  convection. 
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38.  A  scheme  of  general  circulation  of 
water  masses  in  the  surface  of 
the  Sea  of  Okhotsk. 
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The  Sea  of  Okhotsk  and  the  Pacific  are  connected 
by  14  channels  along  Kurile  Islands.  Inflow  mainly 
occurs  at  Bussol  Strait  and  Kruzenshtern  Strait; 
outflow  at  Etorofu  Strait,  Bussol  Strait,  and 
shallow  western  Kurile  channels. 


Ocean  model  (a  modified  Killworth  and  Smith  model): 


(A(z)W(z)),  =  Qiz)  (1) 

ATt  +  (AWT),  =  k(AT,),  +  Qt  +  mix  (2) 

ASt  +  (AWS)z  =  Ki^Sz)z  +  Qs  +  mix  (3) 

where: 

Q{z,  t)  =  Qpac  —  Qout  +  Qsoya  +  Qshf  +  Q surface 

Qt(z,  t)  =  QpacTpac  —  QoutT  +  QsoyaTsoya  +  Qshf'I'shf  +  Q^surface 

(^si^Z ^  —  QpacTpac  Qout^  "f”  Q soya^soya  ”1”  Q shf^shf 

Qmrface  =  ■A(0)(-f’(*)  ~  “  ^<)  +  Qnmcff(t) 

Qt..rf.c,  =  4(0)(P(<)Tp  -  E(t)%  -  6,Tf)  +  Q  runo  ffij'^Trxinoff 
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Future  work: 


•  Modelling  wind-stress-driven  and  buoyancy-driven  circulations  in  the  Sea 
of  Okhotsk; 

•  coupling  a  three-dimensional  OGCM  to  a  dynamic  and  thermodynamic 
sea-ice  model; 

•  to  understand  the  physical  processes  that  maintain  the  unique  oceano¬ 
graphic  state  in  this  large  marginal  sea; 

•  to  examine  various  possible  mechanisms  that  contribute  to  a  strong  in¬ 
terannual  variability. 
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ABSTRACT 


Three-dimensional  features  of  dense  bottom  plumes  and  deep  water  formation  over 
continental  slopes  with  or  without  along-slope  topographic  variations  are  investigated  by 
simulating  the  evolution  of  a  density  front  at  a  southern  ocean  continental  shelf  break,  using  a 
three-dimensional,  primitive  equation  numerical  model  with  a  second-order  turbulence 
closure  scheme  embedded.  The  focus  of  our  investigation  is  the  role  of  topography  in 
determining  mixing  and  deep  water  formation  during  the  transient  adjustment  process  of  a 
density  front  over  a  continental  slope.  We  compare  and  discuses  the  numencal  simulations 
for  two  cases:  a  uniform  shelf  and  slope  case  and  a  case  with  a  canyon  that  leads  from  the 
coast  to  the  deep  ocean  crossing  the  shelf  and  the  slope. 

The  numerical  simulations  show  that  baroclinic  instability  together  with  the  rotation 
deform  the  front  into  surface  and  bottom  plumes  along  the  shelf  break.  Correspondingly,  a 
street  of  counter-rotating  eddies  form  both  at  the  surface  and  the  bottom  near  the  shelf  break. 
The  buoyant  surface  plumes  are  deflected  to  the  left  while  intruding  shoreward  in  this 
southern  hemispheric  case,  with  an  intrusion  scale  controlled  by  rotation  and  buoyancy.  The 
dense  bottom  plumes,  which  are  dominated  by  the  bottom  friction  and  buoyancy,  however, 
can  penetrate  downslope  to  much  larger  extent  than  do  the  surface  plumes.  Ultimately 
diffusion  will  Umit  the  downward  penetration  of  the  shelf  water.  Eddy  transport  of  shelf 
water  is  found  to  be  very  important  in  deep  water  formation  on  the  continental  slope. 

The  presence  of  a  canyon  in  the  shelf  and  slope  enhances  the  penetration  depth  of  the 
dense  bottom  currents  in  the  canyon  and  thus  increases  deep  water  production,  due  to  partial 
cancellation  of  the  rotation  effect  in  the  lower  part  of  the  water  in  the  canyon.  A  deep  canyon 
has  more  significant  effects  in  enhancing  the  deep  water  formation  by  draining  more  shelf 
water  into  deep  ocean  through  the  canyon,  while  a  narrower  canyon  has  reduced  drainage. 
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Three-dimensional  features  of  dense  bottom  plumes  and  deep  water  formation  over 
continental  slopes  with  or  without  along-slope  topographic  variations  are  investigated  by 
simulating  the  evolution  of  a  density  front  at  a  southern  ocean  continental  shelf  break, 
using  a  three-dimensional,  primitive  equation  numerical  model  with  a  second-order 
turbulence  closure  scheme  embedded.  The  focus  of  our  investigation  is  the  role  of 
topography  in  determining  mixing  and  deep  water  formation  during  the  transient 
adjustment  process  of  a  density  front  over  a  continental  slope.  We  compare  and  discuses 
the  numerical  simulations  for  two  cases:  a  uniform  shelf  and  slope  case  and  a  case  with  a 
canyon  that  leads  from  the  coast  to  the  deep  ocean  crossing  the  shelf  and  the  slope. 

The  numerical  simulations  show  that  baroclinic  instability  together  with  the  rotation 
deform  the  front  into  surface  and  bottom  plumes  along  the  shelf  break.  Correspondingly, 
a  street  of  counter-rotating  eddies  form  both  at  the  surface  and  the  bottom  near  the  shelf 
break.  The  buoyant  surface  plumes  are  deflected  to  the  left  while  intruding  shoreward  in 
this  southern  hemispheric  case,  with  an  intrusion  scale  controlled  by  rotation  and 
buoyancy.  The  dense  bottom  plumes,  which  are  dominated  by  the  bottom  friction  and 
buoyancy,  however,  can  penetrate  downslope  to  much  larger  extent  than  do  the  surface 
plumes.  Eddy  transport  of  shelf  water  is  found  to  be  very  important  in  deep  water 

formation  on  the  continental  slope.  .  .  u  r  u 

The  presence  of  a  canyon  in  the  shelf  and  slope  enhances  the  penetration  depth  of  the 
dense  bottom  currents  in  the  canyon  and  thus  increases  deep  water  production,  due  to 
partial  cancellation  of  the  rotation  effect  in  Ae  lower  part  of  the  water  in  the  canyon.  A 
deep  canyon  has  more  significant  effects  in  enhancing  the  deep  water  formation  by 
draining  more  shelf  water  into  deep  ocean  through  the  canyon,  while  a  narrower  canyon 
has  reduced  drainage. 


1.  Introduction 

The  continental  shelves  in  southern  oceans  comprise  a  vast  area  and  play  an 
important  role  in  water  mass  formation  and  biological  productivity.  Antarctic  Bottom 
Water  forms  in  the  southwestern  Weddell  Sea  as  a  result  of  frontal  mixing  at  the 
continental  shelf  break  between  the  Warm  Deep  Water  and  the  Western  Shelf  Water,  as 
well  as  the  result  of  deep  open  ocean  convection  [Foster  and  Carmack  1976;  Killworth 

1983;  Foldvik  et  al.  1985].  ...  v 

Here  we  focus  our  attention  on  the  deep  water  formation  by  dense  shelf  water 
flowing  down  a  continental  slope  and  mixing  with  the  deep  ocean  interior  water.  Cold 
and  dense  shelf  water  is  an  important  source  of  deep  water  even  in  the  absence  of  deep 
open  ocean  convection  [Foster  and  Carmack  1976;  Foldvik  et  al.  1985].  Deep  water 
formation  on  a  continental  slope  has  unique  dynamic  characteristics  that  differentiate  it 
from  open  ocean  deep  water  formation.  The  presence  of  a  continental  slope  introduces  a 
downslope  gravitational  component  that  enables  cross-isobath  transport  of  dense  shelf 
water.  Frictional  bottom  boundary  layers  also  contribute  to  cross-isobath  transport  of 
mass,  heat,  salt.  The  shelf  break  constrains  the  flow  around  it  and  allows  possibility  of 


trapped  currents.  Submarine  canyons  iii  the  continental  shelf  and  slope  may  facilitate  the 

offshore  transport  of  dense  shelf  water  into  deep  ocean.  .  i  ,  t  a 

Earlv  analytical  studies  on  dense  bottom  plumes  on  continental  slopes  wei-e  steady 
stream-tube  model  simulations  [Smith  1975;  Killworth  1977]  which  exclude  many 
nrocesses  that  we  find  to  be  important.  Smith  [1975]  descnbed  a  rotating  plume  model 
that  explains  well  the  behavior  of  the  Norwegian  and  Mediterranean  outflows.  Killwort 
[1977]  considered  both  two-dimensional  and  three-dimensional  plume  inodels  to 
investigate  mixing  processes  on  the  Weddell  Sea  continental  slope  and  found  that  the 
stream  tube  models  gave  results  incompatible  with  observations,  until  the  themaobancity 
Swas  incorporated  in  the  model.  The  stream  tube  model  has  several  limitations  in 
dealing  with  outflow  along  very  wide  shelves  with  complex  topography  (e.g  the 
Antarctic  continental  shelves).  Firstly  it  assumes  that  the  ambient  water  is  quiescent  and 
does  not  allow  baroclinic  instability  and  cross-stream  vanations,  and  thus  it  neglects 
some  important  dynamic  features  of  three-dimensional  plumes.  Secondly  it  can  not  deal 
with  plumes  flowing  on  a  slope  with  abrupt  bottom  variations,  such  as  a  canyon  or  a  sea 

Recently  more  complicated  numerical  models  have  been  used  to  investigate  dense 
water  formation  and  transport  on  continental  shelves  and  slopes.  Jungclaus  and 
Backhaus  [1994]  employed  a  hydrostatic,  reduced  gravity,  two-diraensmnal  priimtive 
equation  model  to  simulate  the  nonlinear  and  non-steady  features  of  the  Denm^k  Strmt 
rWerflow  Hsu  [1993]  used  a  two-dimensional  primitive  equation  model  to  study 
negative  buoyancy  forcing  adjacent  to  a  coast  and  found  that  the  dense  water  did  not 
move  offshore  beyond  the  forcing  region  unless  an  along-shelf  current  was  included. 
The  model  did  not  allow  along-shelf  variation  in  the  flow  field  and  thus  could  generate 
Suher  pWs  nor  eddies,  olwarkiewicz  and  Chapman  [1994]  modeled  dense  water 
formation  and  transport  generated  by  an  idealized  coastal  poly nya  over  a  gently  slopi  g 
shallow  continental  shelf.  The  flow  response  to  a  constant  negative  surface  buoyancy 
flux  imposed  over  a  half-elliptical  region  adjacent  to  the  coast  was  discussed.  The 
numerical  results  suggest  that  instability  and  eddy  fluxes  are  the  important  mechanisms 
in  the  transport  of  dense  coastal  water  off  continental  shelves  and  into  marginal  seas.  To 
investigate  offshore  transport  of  dense  bottom  water  and  deep  water  formation  on 
southern  ocean  continental  slopes,  we  pursue  a  similar  approach  but  focus  on  t 
transient  response  to  topographic  influences  on  bottom  plume  formation  and  penetration. 
^  In  this  study  we  simulate  the  evolution  of  a  density  front  at  a  southern  ocean 
continental  shelf  break  with  a  three-dimensional  pnmitive  equation  model  m  which  a 
Scond-order  turbulence  closure  scheme  is  embedded.  We  consider  a  wide  continental 
shelf  and  a  continental  slope  for  two  cases:  a  uniform  shelf  and  slope  c^e  and  a  case 
with  a  canyon  in  the  shelf  and  the  slope.  The  traditional  "dam  b^^  techmijue  is 
adopted  to  simulate  the  adjustment  processes  of  a  density  front  at  Ae  shelf  break  in  the 
absSice  of  external  forcing.  In  the  dam-break  expenment  a  vertical  bamer  ^at  sep^ates 
different  salinities  and  temperatures  of  water  at  the  shelf  break  is  removed  at  JJJ”® 
allow  the  separated  waters  to  exchange.  This  technique  gives  a  good  sense  of  fb^ 
development  of  the  response  of  the  front  over  the  continental  slope  from  an  initial  rest 

state. 

2.  Model  Description 

To  investigate  the  topographic  effects  on  the  dynamic^  processes  on  the  shelf  an 
slope  two  idealized  topographies  are  introduced  in  the  model  domain.  The  first  one  is  a 
uniform  shelf  and  slope  case  which  has  no  along-shelf  topographic  vanations,  as  show 
in  Figure  1.  The  second  case  has  a  100-km-wide  and  200-m-deep  canyon  onented  in  the 
East- West  direction,  halfway  between  the  southern  and  the  northern  boundanes,  as 


shown  in  Figure  2.  In  this  model  study,  the  water  mass  system  is  greatly  simplified,  and 
we  consider  only  the  case  for  the  adjustment  of  a  density  front,  located  at  the  shelf 
break  separating  a  well-mixed  shelf  water  and  a  well  mixed  deep  ocean  intenor.  We  do 
this  in  order  to  gain  physical  insight  into  the  dynamics  of  the  plumes  and  the 
topographic  influences  on  die  plume  formation  and  penetration  and  as  a  Prerequisite  for 
ui^erstanding  the  more  complicated  case  with  stratification  m  both  the  shelf  and 
slope  water.  The  temperature  and  salinity  ranges  used  here  are  those  from  the  wa^er 
masses  involved  in  forming  deep  water  in  the  Antarctic,  i.e.,  temperature  from  -2  to  +1 
C  S^inity  from  34.40  to  35.00  %o  [Foster  and  Carmack  1976].  Figure  3  shows  the 
idealized  croL-shelf  topography  of  the  model  domain  without  along-shelf  topographic 
variations  and  the  hypothetical  initial  temperature  and  salinity  conditions  from  which  the 

evolution  of  the  density  front  starts.  _  .  i,  „i, 

The  model  used  here  is  a  three-dimensional  primitive  equation  model  which  is 
hvdrostanc  fully  nonlinear,  and  applies  the  Boussinesq  approximation.  Vertical 
turbulent  mixing  processes  are  parameterized  with  the  second-order  turbulent  closure 
Sterne [1982],  Smagorinsky's  [1963]  lateral  eddy  mmng 
coefficient  is  used  along  sigma-coordinate  surfaces  with  a  purpose  of  purely  removing 
pSs  compu^^^^^^^^  The  model  uses  a  sigma-coordinate 

model  equations  and  the  method  of  solution  are  given  in  Blumberg  and  Mellor  [1983]. 

High  resolution  in  the  horizontal  and  the  vertical  are  used  to  ensure  that  instability 
processes  due  to  both  the  horizontal  and  vertical  velocity  shears  and  density  padients 
Le  properly  included.  Altogether,  41  sigma  levels  are  used  in  the  vertical,  with  a 
logarithmic^sigma  distribution  at  the  surface  and  the  bottom  in  ^er  to  resolve  the 
buoyant  near-surface  plumes  and  the  bottom  boundary  currents  The  honzontal  grid 
sizes  are  Ax=Ay=lO  km  throughout  the  model  domain  which  is  able  to 
0(i?o)=ll  km  for  the  case  considered.  The  value  of/o  is  set  a  latitude  of  72  S,  so  tha 

At  the  free  surface,  there  is  no  surface  wind  stress,  and  no  heat  and  salinity  fluxes.  At 
the  bottom  zero  heat  and  salinity  fluxes  are  specified.  At  land  bounces  we  use  the 
condition  of  no  diffusive  fluxes  of  any  property  across  the  interface,  pie  extern  open 
boundary  is  set  such  that  interior  energetic  bores  are  allowed  to  pass  through  with  httle 
reflection  to  the  interior.  The  model  is  integrated  with  an  external  time  step  of  10  s  and  a 
tatS  time  step  of  500  s,  which  is  set  hy  the  CFL  condition  of  the  fast-moving  external 

mode. 

3.  A  uniform  slope  case 

The  uniform-slope  case  is  taken  as  the  benchmark  experiment,  whose  results  wil^^e 
compared  with  the  those  for  the  cases  with  along-shelf  topographic  yanattons.  The 
salimtv  fields  at  sigma  level  1  (the  surface),  where  the  surface  outflow  is  maximum,  at 
da^M  18  and  I2,  are  shown  in  Figure  4.  After  the  sudden  removal  of  the  -dam 
fresher  and  warmer  deep  ocean  water  intrudes  onto  the  shelf  ne^  the  surface,  while 
more  saline  and  colder  water  penetrates  downslope  from  below.  The  initially  straight 
front  meanders  and  deforms  into  a  number  of  light  plumes  at  the  surface,  with  a  length 
scale  of  about  50  km  in  the  along-shelf  direction.  The  surface  salinity  fields  display 
strone  along-shelf  asymmetry  in  the  plumes.  The  plumes  are  almost  uniformly  aligned 
^ng^the  sLlf  break^ey  are  deflected  to  their  left-hand-side  as  they  flow  onto  the 
shelf  aoparentlv  dominated  by  the  planetary  rotation.  The  surface  plumes  continue  to 
spread  ^oreward,  but  the  spreading  speed  slows  and  the  intrusion  is  limited  as  the  front 
becomes  diffused.  Eventually  a  quasi-steady  state  is  reached  where  the  surface  Plumes 
are  trapped  near  the  shelf  break  with  the  extent  of  shoreward  intrusion  of  about  50-60 
km  in  the  interior  of  the  model  domain. 


The  salinity  fields  at  level  40  (the  bottom  level)  at  days  14,  18,  a*"®  illustrated  by 

Fipiire  5  resoectively.  Again  the  salinity  field  and  the  temperature  field  have  alinost 
tso  later  on  only  one  field  will  be  shown).  The  dense  cold,  and  saline 
DluSo™  0Mh°  S  and  penetrate  away  from  the  shelf  br^  down  the 
cinnp  Diirine  the  initial  stage  of  the  bottom  plume  formation,  a  street  of  dense  plumes 
anoears  in  the  form  of  individual  plumes  with  a  width  of  about  50  km,  uniformly  ^l^gued 
afone  the  shelf  break  while  flowing  downslope.  They  are  deflected  to  the  left  slightly 
under  the  influence  of  rotation.  Bottom  friction  increases  with  the  increase  of  the  current 
soeed  and  it  balances  the  Coriolis  force  in  dragging  the  plumes  further  down  the  slo^. 
The  bottom  plumes  continue  to  penetrate  downslope,  but  the  spreading  speed  slo 
eraduallv  with  time  due  to  mixing  and  entrainment  with  the  surrounding  water.  Later  at 
fcl8  davs  the  bottom  plumes  grow  larger  and  wider  as  a  result  of  mixing  and  they  sta 
^  with  each  other^  Gradually  the  bottom  plumes  lose  their  individual 

characS^^^^  ^ith  a  wavy  foreface  which  is  almost 

shelf  break.  The  horizontal  penetration  extent  of  the  bottom  plumes  is 
about  140  km  from  the  shelf  break,  corresponding  to  a  vertical  penetration  depth  of 

^'^^To^understand  the^dynamics  of  the  spreading  plumes,  we  examine  the  evolution  of 
the  lateral  velocity  fields  at  the  surface.  Figure  6  shows  the  velocity  fields  at  days  14  18 
22  at  sigma  level^  1.  Baroclinic  instability  causes  a  street  of  eddies  in  the  form  of  P^rs 
with  oDDOsite  rotating  directions,  almost  uniformly  aligned  along 
pddies^on  the  shelf  side  rotate  anti-cyclonically,  apparently  controlled  by  Ae  rotation, 
m  eddL  on  thi  s^  are  larger  in  the  cross-shelf  direction  than  those  on  he 
cKpif  and  thev  rotate  cvclonically.  They  become  elongated  in  the  cross-slope  dire^on 
a^fntmm^lumes^^^^^^^^  downslope  and  at  the  same  time  drag  the  surface 

water  offshore.  The  lateral  velocity  scale  is  about  the  same  at  both  the  surface  and  t  e 

^°^Figure^7^shows°lhe  velocity  fields  at  sigma  level  40  at  days  14,  18,  ^d  22.  Agm 
harnrlmic  instabilitv  results  in  a  train  of  eddies  near  the  bottom  along  the  shelf  bre^. 
The  eddies  on  the  slope  rotate  anti-cyclonically,  while  the  coupling  surface  eddies  on  the 

The  bottom  eddies  on  the  shelf  rotate  cyclonically,  and  he 
Sine  surface  eddies  on  the  shelf  rotate  anti-cyclonically.  They  are  trapped  ne^  the 
3  brl^  anparently  under  more  significant  influence  of  planetary  rotation  and  less 
sScam  iXe3  0?  bottom  friction.  The  velocity  field  near  the  bottom  a^n 
the  dense  bottom  plumes  are  dominated  by  the  friction  and  buoyancy.  This 
friction  allows  them  to  penetrate  further  downslope  before  being  held  back  and  flowing 
Sg  STSottom  plumes  also  turn  to  the  left-hand-side  of  th^r  flow  direction 
due  fo  the  Coriolis  effect  as  bottom  friction  enables  a  downslope  flow.  TJe  magnitude  of 
fcSerS  Sii?  in  the  interior  of  the  model  domain  is  on  order  of  10-20  cm/s. 

4.  A  canyon  in  Shelf  and  Slope 

To  examine  the  influence  of  a  cross-shelf  canyon  on  the  plume  dynamics  and  deep 
water  formation  we  include  a  case  with  a  100-km  wide  and  200-m  deep  canyon  tto 
leads  from  the  coast  to  deep  water,  crossing  the  shelf  and  the  slope.  Ex^mation  of  the 
cniinitv  fields  at  the  bottom.  Figure  8,  reveals  that  a  large  plume  forms  in  the  canyon  at 
-14  day^^^^^  quicW  and  spills  over  the  north  wall  of  the  canyon  as 

it  nrnorLses  downsLpe  The  large  plume  continues  to  penetrate  deeper  along  the  north 
wSl  If  the  cannon,  with  a  considerably  diffused  front  due  to  intense  mixing  and 
t  o;nm7tft  with  the  surrounding  water.  The  plumes  about  two  plume-sizes  from  the 

^  V  cold  saline  shelf  water  progressed 

Slope  mS:h  deeper  (to  about  2050  m)  along  the  north  wall  of  the  canyon  than  for 


cinnp  ra<!P  fto  about  1600  m)  at  day  22.  The  canyon  has  an  effect  similar  to 
thL^f  the  north  wall  of  the  model  in  draining  shelf  water  into  the  deep  ocean.  A  canyon, 
nottSuy  act  as  a  wall  because,!,  also  affects  .he  flows  on  .he  adjacen. 

continental  shelf  and  slope  within  a  two  plume-size  region.  ,  i  nn 

To  studv  the  effects  of  the  depth  of  a  canyon,  an  expenment  is  conducted  for  a  1(^ 

1,  anH  ?m.Tri  deen  canvon  Figure  9  shows  the  salinity  fields  in  the  presence  of  a 

m  m  deep  canyon  and  300-m  deep  canyon  at  days  22  JJ^n 

fr^nnH  that  poldef  shelf  water  penetrates  much  deeper  downslope  (to  about  2800  m)  than 
"hatot“^yr?as?  (.0  abou.  2050  m),  U.us  deep  waKr  genera..on  .s  more 

10  iUusmates  die  velociiy  fields  in  die  presence  of  die  200-.n  deep  canyon  and 
m  dMo  C^OT  fflday  22  at  level  40,  respecdvely.  The  velocity  m  die  canyon  is 
the  3(X)  m  ^  north  wall  of  the  canyon  because  more  available  potential  energy 

snengdiened  dong  not*  w^l  ^  the  Ww  canyon 

Agato  ?he  vetocity  field  lear  die  coastS  end  of  die  canyon  shows  dial  dense  shelf 
™  To'’SgTte"  Xtnlfif canyon,  we  study  a  case  with  a  ^km 

Sn  enhLes  penetration  along  it,  while  plume  penetration  ^own^ope  is  hm^^^^^^ 
flinng  the  south  wall  When  the  walls  are  far  apart,  they  act  as  a  step,  their  effects 
SJinteJ  act  when  the  canyon  width  is  close  to  order  of  the  local  deformation  radius. 

5  Conclusions 

The  numerical  simulations  show  that  baroclinic  instability  together  with  the  rotation 
detom  dirStoS^  and  bodom  plumes  along  die  shelf  break.  Corresponding^ 

Sef  UlSS’  dSuS  wiTl  Ste  of'Sie  shelf  water. 

FHftv  transnort  of  shelf  water  is  found  to  be  very  important  in  deep  water  fomation  on 
Eddy  tran  p  Th#>  nrf^cpnrp  of  a  canvon  in  the  shelf  and  slope  enhances  the 

t  “^p  A  “<•  <?“  "is 

wateJ  moduction  due  to  partial  cancellation  of  the  rotation  effect  in  the  lower  of 
tor  ?n  thp  canvon  A  deeo  canyon  has  a  more  significant  effect  in  enhancing  the  deep 
water  formation ^  draining  more  shelf  water  into  deep  ocean  through  the  canyon,  whi  e 
a  narrower  canyon  has  reduced  drainage. 
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fcQ.2.  Perspective  View  of  the  Model  Domain 
»  A  Canyon  in  the  shelf  and  slope 
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The  velocity  fields  at  sigma  level 
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The  salinity  fields  in  the  presence  of  a  200-nn  deep  and  300-nn  deep  canyon  at  sigma  level  40  at  day  22, 


The  salinity  fields  in  the  presence  of  a  300-km  wide  and  40-km  wide  canyon  at  sigma  level  40  at  day  22. 


Summary  and  Conclusions 


Ruth  H.  Preller  and  Robert  Edson 


The  purpose  of  this  workshop  was  to  familiarize  ANWAP  scientists,  both  modelers 
and  observationalists,  with  work  presently  ongoing  in  the  ANWAP.  It  was  intended  that  such 
interaction  would  foster  cooperative  work  and  a  productive  flow  of  data  between  the 
observationalists  and  the  modelers.  The  workshop  initiated  that  interaction  and  it  is  hoped 
that  such  interactions  continue  to  flourish  in  the  future. 

This  workshop  made  it  clear  that  a  large  amount  of  data  pertaining  to  processes 
important  in  assessing  the  Arctic  nuclear  waste  problem  already  exist.  Large  amount  of  data 
available  in  the  literature  as  weU  as  new  data  being  collected  by  ANWAP  scientists  need  to 
be  compiled  and  made  available  to  ANWAP  scientists.  This  task  is  presently  being  addressed 
via  development  of  a  GIS  data  base  and  the  compilation  of  an  atlas  developed  from  this  data 
to  appear  in  publication  in  1995  (Crane). 

Real  time  data,  presently  being  gathered  and  analyzed,  will  be  needed  to  understand 
the  magnitude  of  the  existing  problem  (Brooks).  Data  on  existing  source  sites  and  possible 
future  sites  need  to  be  validated,  compiled  and  made  available  to  the  community. 

Past  and  present  field  experiments  examining  sea  ice  as  a  mechanism  for  the  uptake 
and  transport  of  contaminants  via  sediments  raised  many  important  questions.  Suspension 
freezing  appears  to  be  a  key  process  in  sediment  uptake  (Reimnitz).  Is  there  any  proof  that 
contaminants  exist  in  sea  ice?  Results  of  a  recent  summer  field  experiment  presently  being 
analyzed  may  shed  some  light  on  that  question  (Tucker).  Another  important  issue  was  once 
sediment  is  incorporated  into  the  ice,  what  are  the  mechanisms  for  the  release  of  those 
sediments?  Does  sediment  remain  in  the  ice  until  it  reaches  the  ice  edge?  Or  are  there  other 
mechanisms  that  can  cause  the  release  of  sediment  from  sea  ice  into  the  water  column  before 
it  reaches  the  ice  edge?  These  were  raised  as  key  issues  for  understanding  the  importance  of 
sea  ice  as  a  transport  mechanism  for  nuclear  contaminants. 

It  was  also  pointed  out  that  there  are  many  available  historical  data  sets  dealing  with 
sea  ice.  Data  from  Arctic  buoys  exist  and  have  been  analyzed  to  reveal  prevailing  patterns  of 
both  wind  motion  and  ice  drift  for  the  Arctic  (Colony).  In  addition  to  Arctic  buoy  data,  a 
Russian  historical  data  set  (1967-1990)  of  sea  ice  conditions  now  exists  (Rigor).  This  data  set 
exists  in  Sea-ice  Data  in  Digital  Form  (SIGRID)  format  and  emphasizes  conditions  in  the 
Kara,  Laptev  and  East  Siberian  Sea. 

Not  only  is  Russian  sea  ice  data  becoming  readily  available,  but  other  types  of  data  of 
interest  to  the  ANWAP  program.  Ocean  circulation,  wind  and  tidal  data  is  available  and  has 
been  used  in  the  formation  of  numerical  models  (Proshutinsky).  In  addition,  Russian 
historical  data  is  also  becoming  available  for  Russian  rivers.  Of  particular  interest  to  the 
ANWAP  program  are  data  available  on  the  Ob,  Yenisei  and  Pechora  rivers  which  empty  into 
the  Kara  and  Barents  Seas  (Becker). 
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Another  data  set  of  great  interest  to  the  nuclear  contamination  problem  is  data  taken  in 
the  vicinity  of  the  seabed  (Sternberg)  Understanding  the  motion  near  the  ocean  bottom  is 
important  to  understanding  the  processes  of  suspension  and  settling  of  particles.  A  lengthy 
time  series  of  measurements  of  deep  currents,  temperatures  and  light  transmission  have  been 
taken  in  the  Barents/Norwegian  seas.  Preliminary  results  were  presented  and  a  detailed 
analysis  of  these  data  is  presently  underway. 

A  final  data  presentation  looked  at  the  use  of  AVHRR  data  in  studying  river  outflow 
effects  and  characteristic  albedo  and  temperature  changes  in  ice  covered  regions  near  nver 
mouths  (Fett)  Successful  past  studies  in  the  Beaufort  Sea  region  suggest  that  a  similar  study 
be  conducted  in  the  Kara  Sea.  Data  for  the  fall  freeze  up  period  of  1994  is  being  collected 
and  analyzed  to  determine  if  the  river  plume  of  the  Ob  or  Yenisei  rivers  can  be  seen  and 
described. 

Numerical  model  studies  within  the  ANWAP,  range  from  large  scale  to  regional/basin 
scale  to  small  scale  studies  represented  by  river  models.  Large  scale  models,  based  on  the 
Cox  ocean  model,  covered  the  Arctic,  its  marginal  seas  and  subarctic  marginal  seas.  These 
models  were  run  using  source  locations  mainly  in  the  Kara  Sea  with  values  derived  from  the 
Yablokov  report.  Ten  year  simulations  using  various  sources  showed  that  levels  of  radiation 
along  the  Alaskan  coast,  in  all  cases,  were  approximately  5  orders  of  magnitude  smaller  than 
those  in  the  Kara  Sea  itself.  Results  showed  that  using  Sellafield  as  the  only  source  does  not 
appear  to  account  for  radiation  levels  observed  in  the  Kara  Sea  (Preller).  Ongoing  and  future 
work  in  the  large  scale  modeling  effort  is  focusing  on  more  detailed  model  to  data 
comparison,  particularly  in  the  Kara  Sea  (Cheng).  Future  work  wUl  use  a  large  scale  coupled 
ice-ocean  model  to  study  ice  as  a  transport  mechanism  for  contaminated  sedments.  This 
work  will  have  to  be  done  in  close  cooperation  with  those  scientists  conducting  the  field  work 
on  sea  ice  as  well  as  laboratory  experiments  of  sediment  uptake  into  sea  ice,  in  order  to 
properly  parameterize  the  uptake  of  these  sediments  from  and  release  back  into  the  water 

column. 


A  complementary  study  (Rigor)  uses  ice  charts,  ice  motion  and  geostrophic  wind  to 
look  at  the  net  exchange  of  ice  area  from  the  Laptev  Sea  into  the  ice  pack.  This  study 
attempts  to  determine  the  mean  and  interannual  variability  of  ice  production  in  the  Laptev  and 
to  associate  pollution  transport  with  the  exchange  of  ice  from  this  area  into  the  pack. 

A  second  complementary  study  involves  the  development  of  a  new  ice  model,  using 
an  anisotropic  constitutive  law  for  sea  ice  (Coon).  The  new  ice  model  will  be  designed  to 
account  for  the  physical  processes  of  sea  ice  as  a  pathway  for  the  transport  of  radionuclides. 
This  new  model  could  be  adapted  by  the  large  and  regional  scale  models  if  proven  an 
improvement  over  existing  models. 

An  ongoing  cooperative  study  also  exists  between  the  river  modeling  group  and  the 
large  scale  modeling  group.  Numerous  sources,  including  the  Yablokov  report  indicate  that 
rivers  may  act  as  a  major  source  of  radioactive  contamination  of  the  Arctic.  The  river  and 
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large  scale  modelers  are  at  using  both  newly  obtained  river  data  (Becker)  and  river  model 
results  (Paluszkiewicz  and  Hibler)  as  river  source  data  for  the  large  scale  models.  Studies 
already  underway  (Allard)  using  Russian  river  data  in  the  large  scale  coupled  ice-ocean  model 
were  presented. 

Regional  models  are  also  investigating  the  importance  of  rivers  as  a  source  of  nuclear 
contamination  into  the  Arctic.  A  regional  model  of  the  Kara  Sea  is  presently  being  developed 
to  look  at  processes  which  might  contribute  to  the  redistribution  of  nuclear  contaminants  from 
rivers.  Initial  efforts  are  testing  the  responses  of  such  a  coastal  ocean  to  fresh  coastal  ocean 
buoyancy  forcing  (Smith).  Previous  work  in  the  Gulf  of  Maine  has  shown  that  multiple  river 
sources  can  interact,  merge  and  form  a  single  source  for  the  formation  of  a  coastal  current 
(Brooks).  The  Ob  and  Yenisei  rivers,  both  with  large  summer  runoff,  are  speculated  to 
interact  in  a  similar  way.  Model  studies  are  planned  to  look  at  these  interactions.  It  is 
important  to  note  that  the  large  scale  models  presently  use  grid  resolutions  too  coarse  to 
resolve  these  coastal  currents.  Regional  models  will  have  to  be  used  to  understand  these 
river/multi-river  interactions  with  the  regional  circulation. 

At  the  smaller  scales,  an  effort  to  model  the  entire  system  from  land-based  sources 
through  the  rivers  into  the  estuaries  is  presently  being  performed  (Paluszkiewicz).  The 
approach  is  to  use  a  hierarchy  of  simple  models  to  quantify  the  transport  of  contaminated 
materials  through  the  system.  Two  scenarios  are  being  evaluated  for  both  river  systems:  a 
steady  release  and  a  catastrophic  release  into  the  system.  A  study  of  a  catastrophic  release 
scenario  from  three  sites  (Krasnoyarsk,  Tomsk  and  Mayak)  was  presented  (Hibler).  Using  a 
dynamic,  one-dimensional  model,  results  showed  that  the  interaction  of  the  following  key 
parameters,  critical  shear  stress,  erodibility  and  partitioning  and  rate  constants  had  a 
pronounced  effect  on  the  contaminant  transport. 

A  somewhat  different  but  very  important  small  scale  modeling  effort  was  presented. 
This  effort  looked  at  modeling  various  release  mechanisms  into  the  environment  from  the 
radioactive  waste  source  of  the  dumped  reactor  compartment  of  the  Icebreaker  Lenin  (Mount). 
The  results  of  five  scenarios,  each  depending  on  the  integrity  of  the  levels  of  contaminant 
were  shown.  This  model  is  being  applied  to  other  forms  of  radioactive  waste  dumped  in  the 
Kara  and  Barents  Seas.  These  type  of  release  rates  may  be  used  by  the  regional  and  large 
scale  models  to  determine  more  realistic  results  than  those  obtained  using  constant  release 
rates. 


A  presentation  on  atmospheric  modeling,  another  important  factor  in  the  atmospheric- 
ice-ocean  system,  was  given  (Thompson).  An  atmospheric  model  for  the  Kara  Sea  region  has 
been  run  for  the  fall  1994  period.  This  model  uses  20  km  grid  resolution  and  is  being 
compared  to  the  model  results  of  a  corresponding  global  circulation  model  presently  run  by 
the  U.  S.  Navy.  Comparison  of  these  two  models  should  determine  whether  important 
atmospheric  effects  (wind  features)  are  lost  in  the  coarse  resolution  of  the  global  model.  The 
inability  of  the  global  models  to  calculate  these  finer  scale  features  may  influence  model 
results  of  the  transport  of  radioactive  contaminants  in  the  source  regions  (the  Kara  Sea). 
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Two  additional  modeling  presentations  (not  part  of  ANWAP)  were  presented  for 
information.  The  first  presentation  showed  some  initial  results  of  a  sea  ice-ocean  model  in 
the  Sea  of  Okhotsk  (Yang).  The  Sea  of  Okhotsk  is  also  a  source  location  for  dumped 
radioactive  waste  and  is  composed  of  a  complex  ice-ocean  system.  The  final  presentation 
showed  the  modeling  of  deep  water  formation  over  continental  slopes  (Jiang).  These 
processes  may  be  important  for  the  transportation  of  contaminants  within  the  water  column 
but  at  depth. 

The  abstracts  included  in  this  proceedings  indicate  that  a  large  effort  presently  exists 
to  investigate  the  possible  release,  dispersion  and  consequences  of  radioactive  contaminants  in 
the  Arctic.  Despite  this  breadth,  several  issues  stand  out  as  problem  areas.  Most  of  these 
areas  have  components  in  each  of  the  following  three  categories:  observational  evidence, 
process  understanding,  and  model  treatment.  We  simply  do  not  have  enough  observational 
data  to  accurately  quantify  the  important  processes,  do  not  understand  the  processes  due  to 
this  lack  of  evidence,  and  are,  therefore,  ignoring  or  poorly  treating  these  processes  in  the 
modeling  efforts. 

Principal  amongst  these  problems  is  sea  ice  as  it  relates  to  contaminant  transport.  The 
following  areas  are  of  note: 

a)  coastal  dynamics  of  ice; 

b)  methods  of  sediment  entrainment  into  the  ice; 

c)  specific  characteristics  of  the  ice  in  each  geographic  area  of  major  consequence  in 
the  contamination  problem; 

d)  characteristics  and  movement  of  the  ice  in  the  zone  between  pack  ice  and  fast  ice; 

e)  ice  contaminant  transportation  as  a  mechanism  for  contaminant  source  term 
repositioning;  and 

f)  treatment  of  ice  contaminant  transport  processes  in  large  scale  models. 

All  of  these  points  need  to  be  quantified,  understood  and,  where  necessary,  integrated  into  the 
modeling  efforts  to  accurately  quantify  contaminant  transport  in  the  Arctic. 

Of  equal  importance  to  ice  dynamics  are  the  characteristics  of  the  sediment  in  the 
Arctic  region,  and  the  Kara  Sea  in  particular.  Many  radionuclides  strongly  adhere  to 
sediments,  making  radionuclide  transport  and  sediment  dynamics  strongly  related.  The  items 
listed  below  were  highlighted  as  areas  requiring  additional  information: 

a)  sediment  transport  in  the  bottom  boundary  layer; 

b)  mid-column  sediment  transport; 

c)  dynamics  of  the  river-sediment  plume; 

d)  interactions  between  sediments  and  ice; 

e)  sediment  as  a  transport  matrix  for  contaminants;  and 

f)  treatment  of  sediment  contaminant  transport  processes  in  large  and  small  scale 
models. 
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While  many  of  these  issues  are  being  addressed  in  the  Arctic,  in  many  cases,  they  must  be 
investigated  in  the  specific  areas  of  high  contamination.  The  status  of  sediment  research  in 
the  Kara  Sea  must  be  better  understood. 

Several  additional  data/process/model  deficiencies,  while  not  specifically  noted  or 
discussed  in  detail  during  the  workshop,  should  be  briefly  listed: 

1)  data  from  periods  of  maximum  river  flow; 

2)  data  during  winter  conditions;  .  .  r 

3)  complete  annual  records  for  the  sediment,  ice  and  hydrographic  charactenstics  of 

the  regions  of  interest; 

4)  migration  patterns  of  animals  that  frequent  the  zones  of  elevated  radionuclide 
activity  on  the  Russian  land-mass  and  in  the  localized  zones  close  to  Novaya  Zemlya  with 
relatively  high  radionuclide  concentrations  in  the  marine  sediments; 

5)  modeling  pollutant  trajectories,  from  release  to  the  atmosphere  to  deposition  on 

sea-ice  or  into  ocean  waters;  and  ...  .  •  u-u- 

6)  understanding  of  geochemical  processes  in  Siberian  Rivers  that  may  inhibit 

radionuclide  transport  towards  the  ocean. 

Resolution  of  these  deficiencies  is  actively  being  pursued  through  the  ANWAP  investigations. 

Relative  to  the  modeling  efforts  as  a  whole,  several  points  are  important  and  should  be 
highlighted.  The  relevant  time  and  space  scales  of  the  contaminant  transport  problem  need  to 
be  investigated.  Currently,  a  comprehensive  approach  addressing  all  scales  has  been  adopted, 
but  this  approach  may  not  be  necessary.  To  efficiently  answer  the  question  of  contaminant 
transport  through  the  Kara  Sea  and  the  Arctic  basin  as  a  whole,  larger  regional  and  basin 
scale  modeling  may  be  all  that  is  needed.  Additional  concerns  lay  in  model  validation. 

Several  of  the  ANWAP  modeling  efforts  have  reached  the  level  of  maturity  where  validation 
is  an  important  and  necessary  next  step.  Those  modeling  efforts  still  just  beginning  should 
have  a  thorough  model  validation  effort  as  an  integral  part  of  the  research  plan.  Finally,  al 
model  efforts  should  be  initiated  with  a  through  review  of  the  observational  evidence  covenng 
the  modeled  area.  The  Kara  Sea  is  not  a  data  rich  area.  However,  recent  studies,  coupled 
with  an  increasing  amount  of  newly  released  Russian  data,  are  allowing  for  a  clearer  picture 
of  Kara  Sea  circulation  to  be  developed.  This  picture  should  not  be  ignored.  In  the  context 
of  ANWAP  goals,  Russian  data  will  continue  to  be  sought  to  fill  holes  in  our  observational 

knowledge. 

Several  generic  concerns  are  also  worth  noting.  First,  it  is  important  that  the  GIS 
database  efforts  ongoing  at  NRL  be  supported  in  a  timely  and  meaningful  manner.  This 
database  serves  as  an  archive  of  data  specific  to  this  program,  a  tool  to  assist  in  the  risk 
assessment  evaluations,  and  as  a  means  of  communication  results  to  the  international 
community.  Data  and  cruise  information  should  be  forwarded  as  quickly  as  feasible,  even  if 
it  is  caveated  to  prevent  dissemination  outside  ANWAP.  Secondly,  to  assist  in 
communication  of  the  results  and  inter-comparison  of  analysis  within  the  program,  all 
principal  investigators  must  transition  to  SI  units.  All  future  papers,  reports  and  figures 
should  be  in  SI  units. 
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A  great  deal  of  information  and  understanding  about  the  environmental  effects  of  the 
release  of  nuclear  material  in  the  Arctic  can  be  obtained  from  the  combined  efforts  of 
modelers  and  observationalists.  This  workshop  is  a  starting  point  for  these  types  of  efforts. 
It  would  be  beneficial  for  a  second  modeling  workshop  to  take  place  in  early  1996  to 
determine  the  progress  of  these  efforts. 
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